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ABSTRACT 

Results  of  measurements  from  several  Arctic  field  programs  and  numerical 
models  show  that  clouds  affect  wind  stress  during  the  central  Arctic  winter  by 
changing  the  longwave  cooling  of  the  surface  and  cloud  layers.  The  longwave 
cooling  alters  the  thermodynamic  structure  of  the  lower  atmosphere  which  in  turn 
affects  the  efficiency  of  momentum  transfer  to  the  surface.  For  typical  Arctic 
conditions,  wind  stress  is  changed  by  about  40%  one  hour  after  a  cloud  condition 
change,  due  to  changes  in  both  the  surface  layer  stability  and  surface  layer  wind 
speed.  The  actual  wind  stress  effect  due  to  clouds  during  this  time  is  a  function  of 
wind  speed,  thermal  wind,  atmospheric  boundary  layer  depth,  magnitude  of 
radiation  change,  snow  age  and,  sometimes,  snow  depth.  After  several  hours, 
surface  heat  fluxes  are  no  longer  important,  but  the  structure  of  the  atmosphere  has 
been  permanently  altered.  This  affects  the  wind  stress  by  about  10%  to  20%  during 
certain  situations,  but  can  vary  depending  on  the  initial  atmospheric  structure. 
Measurements  of  these  effects  show  variations  in  wind  stress  associated  with  clouds. 
Operational  and  research  studies  of  ice  and  ocean  dynamics  will  benefit  from 
consideration  of  cloud  effects  on  wind  stress. 
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I.  INTRODUCTION 

A.       BACKGROUND 

This  doctoral  dissertation  is  the  result  of  an  investigation  of  the  effect  of 
clouds  on  surface  wind  stress  in  the  Arctic.  The  two  physical  phenomena  which  are 
the  focus  of  this  research,  clouds  and  surface  wind  stress,  represent  very  different 
physical  processes,  but  both  play  important  roles  in  the  Arctic  marine  atmosphere 
(AMA)  and  their  interactions  are  poorly  understood. 

Clouds  affect  the  AMA  in  many  different  ways.  The  focus  of  this  research  is 
on  those  aspects  of  Arctic  clouds  which  are  most  relevant  to  modeling  and  physical 
understanding  of  surface  momentum  flux  or  wind  stress.  Wind  stress  has  special 
importance  in  the  Arctic  because  it  is  the  primary  force  driving  ice  movement 
(Thorndike  and  Colony,  1982). 

The  goal  of  this  research  is  to  answer  the  following  questions,  which  were 
originally  posed  in  the  doctoral  dissertation  research  proposal: 

(1)  Are  clouds  important  to  wind  stress  in  the  Arctic? 

(2)  What  are  the  situations  when  clouds  are  the  most  important? 

(3)  What  is  the  quantified  effect  of  clouds  on  wind  stress? 

Cloud  particles  do  not  directly  affect  wind  stress.  Latent  heat  and  radiation 
processes  within  and  outside  of  clouds  cause  diabatic  heating/cooling  of  the 
atmosphere  and  surface.  The  resulting  temperature  changes  affect  surface  heat 
fluxes  and  stability  in  the  lower  atmosphere,  which  are  crucial  factors  in  determining 
the  geostrophic  drag  coefficient  in  the  Arctic  (Overland,   1985;  Overland  and 


Davidson,  1992).  Therefore,  the  answers  to  the  above  questions  are  closely  related 
to  the  thermal  structure  of  the  atmosphere  and  snow/ice  surface. 

The  very  stable  atmospheric  boundary  layer  (ABL),  which  often  is  present 
during  clear  weather  in  the  central  Arctic,  tends  to  suppress  turbulence  and 
dynamically  de-couple  the  surface  from  the  rest  of  the  AMA.  The  surface 
temperature  is  free  to  drop  until  a  balance  is  reached  between  the  upward  and 
downward  longwave  radiation,  with  a  small  contribution  from  conductive  heat 
transfer  to  the  surface.  If  ABL  clouds  exist,  the  longwave  surface  cooling  is  mostly 
counteracted  by  radiation  from  the  cloud  base.  The  heat  originating  from  the  ice 
surface  will  be  turbulently  transported  to  the  top  of  the  clouds  where  it  will  be  lost 
from  cloud  top  longwave  radiation.  By  changing  the  location  of  the  major  heat  loss 
from  the  surface  to  the  top  of  the  ABL,  the  clouds  have  fundamentally  altered  the 
ABL. 

After  a  change  in  cloud  conditions,  the  surface  of  the  snow  is  subjected  to  an 
energy  imbalance.  This  will  cause  the  snow  surface  temperature,  hereafter  referred 
to  as  surface  temperature,  to  change  and  generate  turbulent  heat  fluxes  which  affect 
both  surface  layer  stability  and  overall  ABL  stability.  After  this  initial  phase,  which 
lasts  a  few  hours,  surface  layer  stability  is  no  longer  important,  but  clouds  will 
continue  to  affect  wind  stress  by  changing  the  overall  thermodynamic  structure  of 
the  ABL. 

B.       MOMENTUM  FLUX  AT  THE  AIR/ICE/SEA  INTERFACE 

The  transfer  of  momentum  between  the  atmosphere  and  the  ocean  directly 
affects  both  the  atmospheric  and  oceanic  boundary  layers  (OBL).  Momentum 
transfer  and  the  associated  shear  generate  turbulent  kinetic  energy  (TKE)  in  the 
ABL  and  the  OBL.  An  increase  in  TKE  can  enhance  mixing  and  entrainment,  thus 


thickening  the  ABL  and  OBL.  The  mean  currents  in  the  upper  ocean  are  also 
driven  primarily  by  the  momentum  flux  from  the  atmosphere. 

The  formation,  destruction  and  movement  of  sea  ice  is  controlled  by  heat  and 
momentum  fluxes.  Campbell  et  al  (1987)  showed  that  Fram  Strait  ice  movement 
and  deformation  were  primarily  controlled  by  the  atmospheric  momentum  flux 
when  surface  wind  speeds  are  greater  than  6  m/s.  Ocean  currents,  internal  stresses 
and  surface  tilt,  which  also  drive  ice  motion,  are  related  to  wind  stress  fields. 

Formation  of  sea  ice  is  a  direct  result  of  atmospheric  factors  such  as  high 
winds  and  low  temperatures.  Melting  also  can  be  directly  caused  by  the 
atmosphere,  but  in  the  Greenland  Sea  Marginal  Ice  Zone  (MIZ)  most  melting 
occurs  when  the  ice  is  forced  over  warmer  water  (McPhee  et  al.,  1987).  In  this  case, 
the  atmosphere  causes  ice  melting  indirectly  by  moving  the  ice. 

C.      ATMOSPHERIC  FACTORS  WHICH  AFFECT  MOMENTUM  FLUXES 

Momentum  flux  is  controlled  by  the  surface  wind  speed,  surface  topography 
and  atmospheric  stability.  The  relation  between  wind  speed  and  momentum  flux 
can  be  parameterized  by  the  equation 

=      P   Q  U2,  (1) 

where  r  is  the  surface  wind  stress,  p  is  air  density,  U  is  a  surface  wind  speed  and  Cd 
is  a  drag  coefficient  which  is  a  function  of  surface  roughness,  atmospheric  stability 
and  height  of  U.  Considerable  progress  has  been  made  in  recent  years  in 
developing  schemes  for  estimating  the  value  of  Cd. 

Although  Cd  as  a  function  of  ice  conditions  is  now  well  known,  estimates  of 
wind  stress  in  the  Arctic  are  only  as  accurate  as  the  specification  of  U.     The 


dependence  of  stress  on  the  square  of  U  makes  wind  stress  estimates  particularly 
sensitive  to  U.  Although  surface  wind  data  are  typically  not  available  in  the  Arctic, 
there  are  enough  buoys  with  pressure  measurements  to  enable  the  estimation  of  the 
geostrophic  wind  speed,  UG,  and  direction  with  as  much  accuracy  as  mid-latitude 

locations  (Moritz,  1985). 

The  scalar  relation  between  surface  wind  stress  and  geostrophic  wind  can  be 
parameterized  by  use  of  a  geostrophic  drag  coefficient,  CG,  which  is  defined  as 


u* 
CG  -  \f,  (2) 


where  u. ,  the  friction  velocity,  is  defined  as 


u*  -  It  "  •  0) 


This  is  similar  to  Cd1/2,  except  that  the  geostrophic  wind  rather  than  the 
surface  wind  is  used.  Unlike  Cd,  CG  is  not  easy  to  determine  accurately.  Many 
boundary  layer  effects  such  as  stability,  baroclinicity,  isobar  curvature,  depth  of  the 
boundary  layer  and  horizontal  changes  of  these  quantities,  as  well  as  surface 
conditions  influence  CG.  Therefore,  virtually  any  phenomenon  which  affects  the 
boundary  layer  will  also  affect  the  value  of  CG. 

D.      CLOUDS  AND  LONGWAVE  RADIATION 

One  phenomenon  which  affects  the  ABL  and  CG  is  clouds.   Clouds  affect  the 

atmosphere  in  three  ways.  (1)  They  release  or  store  latent  heat  associated  with 


water  phase  changes;  (2)  they  strongly  influence  long  and  short  wave  radiation 
processes;  and  (3)  they  may  create  hydrometeors  (liquid  or  ice  particles  large 
enough  to  be  affected  by  gravity).  In  the  Arctic,  liquid  and/or  ice  clouds  can  exist  in 
the  boundary  layer.  This  study  will  analyze  only  dark  season  situations  with  no 
shortwave  radiation  effects  considered.  Latent  heat  effects  are  minor  in  the  cold 
Arctic  winters.  Hydrometeors  and  suspended  particulate  matter  affect  radiation 
and  surface  energy  balance,  but  we  do  not  have  good  measurements  of  these 
particles.  Therefore,  the  focus  of  this  study  will  be  on  how  the  longwave  radiation 
characteristics  of  clouds  affect  wind  stress. 

Longwave  radiation  is  affected  primarily  by  (1)  water  vapor,  (2)  clouds,  (3) 
C02  and  other  trace  gases  and  (4)  aerosol.  The  amount  of  water  vapor  in  the  Arctic 

atmosphere  is  limited  by  the  cold  temperatures.  The  other  gases  are  constant, 
except  on  climatic  time  scales.  Therefore,  the  effect  of  changing  cloud  cover  on 
longwave  radiation  is  particularly  important  in  the  Arctic.  The  longwave  radiation 
emitted  from  particles  near  the  top  of  clouds  is  greater  than  the  amount  absorbed, 
particularly  when  the  air  above  the  cloud  is  clear.  Usually  the  bottom  of  a  cloud  will 
have  radiation  flux  convergence,  since  the  radiation  sources  below  the  cloud  are 
usually  warmer  than  the  cloud  base.  Both  these  effects  will  tend  to  destabilize  the 
cloud  layer  and  lead  to  larger  TKE  and  entrainment. 

The  net  effect  of  the  longwave  radiation  from  low  clouds  is  to  cool  the  lower 
atmosphere.  Cooling  of  the  lower  atmosphere  can  also  occur  indirectly  due  to 
turbulent  heat  fluxes  from  a  surface  cooled  by  radiation.  This  is  common  on  clear 
nights  over  land  or  non-moist  ice.  Open  ocean  and  moist  ice  or  snow  surfaces  have 
mixing  (open  ocean)  or  liquid/solid  phase  changes  (moist  ice  or  snow)  which 
prevent  radiation  processes  from  causing  large  temperature  changes  of  the  surface 


material.  Therefore,  in  the  ABL  of  the  AMA,  the  major  loss  of  heat  due  to 
longwave  radiation  is  caused  directly  by  clouds  and  indirectly  by  surface  cooling. 
The  latter  is  only  important  when  the  surface  is  dry  (below  freezing)  and  no  clouds 
are  present  in  the  ABL. 

This  dissertation  is  arranged  as  follows.  After  this  introduction,  the  current 
state  of  knowledge  concerning  wind  stress  and  clouds  in  the  Arctic  is  reviewed.  The 
discussion  of  results  begins  with  an  analysis  of  the  effect  of  clouds  on  the  thermal 
structure  of  the  AMA  using  experimental  data.  Next,  the  short  term  effects  of 
clouds  on  surface  layer  stability  are  analyzed.  Following,  longer  term  and  complete 
ABL  physics  are  modeled  in  order  to  determine  the  effect  of  clouds  on  wind  stress 
for  realistic  situations.  Then,  statistical  relationships  between  clouds  and  wind  stress 
are  examined.  The  conclusion  will  summarize  the  significance  of  the  results  and 
potential  applications. 


II.  PREVIOUS  RESEARCH 

This  chapter  will  examine  the  state  of  knowledge  concerning  wind  stress  and 
clouds  in  the  Arctic  previous  to  this  dissertation.  This  chapter  is  divided  into  three 
sections  (1)  Results  from  wind  stress  studies  before  the  Marginal  Ice  Zone 
Experiment  (MIZEX)  programs,  (2)  Cloud  studies  based  on  programs  before 
MIZEX  and  (3)  Wind  stress  and  cloud  studies  from  MIZEX  and  later  programs. 
Data  from  the  MIZEX  program  and  later  studies  were  obtained  by  the  investigators 
at  the  Naval  Postgraduate  School. 

A.       WIND  STRESS  MEASUREMENTS  AND  MODELS  BEFORE  MIZEX 

A  summary  of  results  from  all  wind  stress  measurements  over  sea  ice  reported 
in  Western  literature  before  1983  is  available  from  Overland  (1985).  Some  of  these 
results  will  be  mentioned,  but  the  reader  is  referred  to  this  source  for  references  to 
all  studies  of  Cd  and  summaries  of  specific  values  obtained  as  a  result  of  these 

studies. 

1.       Early  Studies 

Since  the  earliest  explorations  of  the  Arctic,  there  has  been  interest  in 
the  effect  of  wind  on  ice  movement  (Nansen,  1902;  Sverdrup,  1933;  Shuleikin,  1938). 
These  earlier  studies  estimated  magnitude  ratios  and  angles  between  wind  velocity 
and  ice  drift. 

The  recent  emphasis  has  been  on  measurement  and  modeling  of  the 
ratio  between  wind  speed  and  momentum  transfer.  The  first  estimates  of  surface 
drag  coefficients  over  sea  ice  were  based  on  measurements  of  surface  wind  speed 
profiles  from  towers  (Untersteiner  and  Badgley,  1965,  Doronin,  1969,  Ling  and 


Untersteiner,  1974).  Smith  et  al  (1970)  is  the  first  publication  reporting  Cd 
estimates  based  on  direct  eddy-correlation  measurements  using  a  sonic 
anemometer. 

2.  AIDJEX 

The  first  extensive  program  to  relate  wind  stress  to  surface  and 
atmospheric  factors  in  the  Arctic  occurred  during  the  Arctic  Ice  Dynamics  Joint 
Experiment  (AIDJEX),  which  took  place  in  the  Beaufort  Sea  in  1971-1976.  Results 
of  surface  eddy  correlation  measurements  during  AIDJEX  were  reported  by  Banke 
and  Smith  (1973)  and  Banke  et  al  (1976,  1980),  while  results  from  dissipation 
measurements  were  reported  by  Leavitt  (1980).  The  value  of  Cd  from  these  studies 

ranged  from  1.14  x  10  to  1.90  x  10  .  The  locations  for  these  measurements  were 
chosen  to  avoid  any  wake  effects  from  ice  ridges.  Therefore,  the  Cd  values  were  too 

low  to  be  used  for  area  averaged  wind  stress  estimates. 

AIDJEX  aircraft  eddy  correlation  and  momentum  integral 
measurements  reported  by  Brown  (1977),  Katz  (1979,1980)  and  Carsey  (1980) 
range  from  1.7  x  10  to  2.8  x  10  .  The  higher  values  obtained  by  the  latter  methods 
are  in  good  agreement  with  hindcast  estimates  of  Cd  from  floe  trajectories  in  Coon 

(1980),  Hibler  (1979)  and  Neralla  et  al  (1980). 

3.  Other  Studies 

In  the  late  1970's  and  early  1980's  more  surface  and  aircraft  based 
measurements  of  wind  stress  over  sea  ice  were  obtained  using  a  variety  of  methods. 
These  are  reviewed  by  Overland  (1985)  and  will  not  be  discussed  individually  here. 
Overland's  summary  of  surface  drag  coefficients  is  shown  in  Table  1.   Higher  drag 


TABLE  1 

COMPOSITE  SURFACE  DRAG  COEFFICIENTS  AS  A  FUNCTION  OF  ICE 

AND  METEOROLOGICAL  REGIME 


r, 

,  <  -5°. 

r, 

,  <  -5°. 

Ice  Regime 

Characteristics 

T„  -0° 

z, 

<  300  m 

z, 

>  400  m 

Smooth  ice 

large,  flat  floes 

1.5" 

1.5* 

Arctic  pack 

large  range  of  floe  sizes,  large 
pressure  ridges.  Ct  >0.9 

1.7' 

2.6''r 

Marginal  seas 

broken,  first-year  ice.  C,  = 
0.9.  occasional  big  floes 

2.2' 

2.7* 

.1.0'' 

Inner  MI7. 

small  floes,  rafted.  C»  = 
OR-0.9 

2.6' 

3.0' 

3.7' 

Outer  Ml/, 

C,  =  0.4 

C,  -  0.3,  nibble  field 

2.2' 

2.R'" 

Footnotes  refer  to  the  data  sets  that  form  the  basis  for  the  coefficient.  Numbers  represent  a 
subjective  median  for  a  range  of  values.  Hrror  estimates  are  not  possible. 


coefficients  were  obtained  over  sea  ice  in  the  marginal  seas  than  in  the  central 
Arctic.  The  highest  values  were  found  in  the  MIZs. 

4.       Wind  Drag  Coefficient  Models 

Arya  (1973,1975)  partitioned  the  surface  momentum  stress  into  a  skin 
drag  and  a  form  drag.  The  skin  drag  was  caused  by  the  regular  floe  surface  while 
the  form  drag  was  due  to  pressure  ridges.  The  surface  measurements  from  AIDJEX 
were  considered  to  be  measurements  of  skin  drag  only.  Guest  and  Davidson 
(1984b)  extended  this  model  to  include  the  effect  of  floe  edges,  which  are  usually 
the  dominant  roughness  elements  in  the  MIZ. 

Banke  et  al.   (1980)  related  Cd  to  surface  roughness  elements  after 

filtering  wavelengths  greater  than  13  cm.  Brown  (1981)  modeled  the  geostrophic 
drag  coefficient  using  a  two  layer  ABL  model  containing  a  surface  layer  and  a 
modified  Ekman  layer.  His  model  included  the  effect  of  secondary  circulations  and 
thermal  wind.  Above  a  classic  surface  layer,  the  wind  profile  was  calculated  by  using 
similarity  functions  which  were  empirically  determined  from  AIDJEX 
measurements.  The  value  of  CG  was  sensitive  to  stability  at  near-neutral  conditions. 

The  effect  of  clouds  was  not  considered. 

B.       CLOUD  STUDIES 

The  effect  of  clouds  on  the  ABL  has  been  the  object  of  a  considerable 
research  in  recent  years.  Stull  (1988)  provides  a  good  review  of  our  basic  knowledge 
concerning  clouds  and  the  ABL.  Almost  all  observational  cloud  studies  have 
occurred  in  the  mid-latitudes  and  tropics.  Many  of  the  observed  effects  of  clouds  on 
the  dynamics  of  the  ABL  can  be  applied  to  the  Arctic,  although  differences  in 
temperature,  cloud  condensation  nuclei  (CCN)  size  spectra  and  other  factors  must 
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be  considered.  The  most  common  type  of  clouds  in  the  Arctic  are  stratus, 
stratocumulus  and  fog  (fog  will  be  considered  to  be  a  type  of  cloud).  These  types  of 
clouds  are  the  focus  of  this  dissertation. 

1.  Early  Studies 

Vowinckel  and  Orvig  (1970)  compiled  cloud  statistics  for  the  marine 
Arctic  which  showed  a  high  percentage  of  stratus  in  the  summer.  Jayaweera  and 
Ohtake  (1973)  observed  that  Arctic  stratus  usually  occurred  in  several  layers,  unlike 
mid-latitude  stratus.  This  layering  was  modeled  by  Herman  and  Goody  (1976). 
They  show  that  the  layering  can  be  caused  by  shortwave  radiation  which  is  trapped 
inside  a  cloud  that  is  opaque  to  longwave  radiation.  They  propose  that  diurnal 
effects  prevent  similar  layering  in  mid-latitude  clouds. 

2.  AIDJEX 

The  stratus  cloud  coverage  during  AIDJEX  was  anomalously  low  in 
1975,  due  to  unusual  northerly  winds  (Jayaweera,  1977).  Jayaweera  also  describes 
results  from  May  1976  Cessna  180  aircraft  missions  which  found  that  the  inversion 
would  lower  to  the  center  of  the  stratus  clouds  after  a  few  days. 

AIDJEX  included  two  radiation  missions  by  the  NCAR  Electra  aircraft 
described  by  Herman  (1977).  He  estimated  Arctic  stratus  shortwave  radiation 
parameters  such  as  bulk  values  for  reflectance,  transmittance  and  absorbance,  and 
other  parameters  such  as  single  scattering  albedo,  absorption  optical  depth  and 
multiple  scattering  parameters. 

Using  data  from  the  same  flights,  Herman  (1980)  determined  various 
longwave  parameters  for  the  Arctic  stratus,  including  mass  absorption  coefficients 
for  selected  liquid  water  distributions.  He  concluded  that  the  longwave  radiative 
properties  of  Arctic  stratus  are  similar  to  mid-latitude  stratus. 
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3.       Arctic  Stratus  Cloud  Experiment 

During  June  1980,  six  NCAR  Electra  missions  were  flown  over  the 
Beaufort  Sea  as  part  of  the  Arctic  Stratus  Cloud  experiment  (ASC).  This  was  the 
most  comprehensive  program  to  specifically  study  Arctic  clouds.  The  first  direct 
measurements  of  cloud  liquid  water  and  droplet  spectra  were  made.  There  were  no 
surface  based  measurements  during  ASC. 

Tsay  and  Jayaweera  (1984)  described  the  cloud  morphology  and  drop 
size  spectra  in  the  clouds.  Tsay  and  Jayaweera  (1983)  emphasized  the  different 
radiative  properties  resulting  from  wide  variations  in  droplet  spectra.  Herman  and 
Curry  (1984)  analyzed  the  effect  of  the  clouds  on  shortwave  radiation  with  the  aid  of 
a  theoretical  model.  Curry  and  Herman  (1985a)  analyzed  the  longwave  radiation 
properties.  These  papers  contain  detailed  tables  and  figures  of  cloud  characteristics 
and  radiation  measurements.  There  was  a  large  variability  observed  in  the  cloud 
characteristics  on  different  days  so  that  a  "typical"  summer  Arctic  stratus  cloud  could 
not  be  described.  Large-scale  factors  which  affect  the  occurrence  of  Arctic  stratus 
were  examined  by  Curry  and  Herman  (1985b). 

Curry  (1986)  examined  the  interactions  between  turbulence,  radiation 
and  microphysics  in  Arctic  stratus  by  examining  four  case  studies  of  the  ASC  flights. 
She  has  several  conclusions  which  are  important  to  understanding  how  Arctic  clouds 
might  affect  surface  wind  stress.  (1)  The  cloud  layer  is  often  de-coupled  from  the 
surface,  and  several  cloud  and/or  fog  layers  may  exist.  Therefore  the  entire  ABL 
cannot  be  modeled  as  a  mixed  layer.  (2)  The  clouds  themselves  are  well  mixed  due 
to  cloud  top  radiative  cooling  and  latent  heat  but  not  due  to  surface  fluxes.  (3)  The 
longwave  cooling  in  the  clouds  always  exceeds  the  warming  from  shortwave 
absorption.    The  magnitude  and  vertical  distribution  of  cooling  is  sensitive  to  size 


12 


spectra  and  liquid  water  amount  at  the  top  of  the  clouds.  (4)  The  direct  effects  of 
entrainment  on  temperature  and  liquid  water  are  not  obvious  lower  than  50  meters 
below  the  cloud  top.  (5)  There  is  no  evidence  of  inhomogeneous  mixing  as 
described  by  Baker  and  Latham  (1979).  (6)  Surface  fluxes  of  moisture  do  not 
contribute  significantly  to  maintaining  the  clouds.  Gravitational  settling  has  a  strong 
effect  on  the  particle  size  spectra  in  different  regions  of  the  clouds.  (7)  The  spectral 
dispersions  of  droplet  radii  are  very  large  due  to  a  variety  of  mechanisms.  (8)  Only 
10%  of  the  cloud-top  cooling  is  balanced  by  turbulent  sensible  and  latent  fluxes 
from  below.  The  rest  of  the  heat  comes  from  entrainment  and  droplet  growth. 
Curry  does  not  mention  subsidence. 

4.       Other  Cloud  Studies 

A  major  current  program  to  study  marine  stratocumulus  clouds  is  called 
the  First  ISCCP  Regional  Experiment  (FIRE;  ISCCP  =  International  Satellite 
Cloud  Climatology  Project).  Much  of  the  research  related  to  this  program  is  based 
on  a  field  program  which  occurred  off  the  coast  of  southern  California  during  the 
summer  of  1987  (Kloessel  et  al,  1988). 

Borisenkov  et  al.  (1985)  numerically  modeled  the  influence  of  Arctic 
clouds  on  large  scale  meteorology.  They  found  Arctic  stratus  clouds  would  increase 
precipitation  in  Europe  and  other  places  far  from  the  Arctic. 

C.      WIND  STRESS  AND  CLOUDS:  MEASUREMENTS  AND  MODELS  SINCE 
MIZEX 

1.       Surface  Wind  Stress  Measurements 

Monin-Obukhov  surface  layer  similarity  theory  (reviewed  by  Dyer,  1974) 
applies  to  surfaces  that  are  horizontally  homogeneous.    But  recent  evidence  (see 
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Stull,  1988)  shows  that  the  theory  can  also  be  used  when  some  horizontal  variability 
is  present.  Therefore,  more  recent  Arctic  surface  measurements  were  obtained  with 
ridges,  ice  floe  edges,  or  other  roughness  elements  upwind,  instead  of  in  the  center 
of  large  flat  floes  as  in  earlier  studies. 

During  the  Marginal  Ice  Zone  Experiments  of  1983  and  1984,  (MIZEX- 
83  and  MIZEX-84)  many  geophysical  studies  were  carried  out  (MIZEX  Group, 
1986).  Results  from  wind  stress  measurements  were  reported  by  Guest  and 
Davidson  (1984a,  b,  c,  d,  1985,  1987a),  Davidson  and  Guest  (1986,1987),  Davidson 
and  Geernaert  (1984a,  b,  1985),  Fairall  and  Markson,  (1987)  and  Anderson,  (1987). 
All  the  studies  show  that  Cd  increased  with  ice  concentration  in  the  MIZ.  Guest  and 

Davidson  (1987a)  and  Anderson  (1987)  measured  higher  drag  coefficients  over 
rough  ice  than  had  been  previously  reported.   The  highest  MIZ  values  of  Cd,  5-6  x 

10" ,  were  measured  in  regions  of  very  rough  ice  that  had  been  broken  up  and  rafted 
due  to  swell  action. 

During  the  MIZEX-87  spring  field  program  (MIZEX  "87  Group,  1989) 
the  air  temperature  was  well  below  freezing  and  new  ice  was  forming.  Many  stages 
of  ice  formation  were  encountered.  Guest  et  al.  (1988),  Davidson  and  Guest  (1988) 
and  Guest  and  Davidson  (1991a)  determined  values  of  Cd  as  a  function  of  all  ice 

types  (Table  2). 

The  highest  values  of  Cd  ever  measured  over  sea  ice  occurred  during  the 

Coordinated  Eastern  Arctic  Experiment  (CEAREX)  drift  phase.  Multi-year  ice  had 
undergone  extensive  deformation  as  it  was  pushed,  along  with  the  vessel 
Polarbjoern,  toward  Kvitoya  island.  Pressure  ridges  as  high  as  4  m  were  common 
around  the   ship.   An  average  Cd  value  of  8  x  10"3  was  obtained  during  a  two-day 
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TABLE  2 
THE  ROUGHNESS  LENGTH,  z(),  AND  NEUTRAL  DRAG  COEFFICIENT,  Cdn 

FOR  VARIOUS  ICE  AND  SEA  SURFACES 


_  - 



— 

70  x  I01  m.  Median 

C*  *  I01 

Ice  Type 

Median 

Minimum 

Maximum 

Grease 

00027 

0.7 

0.6 

II 

Nilas 

0.45 

1.6 

1.4 

1.9 

Pancake 

Diameter  <0.75  m 

0.016 

0.9 

0.7 

13 

Diameter  0.75-1. 5  m 

045 

1.6 

II 

2.2 

Diameter  >l.5  m 

2.8 

2.4 

1.9 

2.9 

Fused 

1.0 

1.9 

1.5 

2.6 

Young 

Smooth 

24 

2.3 

1.9 

2.7 

Rough 

7.5 

3.1 

2.6 

3.6 

First  year 

Very  smooth* 

0.33 

1.5 

1.2 

1.9 

Smooth 

13 

2.0 

1.6 

2.4 

Rough 

7.5 

3  1 

2.2 

40 

Very  rough 

21.0 

4.2 

3.1 

50 

Mulliyear 

Very  smooth* 

0  33 

15 

12 

19 

Smooth 

20 

2.2 

1.9 

2.5 

Rough 

10  0 

3.4 

2.5 

4.1 

Very  rough 

27.0 

4.6 

3.6 

55 

Extremely  rough 

110  0 

8.0 

6.7 

9.1 

Ice-free  water 

(Steady  state  open  ocean)t 
U  =  10  ms' 

0  15 

13 

na 

na 

U  *  25  ms  * ' 

1.3 

2.0 

na 

na 

East  Greenland  Sea 

[V  <  12  m  «') 

All  wind  directions 

080 

1.8 

0.7 

3.0 

tee  upwind  2-10  km 

0.23 

1.4 

1.1 

1.8 

Values  are  based  on  author's  measurements  unless  otherwise  noted  and  are  only  from  periods  when 
lowest  inversion  was  higher  than  150  m.  At  least  85%  of  the  stated  ice  type  was  upwind. 
'Includes  studies  summarized  by  Overland \\985). 
tFrom  Smith  (1988).  No  ranges  given. 
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period.    These  recent  measurements  show  that  Cd  can  have  higher  values  than 
reported  by  Overland  (1985)  in  certain  regions. 

As  a  result  of  these  studies,  the  value  of  the  surface  drag  coefficient  can 
be  accurately  estimated  for  variations  in  ice  condition.  For  this  to  be  useful, 
methods  must  be  developed  to  identify  types  and  concentrations  of  ice  in  the  Arctic. 
The  ice  maps  currently  produced  by  NOAA  give  only  general  descriptions  of  ice 
types  and  concentrations. 

2.       Observations  of  Clouds  During  MIZEX  and  CEAREX 

During  the  MIZEX  or  CEAREX  programs  continuous  surface 
measurements,  observations  and  rawinsonde  profiles  provide  a  large  data  base  to 
study  many  aspects  of  clouds  in  the  Arctic  (Davidson  et  ai,  1984;  Lindsay,  1985; 
Guest  and  Davidson,  1988).  During  certain  periods,  there  were  SODAR 
measurements,  aerosol  measurements  and  aircraft  missions  with  cloud  physics 
instrumentation. 

Guest  (1985)  developed  techniques  for  forecasting  fog  in  the  MIZ. 
Stratus  or  low  fog  is  virtually  always  present  during  on-ice  winds  in  the  MIZ  and  fog 
forms  from  stratus  lowering.  With  parallel  wind  flow,  boundary  layer  fronts  are 
common,  and  the  onset  of  fog  or  stratus  will  be  abrupt. 

Statistical  studies  by  Guest  and  Davidson  (1987b,  1988)  show  that  the 
presence  and  thickness  of  clouds  were  correlated  with  the  height  of  the  inversion 
base,  Zj.   When  stratus  clouds  were  present,  Zj  was  higher  and  was  less  correlated 

with  wind  speed  than  clear  sky  cases  (Figure  1). 

Guest  et  at  (1988)  compared  median  Zj  values  for  on-ice  and  off-ice 

winds  in  the  MIZ  (Figure  2).   At  all  relative  locations,  Zj  medians  were  higher  for 

on-ice  winds  than  off-ice  winds.    The  difference  in  Zj  becomes  greatest  over  the 
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Figure  1  Composite  Median  Inversion  Base  Heights,  Z„  from  MIZEX-84  Divided 
into  Wind  Speed  and  Cloud  Thickness  Categories.  Note  that  during  clear  (no 
cloud)  conditions  the  inversion  was  lower  and  more  greatly  affected  by  wind  speed 
than  during  cloudy  conditions. 
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Figure  2  Composite  Median  Zj  Values  Across  the  MIZ  for  Different  Wind 
Directions  during  MlZEX-87.  Triangles,  circles  and  squares  represent  off-ice  (left- 
to-right),  parallel-right  (out-of-paper)  and  on-ice  wind  conditions,  respectively. 
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pack  ice.  During  on-ice  winds,  fog  or  stratus  inevitably  would  exist  over  the  ice. 
The  destabilizing  effect  of  the  clouds  may  be  a  major  reason  for  the  observed 
difference  in  Z;  between  wind  regimes.  (There  are  also  influences  due  to  advection 
and  subsidence.)  If  the  clouds  affect  Z-,  they  will  also  affect  surface  momentum  flux, 
because  a  higher  Zj  is  associated  with  more  efficient  transport  of  momentum  to  the 
surface. 

3.       Boundary  Layer  Models 

Although  surface  drag  coefficients  are  well  known  and  straightforward  to 
use,  their  use  requires  measurements  of  surface  wind  speed.  Estimating  wind  stress 
from  a  pressure  field  is  much  more  complicated.  This  involves  the  specification  of 
the  geostrophic  drag  coefficient,  CG,  which  requires  some  type  of  ABL  model. 

Overland  (1985)  used  a  steady-state  one-dimensional  ABL  model  to 
determine  CG  and  the  turning  angle  between  the  geostrophic  and  surface  wind,  a,  as 

a  function  of  a  mechanical  mixing  stability  parameter.  He  found  that  when  the 
atmosphere  became  very  stable,  the  surface  wind  tends  to  de-couple  from  the  upper 
winds,  and  CG  is  no  longer  affected  by  surface  roughness.     These  results  are 

applicable  in  the  Arctic  pack  ice  away  from  the  MIZs.  Clouds  effects  were  not 
considered. 

In  the  1980's,  the  focus  of  Arctic  studies  was  on  MIZs.  The  few  ABL 
models  for  the  MIZ  that  have  been  published  will  be  reviewed.  Overland  et  at 
(1983)  used  a  primitive  equation  slab  model  to  explain  observations  of  the  Bering 
sea  MIZ  in  March  during  off-ice  winds.  Reynolds  (1984)  examined  the  same  case 
with  a  similar  model,  but  included  radiative  cooling  from  clouds.  He  attributed  the 
observed  variations  in  temperature,  wind  and  ABL  height  to  drag  variation  over  the 
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ice  and  heat  flux  over  the  open  ocean.  Clouds  and  heat  flux  through  the  ice  were 
not  important. 

Andreas  et  al  (1984)  described  a  MIZ  ABL  based  on  a  rawinsonde 
cross-section  of  the  Antarctic  MIZ  during  off-ice  winds  in  the  Weddell  sea.  Bennett 
and  Hunkins  (1986)  simulated  this  case  with  a  two-dimensional,  multi-level  model 
which  includes  radiation.  They  concluded  that  adiabatic  lifting  due  to  convergence 
and  longwave  cooling  from  clouds  dominate  the  cooling  in  the  ABL.  Although 
there  was  some  discussion  concerning  model  details  (Andreas,  1987;  Bennett  and 
Hunkins,  1987),  all  authors  agreed  that  there  were  large  changes  in  the  wind  stress 
across  the  MIZ  due  to  roughness  and  cloud  effects. 

Chu  (1986a,  b,  c,  1988a,  b,  c)  has  modeled  several  mechanisms  by  which 
air-ice-sea  interactions  may  cause  ice  features  in  the  MIZ.  Chu  et  al  (1990)  and 
Chu  and  Garwood  (1990,  1991)  investigated  feedback  mechanisms  between  clouds 
and  surface  fluxes  over  ice-free  oceans. 

Brown  (1986)  applied  his  one-dimensional,  two-level  ABL  model  to  an 
off-ice  wind  case  in  the  Fram  Strait  MIZ.  A  realistic  surface  roughness  and 
temperature  field  was  used.  The  greater  roughness  in  the  outer  MIZ  slows  and 
backs  the  surface  wind.  A  greater  upward  heat  flux  at  the  surface  has  the  opposite 
effect. 

Recent  results  of  two  dimensional  multi-level  MIZ  ABL  models  by 
Glendening  (1992)  and  Kantha  and  Mellor  (1989)  show  the  same  general  results  for 
off-ice  winds.  Kantha  and  Mellor  also  examined  the  MIZ  ABL  during  on-ice, 
parallel-left  and  parallel-right  wind  regimes.  The  above  studies  show  the  surface 
stress  field  is  very  complicated  in  the  MIZ  because  of  the  changes  in  the  ABL  due  to 


20 


surface  roughness  and  horizontal  temperature  variations.  Cloud  effects  have  not 
been  modeled  by  these  researchers. 

Glendening's  model  of  the  MIZ  ABL  during  off-ice  winds  shows  that  Zj 

is  controlled  by  surface  wind  stress  over  the  ice  regions  and  surface  temperatures 
(heat  fluxes)  over  open  ocean  regions.  He  has  not  examined  very  stable  cases  when 
surface  conditions  become  de-coupled  from  the  upper  ABL. 

Modeling  stable  ABL  situations  is  more  difficult  because  of  the 
intermittent  nature  of  the  mixing  events.  Overland  (1988)  approaches  the  problem 
with  a  one-dimensional,  multi-level,  TKE-mixing  length  model  which  simulates  the 
ABL  over  ice  in  the  winter.  When  cloud  radiation  effects  are  introduced,  the 
boundary  layer  becomes  unstable  and  surface  stress  increases  by  approximately 
20%. 

4.       Wind  Stress  in  the  MIZ 

The  above  models,  as  well  as  observations  (Davidson  and  Guest,  1988) 
show  that  the  surface  wind  stress  field  does  not  match  the  surface  drag  coefficient 
field  because  there  are  changes  in  the  surface  wind  speed  across  the  MIZ.  The  wind 
speed  at  the  surface  is  affected  by  changes  in  stability  and  surface  roughness.  By 
using  observations  and  Overland's  (1985),  Brown  and  Liu's  (1982)  and  Brown's 
(1986)  ABL  models,  Campbell  et  al.  (1987)  and  Guest  (1988)  showed  that  the 
combined  effects  of  horizontally  varying  surface  wind  speeds  and  roughness  result  in 
a  wind  stress  field  which  is  complicated  and  dependent  on  the  wind  direction 
relative  to  the  ice.  Typically,  the  wind  stress  is  greater  over  the  open  ocean  than  the 
pack  ice,  even  though  Cd  has  an  opposite  distribution.  This  fact  has  been  neglected 

by  several  MIZ  ice  movement  models  (Roed,  1983;  Roed  and  O'Brien,  1983; 
Hakkinen,  1986a,b;  Smith  et  al.  1988).    These  models  use  an  unrealistic  constant 
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surface  wind  speed  and  direction  across  the  MIZ,  which  results  in  greater  stress  over 
the  ice  regions. 

Guest  (1988)  also  shows  that  clouds  can  have  a  large  effect  on  surface 
wind  stress  if  they  change  the  stability  of  the  surface  layer  from  stable  to  unstable. 
This  effect  of  clouds  will  be  further  examined. 
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III.    AN  EXAMINATION  OF  THE  RELATION  BETWEEN  CLOUDS  AND 
ATMOSPHERIC  BOUNDARY  LAYER  THERMODYNAMIC  STRUCTURE 

Clouds  influence  wind  stress  by  their  thermodynamic  effects  on  the 
atmosphere.  This  chapter  will  examine  the  thermodynamic  effects  of  clouds,  which 
form  the  basis  for  understanding  the  effect  of  clouds  on  wind  stress. 

Divergence  of  longwave  radiation  at  cloud  tops  and/or  at  the  snow  surface  has 
a  significant  effect  on  the  thermodynamic  structure  of  the  lower  atmosphere  over 
pack  ice  during  the  Arctic  winter.  Overland  and  Guest  (1991),  hereafter  OG, 
showed  how  the  temperature  at  the  surface  was  primarily  controlled  by  longwave 
radiation.  Heat  conduction  through  the  ice  and  leads  are  less  important  but  do 
prevent  extremely  cold  temperatures  ( <  -40  C)  from  existing  over  sea  ice.  The  layer 
of  air  just  above  the  inversion,  or  isothermal  layer,  contains  the  warmest 
temperatures  in  the  atmosphere  (Figure  3).  It  will  tend  to  experience  more 
radiational  cooling  than  the  rest  of  the  atmosphere.  Warm  horizontal  advection 
must  counteract  the  longwave  cooling  in  the  isothermal  layer,  on  the  average. 

The  actual  location  of  the  maximum  longwave  cooling  in  the  lower 
atmosphere  depends  on  the  cloud  conditions.  The  location  of  the  cooling  affects  the 
stability  and  therefore  the  dynamics  of  the  lower  AMA.  When  conditions  are  clear, 
the  maximum  longwave  cooling  occurs  in  a  thin  layer  at  the  top  of  the  snow/ice. 
When  clouds  exist,  the  maximum  radiational  cooling  occurs  at  the  top  of  the  upper 
cloud  layer.  There  is  a  strong  correlation  between  amount  of  cloud  cover  and 
downward  longwave  radiation  (Figure  4).  When  low  overcast  is  present,  all  of  our 
Arctic  data  indicate  a  clear  trend:  the  downward  longwave  radiation  at  the  surface  is 
within  ±  10  Wm"2  of  the  blackbody  radiation  at  cloud  bottom.  Therefore,  it  can  be 
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Figure  3  A  Typical  Temperature  Sounding  From  the  CEAREX  Drift.  The  mixed 
layer,  lower  200  m,  inversion  layer,  200  to  500  m,  and  the  isothermal  layer  (500  to 
1200  m)  are  shown. 
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Figure  4     Frequency  Distribution  of  Downward  Longwave  Radiation  during  the 
CEAREX  Drift  for  Clear  (small  squares)  and  Low  Overcast  (large  diamonds) 
Conditions.    Note  that  the  peak  frequency  for  cloudy  conditions  is  about  75  Wm" 
higher  than  for  clear  conditions. 
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assumed  that  all  but  the  thin  (less  than  100  m)  low  level  cloud  layers  that  exist  in  the 
Arctic  are  essentially  black  to  longwave  radiation.  The  examination  will  be  based 
on  data  obtained  during  the  MIZEX  and  CEAREX  programs. 

A.       CEAREX  DRIFT  MEASUREMENTS 

Most  of  the  observational  data  used  for  this  study  were  from  measurements 
from  the  Polarbjoern  during  the  CEAREX  drift  phase,  15  September  to  12 
December,  1988.  During  the  drift  phase,  the  Polarbjoern  was  moored  to  a  large  ice 
floe  which  drifted  in  the  region  North  and  East  of  Svalbard.  Instrumentation  was 
located  both  on  a  large  ice  floe  (alpha  floe)  to  which  the  Polarbjoern  was  moored, 
and  on  board. 

A  meteorological  station  was  located  on  a  platform  extending  forward  of  the 
Polarbjoern's  bow  mast  at  a  height  of  14  meters  above  sea  level.  This  station 
measured  wind  speed  and  direction,  temperature,  and  relative  humidity. 
Temperature  data  were  accurate  to  within  1.0  C.  The  humidity  measurement  was 
not  reliable.  Wind  speeds  were  accurate  to  within  0.3  m/s  and  direction  to  within  10 
degrees,  although  during  periods  of  unfavorable  wind  direction,  errors  may  have 
been  larger. 

The  radiation  system  measured  both  downward  infrared  and  downward  solar 
irradiance,  and  was  located  on  the  ice  until  the  floe  was  crushed  on  November  15. 
At  this  time  the  system  was  relocated  aboard  the  Polarbjoern.  The  system  remained 
operational  until  December  12.  The  sensors  would  often  frost  quickly;  therefore 
hourly  cleanings  were  performed.  Radio  transmissions  would  interfere  with  the 
measurements.  The  data  were  manually  edited  to  eliminate  poor  measurement 
periods.      The   instruments  were  modified  so  that   errors   due   to   temperature 
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differences  within  the  sensor  could  be  identified  and  corrected.  The  resulting 
irradiances  were  generally  accurate  to  within  5  watts/m  . 

Rawinsondes  provided  vertical  profiles  of  temperature,  humidity,  and  wind 
speed  and  direction  twice  daily  throughout  the  experiment,  with  additional 
soundings  during  periods  of  extreme  or  unusual  weather.  The  sondes  were 
equipped  with  thermistors  that  measure  temperatures  to  within  0.2  C.  The  humidity 
sensor  consisted  of  a  specially  coated  glass  plate,  the  resistance  across  which  varies 
with  humidity.  It  was  accurate  to  about  5%.  Sonde  wind  directions  were  generally 
accurate  to  within  20  degrees. 

A  6  meter  profile  mast,  located  on  an  adjacent  floe,  measured  temperature 
and  wind  speed  at  4  levels.  This  was  used  to  provide  wind  stress  and  heat  flux 
measuirements  when  the  wind  direction  was  favorable.  The  heat  flux  was  estimated 
to  be  accurate  to  5  Wrrr2  and  the  wind  stress  to  20%.  Surface  temperature  was 
measured  with  a  thermistor  placed  on  the  snow  surface  and  was  accurate  to  1.0  C. 

B.       THERMODYNAMIC  COUPLING  OF  THE  SURFACE,  CLOUD  LAYER 
AND  ISOTHERMAL  LAYER 

This  section  will  use  the  results  of  measurements  from  the  CEAREX  drift  to 

examine  the  relationship  between  surface  temperature,   air  temperature,  ABL 

temperature  and  isothermal  layer  temperature  over  ice  in  the  Arctic  winter.    (In  a 

strict  sense,  the  term  "air  temperature"  refers  to  the  temperature  at  10  m  above  the 

surface.  In  practice,  the  temperature  at  the  top  of  the  ice  profile  mast,  5.2  m,  and  at 

the  ship  mast,  14  m,  were  used  to  determine  air  temperature.)    The  differences 

between  these  temperatures  determine  the  stratification  of  the  lower  atmosphere 

which  in  turn  affects  momentum  transfer  to  the  surface.  The  effect  of  clouds  on  the 

stratification  will  be  examined. 
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The  surface  temperature  and  the  air  temperature  are  closely  linked  (Figure  5) 
because  surface  sensible  heat  flux  quickly  counteracts  any  temperature  difference 
between  the  snow  and  near-surface  air.  The  one  period  when  there  was  a  difference 
of  greater  than  2  C,  2-6  November  1988,  occurred  during  very  low  or  zero  winds, 
when  turbulent  surface  fluxes  were  suppressed.  Usually,  enough  turbulence  exists  to 
support  surface  layer  heat  fluxes  and  therefore  the  difference  between  surface 
temperature  and  ABL  potential  temperature  is  usually  very  small. 

OG  showed  how  surface  temperature,  and  therefore  air  temperature,  is 
thermally  coupled  to  a  radiational  boundary  layer  (RBL).  The  RBL  is  characterized 
by  an  isothermal  temperature  layer  above  an  inversion  layer.  The  isothermal  layer 
extends  to  approximately  1.5  km  elevation.  The  surface  is  prevented  from  becoming 
extremely  cold  because  the  heat  that  is  lost  by  upward  longwave  radiation  is 
replaced  by  downward  radiation  from  the  isothermal  layer,  plus  a  small  amount  of 
heat  conduction  upward  through  the  ice/snow. 

OG  did  not  consider  the  effect  of  clouds,  which,  in  the  isothermal  layer  or 
ABL,  would  greatly  increase  the  downward  radiation  at  the  surface  compared  to 
clear  sky  conditions.  This  causes  a  very  close  thermodynamic  coupling  between  the 
clouds  and  the  surface.  This  coupling  is  illustrated  in  Figure  6,  a  time  series  plot  of 
surface  air  temperature  vs.  sky  temperature.  Sky  temperature,  Tsky,  is  defined  as  the 

temperature  of  a  blackbody  corresponding  to  the  measured  downward  radiation  at 
the  surface,  Rdownsfc 


i 
Tsk,    *     {-^T^f,  (4) 


downsfc  *^ 
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Figure  5   Time  Series  of  Surface  Temperature,  5.2  Meter  Air  Temperature,  Wind 
Speed  and  Air-Snow  Temperature  Difference  during  the  CEAREX  Drift. 
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Figure  6     Time  Series  of   14   Meter  Air  Temperature   (dashed   line)   and   Sky 
Temperature  (solid  line)  during  the  CEAREX  Drift. 
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where  cr  =  5.67  x  10  Wm  K  is  the  Stefan-Boltzmann  constant.  The  air 
temperature,  Tair,  in  Figure  6,  is  usually  very  close  to  the  sky  temperature,  except  for 

a  few  periods  of  colder  sky  temperatures. 

The  periods  when  Tair  was  close  to  Tsky  were  cloudy  periods,  while  the  periods 
when  Tsky  was  colder  occurred  during  clear  periods.  This  is  illustrated  by 
scatterplots  of  Tair.  vs.  Tsky  for  overcast  and  clear  sky  conditions  (Figures  7  and  8).  A 
link  between  Tair  and  Tsky  exists  for  all  sky  conditions.  This  is  because  the  dry  snow 
surface  has  a  low  heat  capacity  and  low  thermal  conductivity  and  therefore  the  snow 
surface  temperature  quickly  responds  to  changes  in  surface  energy  balance  caused 
by  changes  in  downward  radiation.  The  snow  surface  is  then  closely  linked  to  the 
ABL  air  temperature  by  turbulent  heat  flux  as  discussed  above. 

During  low  level  overcast  periods,  a  better  correlation  between  air 
temperature  and  surface  temperature  exists  (Figure  8).  Tair  and  Tsky  are  never  more 

than  a  few  degrees  different  from  each  other.  During  these  overcast  periods,  the 
surface  is  not  only  linked  by  radiation  with  the  surface,  it  is  also  turbulently  linked  if 
the  cloud  bottom  is  within  the  ABL.  ABL  turbulence  effectively  transfers  heat 
between  the  surface  and  the  top  of  the  ABL  during  overcast  conditions. 

A  conclusion  of  OG  was  that  longwave  radiational  cooling  of  the  snow  surface 
causes  the  lower  atmosphere  to  cool  during  the  dark  seasons.  The  temperature  of 
the  air  near  the  surface  is  determined  by  a  balance  between  upward  radiation  at  the 
snow  surface  and  downward  radiation  from  the  isothermal  layer  (with  a  small  effect 
due  to  heating  from  the  ocean  through  leads  and  ice).  This  study  proposes  that 
these  conclusions  be  modified  for  situations  when  the  ABL  is  cloud-capped.  During 
these  periods,  important  radiational  exchanges  are  occurring  at  the  cloud  top  rather 
than  at  the  surface. 
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Figure  7  Scatterplot  of  Sky  Temperature  vs.  14  Meter  Air  Temperature  for  Clear 
Sky  Conditions  during  the  CEARBX  Drift.  The  line  represents  surface  temperature 
equal  to  sky  temperature. 
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Figure  8  Same  as  Figure  7  but  for  Low  Overcast  Conditions. 
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Overland  and  Davidson  (1992),  hereafter  referred  to  as  OD,  found  that  wind 
stress  was  affected  by  a  measure  of  the  external  or  background  atmospheric 
stratification  based  on  the  difference  in  temperature  between  the  surface  and  an 
upper  level  (900  or  850  mb),  assumed  to  be  in  the  isothermal  layer.  When  ABL 
clouds  exist,  the  surface  is  no  longer  directly  connected  by  radiation  to  the  upper 
level.  In  this  case,  an  external  stability  parameter  defined  by  the  difference  in 
temperature  between  the  cloud  top  (instead  of  the  surface)  and  an  upper  level  may 
be  appropriate,  since  it  is  at  the  cloud  top  that  the  direct  radiational  linking  between 
the  surface  and  the  upper  levels  is  actually  occurring. 

C.       CLOUDS  AND  BOUNDARY  LAYER  THERMAL  STRUCTURE 

To  determine  the  effect  of  clouds  on  wind  stress,  it  is  important  to  understand 
how  clouds  affect  or  are  associated  with  ABL  structure.  Two  parameters  for 
describing  ABL  structure  were  derived  from  hundreds  of  different  MIZEX  and 
CEAREX  rawinsonde  soundings.  These  were  (1)  height  of  the  lowest  inversion 
base,  Zj,  and  (2)  "type"  of  boundary  layer,  described  below. 

1.       Inversion  Height  -  Average  Values 

Low  cloud  cases  have  a  higher  average  Z,  than  clear  sky  cases.   This  is 

shown  in  Table  3  for  locations  within  the  pack  ice,  in  the  MIZ  and  in  the  adjacent 
open  ocean.  Average  inversion  heights  are  considerably  higher  for  cloudy  skies 
compared  to  clear  skies.  This  is  true  for  the  various  locations  shown  in  Table  3  and 
also  is  true  for  any  data  set  grouped  by  season,  platform,  wind  direction  regime  (off- 
ice, on-ice)  or  wind  speed  regime  (not  shown). 
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TABLE  3 
MEAN  AND  MEDIAN  Z;  FOR  DIFFERENT  LOW  CLOUD  CONDITIONS 

AND  LOCATIONS 


Location 


Cloud  Condition 


Pack  ice 


1 


MIZ^ 

Open  Ocean 

400 

610 

480 

730 

390 

470 

60 

420 

220 

450 

Overcast 


Clear 


median 
mean 

sd4 

median 
mean 

sd4 


220 
270 

210 

100 
160 

240 


370 


340 


i 
50  km  or  more  from  ice  edge  (primarily  CEAREX  drift  data) 

2 
<50  km  from  ice  edge  but  still  over  sea  ice 

3 
over  open  ocean  but  within  200  km  of  ice  edge 

4  •    • 

sd  =  standard  deviation  These  values  represent  natural  variations,  not 

experimental  errors. 
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The  large  standard  deviations  of  the  Zj  values  show  that  other  factors 

besides  clouds  affect  ABL  structure.  An  analysis  of  variance  shows  that  22%  of  the 
variation  in  Z;  values  not  associated  with  distance  from  the  ice  edge  can  be 
attributed  to  cloud  conditions.  The  higher  average  Z;  values  associated  with  low 
clouds  are  not  due  entirely  to  the  direct  effect  of  clouds  on  the  ABL.  Cloud 
occurrence  is  associated  with  warm  advection,  low  pressure  systems,  high  wind 
speeds  and  moist  marine  air  masses.  All  of  these  would  be  associated  with  higher  Zj 

values  even  without  cloud  effects. 

In  this  section,  Zj  vs  cloud  relationships  were  examined,  and  in  the 

previous  section,  the  coupling  between  longwave  radiation  and  ABL  air 
temperature  for  different  cloud  conditions  was  discussed.  Results  from  both  topics 
can  be  combined  by  plotting  the  difference  between  Tsky  and  Tair  vs.  Zj  for  clear  and 

low  overcast  conditions  (Figure  9).  Several  differences  between  overcast  and  clear 
conditions  are  apparent.  Overcast  conditions  are  characterized  by  Tsky  being  nearly 
equal  or  slightly  greater  than  Tajr,  and  by  having  higher  average  inversion  bases 
when  compared  with  clear  cases.  Surface-based  inversions  (Zj  =  0)  are  extremely 
rare  during  cloudy  skies,  but  clear  skies  do  not  guarantee  that  a  surface  inversion 
will  be  present. 

2.       Inversion  Height  and  u.  Scaling 

It  has  been  observed  since  early  explorations  (Sverdrup,  1933)  that  the 
value  of  Z,  is  related  to  surface  wind  speed.    Increasing  wind  speed  generates  more 

turbulent  kinetic  energy,  TKE,  due  to  turbulent  shear  production,  which  is  able  to 
extend  mixing  to  higher  levels.  Attempts  to  scale  or  predict  Zj  or  the  ABL  depth,  h, 

for  stable  or  neutral  ABLs  invariably  use  parameters  related  to  the  mechanical 
production  of  turbulence:  the  friction  velocity,  u.,  an  eddy  viscosity  coefficient,  K, 
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Figure  9  Scatterplot  of  the  Difference  Between  the  14  Meter  Air  Temperature  and 
the  Sky  Temperature  vs.  Inversion  Base  Height  for  Clear  (squares)  and  Low 
Overcast  (plusses)  Sky  Conditions. 
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the  bulk  Richardson  number,  RiB,  or  a  wind  speed  or  wind  shear  parameter.  The 
following  scales  for  equilibrium  ABL  height,  he,  have  been  proposed  for  stable  and 
neutral  ABLs. 

Neutral  ABL  and  Free  Atmosphere: 


2K\  2 
he   =   n  { —  )  Ekman  depth 


u 


he  =  0.2  \   f  j  e.g.  Panofsky  and  Dutton  ( 1984) 

Stable  Surface  Layer,  Neutral  Free  Atmosphere: 

hp   =   L  =  =  Obukhov  length  scale 

kg  w'  ev' 

/u*L\j 
he     a    \-J~J  Zilitinkevich  (1972,  1974) 

RiBUh2  6 

K     a    ~7n — 7~\  Hanna  (1969),  Wetzel  (1982) 

Neutral  Surface  Layer,  Stable  Free  Atmosphere: 

LI  jk  r\  a 

he    oc    ^";  N2    =    f~^  OG,Kitaigorodskii(1988), 

Kitaigorodskii  and  Joffre  (1988) 

The  parameter  vr  was  shown  by  OD  to  be  relevant  for  scaling  the  Arctic 

winter  boundary  layer,  when  surface  fluxes  are  usually  small  and  an  external  or 
background  stability  limits  the  height  of  the  ABL.     OD  used  data  from  the 


38 


CEAREX  drift  to  verify  this  relationship.  A  re-analysis  of  these  data  confirms  that 
Zj  was  correlated  (R2  =  0.34)  with  u.  or  wind  speed  during  the  CEAREX  drift  data 

(Figure  10).  A  similar  correlation  (R2  =  0.21)  occurs  for  the  data  set  containing  all 
MIZEX  and  CEAREX  soundings  made  by  our  group  in  the  Arctic  (Figure  11)  when 
all  sky  condition  cases  are  included.  However,  when  either  data  set  is  divided  into 
categories  depending  on  cloud  conditions,  it  is  found  that  the  wind  speed  (or  u.) 

dependence  completely  disappears  (R2  <  0.05)  for  the  low  overcast  cases  (Figures  12 
and  13). 

The  wind  dependence  noted  by  OD  and  others  apparently  occurs  only 
during  clear  sky  periods.  This  lack  of  wind  speed  dependence  for  Z;  during  overcast 

conditions  was  observed  in  the  MIZ  during  the  summer  in  MIZEX-84  (Guest  and 
Davidson,  1987b)  and  in  the  spring  during  MIZEX-87  (Guest  et  al,  1988)  as  well  as 
over  pack  ice  during  the  fall/winter  CEAREX  drift  (Guest  and  Davidson,  1991a). 
Therefore,  low  clouds  appear  to  reduce  the  relative  influences  of  surface  layer 
mechanical  production  of  turbulence  on  the  Zj  values. 

The  latter  conclusion  seems  to  contradict  some  basic  assumptions  about 
the  factors  that  control  ABL  depth  in  a  stable  atmosphere.  The  lack  of  a  wind 
speed  vs.  Z,  correlation  could  have  several  explanations.   (1)  Zi?  as  defined  here,  is 

not  always  a  good  measure  of  ABL  height,  h.  For  the  modeling  studies,  the  ABL 
height,  h,  is  arbitrarily  defined  as  the  height  at  which  TKE  becomes  10%  of  the 
surface  value.  Clouds  may  create  a  sharp  inversion  which  defines  the  Zj  value  when 

mixing  is  actually  limited  to  a  lower  level  which  cannot  be  identified  from  the 
rawinsonde  data.  (2)  The  boundary  layer  under  clouds  in  the  Arctic  is  usually 
unstable  and  therefore  buoyant  production  of  TKE  may  dominate  over  mechanical 
production  at  times.    (3)  Cloud  top  cooling  and  wind  shear,  both  of  which  often 
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Figure  10  A  plot  of  Inversion  Height  vs.  Wind  Speed  during  the  CEAREX  Drift. 
The  solid  line  is  the  linear  regression  fit,  the  inner  dashed  lines  are  the  95% 
confidence  intervals  of  the  regression  and  the  outer  dashed  lines  are  the  95%  data 
range  limits. 


40 


All  MI2EX  and  CEAREX 


£ 
01 
•H 

■ 
I 

■ 
■ 
II 

CD 

c 

0 
•H 

* 

L 

B 

C 


1200 


1000 


800 


600    - 


400 


200 


0    - 


Figure  11    Same  as  Figure  10  but  using  All  Dntci  from  MIZEX-84,  MIZEX-87  and 
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Figure  12  Same  as  Figure  10  but  for  Low  Overcast  Cases  Only. 
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Figure  13  Same  as  Figure  11  but  for  Low  Overcast  Cases  Only. 
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occur  at  the  top  of  the  ABL,  may  be  more  important  than  surface  processes  in 
controlling  Zj  during  overcast  conditions.  (4)  Advection  of  moisture  at  certain  levels 
may  fix  Z;  at  that  level  regardless  of  surface  conditions.  (5)  Cloud  top  entrainment 
instability,  CTEI,  (Lilly,  1968)  is  another  process  not  directly  related  to  surface  wind 
stress.  A  necessary  condition  for  CTEI  is  that  the  equivalent  potential  temperature 
of  the  air  in  ABL  must  be  greater  than  the  air  just  above  ABL.  This  condition  is 
rare  in  the  Arctic;  therefore  CTEI  cannot  be  common. 

Although  this  sub-section  concerned  the  effect  of  wind  stress  on  an  ABL 
structure  parameter,  Zj?  the  focus  of  this  study  remains  the  more  subtle  reverse 

process:  the  effect  of  cloud-influenced  ABL  structure  on  wind  stress.  These  results 
have  illustrated  how  the  surface  wind  stress  influence  on  boundary  layer  structure  is 
strongly  dependent  on  cloud  conditions.  Scaling  inversion  height  with  u.  is  not  valid 

during  low  overcast  conditions,  which  provides  further  evidence  that  clouds  have  a 
strong  effect  on  ABL  structure. 

3.       Types  of  Atmospheric  Boundary  Layers 

Another  way  to  describe  ABL  thermal  structure  is  to  classify  each 
measured  ABL  according  to  a  scheme  summarized  in  Table  4,  where  the  number 
assigned  each  type  is  arbitrary.    For  types  0,  3,  and  5,  the  specification  of  Z,  is 

straightforward  since  these  represent  "classic"  ABLs  where  a  clearly  defined 
inversion  exists  at  the  surface  or  above  a  surface-based  well-mixed  layer.  For  the 
other  types  the  Zj  specification  is  more  ambiguous  because  either  no  inversion  exists 
or  weakly  stable  layers  exist  below  a  more  obvious  inversion.  In  these  cases,  Zj  is 
still  defined  as  the  lowest  inversion  base,  but  probably  is  not  a  good  measure  of 
ABL  depth,  h. 
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TABLE  4 
TYPES  OF  ATMOSPHERIC  BOUNDARY  LAYERS 


Type 
-1  No  inversion  in  lower  2500  m 

0  Strong  surface-based  inversion  (height  of  top  of  inversion  ,Zt     >  200 
m  or  potential  temperature  increase  within  the  inversion,  0jump  >  3  C) 

1  Weak  sfc-based  inversion  below  mixed  layer  (Zt0    <  200  m  and  Ztop  > 
100  m  and  0Jump  <  3  C) 

2  Elevated  inversion  above  stable  layer;  stable  layer  not  inverted. 

3  Elevated  inversion  above  mixed  layer 

4  Miscellaneous  category  (complicated  structure) 

5  Two  inversions  with  mixed  layer  between 

0.5  Combination  of  0  and  5;  A  strong  surface-based  inversion  below  a 

well-mixed  layer  capped  by  an  upper  level  inversion. 
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All  rawinsonde  profiles  we  have  collected  in  the  Arctic  during  MIZEX 
and  CEAREX  (over  1300)  have  been  classified  according  to  ABL  type  and  amount 
of  low  clouds.  The  percentage  of  total  of  profiles  for  each  type  is  listed  in  Table  5. 

As  indicated  previously  (Figure  9),  surface-based  inversions  are  rare 
during  low  overcast  conditions  and  common  during  clear  skies.  All  the  low  level 
inversions  occurred  over  ice  and  during  clear  skies,  with  the  exception  of  four  cloudy 
cases  when  strong  warm  advection  over  cold  ice  in  the  MIZ  caused  a  surface-based 
inversion  to  form.  During  the  CEAREX  drift,  there  were  no  surface-based 
inversions,  Type  0,  with  overcast  skies.  This  is  a  significant  observation  with  respect 
to  the  wind  stress  because  during  these  very  stable,  low  inversion  periods 
atmospheric  stability  determines  how  wind  stress  is  related  to  geostrophic  forcing 
(Overland,  1985;  and  later  in  this  dissertation).  During  less  stable,  higher  inversion 
periods,  changes  in  surface  roughness  are  more  important  than  atmospheric  stability 
in  affecting  wind  stress.  Clouds  play  an  important  role  in  destroying  or  preventing 
the  formation  of  strong  surface-based  inversions. 

Multiple  upper-level  inversions  (Type  5)  were  twice  as  common  for 
overcast  vs.  clear  skies.  Most  of  the  Type  5  clear  cases  had  high  relative  humidities 
at  the  bottom  of  the  upper  inversions,  which  may  indicate  evidence  of  previous 
clouds.  Arctic  clouds  can  form  mixed  layers  capped  by  inversions  which  are  de- 
coupled from  the  ABL. 

When  surface-based  inversions  existed  below  an  upper-level  inversion 
(Type  0.5,  12  cases),  only  one  case  had  clouds  associated  with  the  upper  inversion. 
Clouds  prevent  surface  inversions  from  forming  even  if  they  are  de-coupled  from 
the  surface. 
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TABLE  5 


A  CROSS  TABULATION  OF  BOUNDARY  LAYER  TYPE  BY  LOW-LEVEL 
CLOUD  CONDITIONS  EXPRESSED  AS  PERCENTAGE  OF  ALL  CEAREX 

DRIFT  SOUNDINGS1 


ABL  Type' 


Low  Cloud  Amount 
Clear  Partial  Overcast 


Surfaced-based  inversion  13.3 

Type  0  or  0.5 

Inversion  above  mixed  layer  34.8 

Type  3  or  Type  5 

Two  inversions  23.7 

Type  5  or  0.5 

Unclear  inversion  location  4.4 

Type  1,  2  or  4 


4.4 


11.1 


12.6 


2.2 


0.7 


26.7 


17.0 


2.3 


Based  on  135  rawinsonde  profiles.  Total  of  all  categories  is  153%.  This  is  because 
14%  of  the  soundings  were  Type  0.5,  which  are  both  surface-based  and  two 
inversions  and  39%  of  the  soundings  were  Type  5  which  is  both  mixed  surface  layer 
and  two  inversions. 

Type  numbers  are  defined  in  Table  4. 
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D.      CLOUDS  AND  SURFACE  FLUXES 

One  way  clouds  can  affect  surface  wind  stress  is  by  their  influence  on  the 
surface  layer  stratification.  This  section  is  an  examination  of  how  observed  sky 
conditions  and  sky  condition  changes  are  related  to  measured  surface  sensible  heat 
fluxes.  Turbulent  heat  flux  at  the  surface  determines  surface  layer  stability  which 
affects  the  value  of  the  surface  drag  coefficient  by  the  well-verified  Monin-Obhukov 
similarity  theory  (Dyer,  1974).  Also,  surface  heat  fluxes  can  affect  wind  stress  by 
influencing  entrainment  into  the  ABL.  Entrainment  brings  momentum 
characteristic  of  upper-level  air  into  the  ABL,  causing  a  short  term  effect  on  surface 
wind  stress.  Entrainment  also  changes  the  ABL  depth,  h,  which  can  have  a  long 
term  effect  on  wind  stress.  Results  from  numerical  models  of  these  processes  are 
presented  in  Chapter  V. 

Thirty-one  periods  from  the  CEAREX  drift  were  identified  when  the  wind 
direction  was  favorable  for  surface  heat  flux  measurements,  meteorological 
conditions  remained  constant  and  several  consecutive  ten  minute  heat  flux  values 
were  similar.  Low  cloud  cover  was  correlated  with  surface  sensible  heat  flux  (Figure 

-2  -2 

14).  The  average  heat  flux  was  -5  Wm  "  and  1  Wm  during  clear  skies  and  low 
overcast  skies  respectively.  A  better  correlation  occurred  between  measured 
sensible  heat  flux  and  downward  longwave  radiation  at  the  surface  (Figure  15)  or 
net  surface  radiation  (Figure  16),  although  scatter  remains.  The  zero-flux  calm  wind 
cases  are  not  shown.  All  the  measured  heat  fluxes  from  the  CEAREX  drift  were 
small  compared  to  typical  radiative  terms  and  MIZ  or  lead  turbulent  heat  fluxes. 
However,  the  largest  changes  in  the  effect  of  stability  on  wind  stress  occurred  in  the 
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Figure  14  Surface  Sensible  Heat  Flux  (positive  upward)  as  a  Function  of  Cloud 
Cover  Octaves  from  Selected  CEAREX  Drift  Periods.  Lines  represents  statistics 
described  in  Figure  9  caption. 
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Figure   15     Surface  Sensible  Heat  Flux  as  a   Function  of  Downward  Longwave 
Radiation  (Wm").  Very  low  winds  cases  are  not  shown. 
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near-neutral  regime  (e.g.  Brown  and  Liu,  1982),  so  even  small  differences  in  heat 
flux  and  surface  stability  may  significantly  affect  wind  stress. 

Most  periods  represented  by  these  data  had  quite  constant  meteorological 
conditions,  including  cloud  cover.  Larger  sensible  heat  fluxes  were  measured 
immediately  after  changes  in  cloud  conditions.  This  is  best  illustrated  in  a  case 
study  (Figure  17)  of  a  period  from  the  CEAREX  drift.  Note  the  close  correlation 
between  cloud  conditions  and  surface  sensible  heat  flux  for  this  case.  Cloud  clearing 
is  closely  linked  to  a  change  from  unstable  to  stable  surface  conditions  while  the 
reverse  occurs  for  cloud  forming.  Several  other  similar  examples  exist  from  the 
CEAREX  drift  period. 

Sixteen  sharp  cloud  change  events  (eight  totally  clear  to  low  overcast  changes 
within  an  hour,  eight  low  overcast  to  clear  changes)  were  identified  during  periods 
when  surface  flux  measurements  were  available.  The  surface  layer  of  the 
atmosphere  changed  from  unstable  to  stable  in  all  the  clearing  cases.  The  surface 
layer  changed  from  stable  to  unstable  in  five  of  the  "clouding"  cases.  The  other 
three  clouding  cases  showed  no  significant  change  in  heat  fluxes  or  stability.  It  is 
obvious  that  surface  fluxes  are  closely  related  to  changes  in  cloud  conditions, 
particularly  immediately  after  clearing  events. 

E.   FACTORS  INFLUENCING  SURFACE  HEAT  FLUX 

Instantaneous  and  integrated  turbulent  surface  heat  fluxes  affect  the  ABL 
structure  and,  therefore,  the  ability  of  the  atmosphere  to  transfer  momentum  to  the 
surface.  It  has  been  shown  that  changes  in  cloud  conditions  are  often  associated 
with  changes  in  the  surface  turbulent  heat  flux.  But  other  factors  can  also  induce 
heat  flux  changes.    In  order  to  evaluate  the  significance  of  the  effect  of  clouds  on 
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wind  stress  in  the  central  Arctic  winter,  the  effect  of  downward  longwave  radiation, 
i.e.  cloud  conditions,  on  surface  heat  fluxes  must  be  compared  with  other  processes 
that  can  generate  a  surface  heat  flux.  This  section  describes  a  method  for  classifying 
a  non-neutral  ABL  according  to  the  mechanism  driving  the  surface  turbulent  heat 
flux.  The  classification  scheme  will  then  be  applied  to  the  CEAREX  drift  data  to 
evaluate  the  relative  importance  of  changes  in  surface  radiation  conditions  in 
generating  surface  heat  fluxes  relative  to  other  mechanisms. 

1.       General  Classification  Scheme 

With  no  outside  forcing  or  advection,  the  turbulent  surface  heat  flux 
drives  the  Tair  -  Tsfc  difference  toward  zero  over  time.  As  shown  in  section  A,  these 

two  parameters  are  closely  linked  unless  no  turbulence  exists  in  the  surface  layer. 
Surface  heat  flux  counteracts  some  "outside"  forcing  mechanism  creating  a 
difference  between  Tair  and  Tsfc.  There  are  two  types  of  outside  forcing 
mechanisms:  Type  R  (radiational)  ones  that  act  directly  to  change  Tsfc  and  Type  A 
(advective)  ones  that  act  directly  to  change  Tair.  In  Type  R  ABLs,  the  surface  is 
controlling  the  temperature  in  the  ABL  through  surface  fluxes;  changes  in  Tsfc 
precede  changes  in  Tair.  In  Type  A  ABLs,  the  temperature  in  the  ABL  is  controlling 
Tsfc  through  surface  fluxes;  changes  in  Tair  precede  changes  in  Tsfc.     These  two 

mechanisms  can  be  distinguished  by  comparing  the  signs  of  the  surface  heat  flux  and 

<5Tair 
_.  .   For  Type  R  they  are  the  same  sign  while  for  Type  A  they  are  opposite  signs. 

This  concept  is  summarized  in  Table  6. 
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TABLE  6 
TYPES  OF  TURBULENT  HEAT  FLUX  EVENTS 

Type  R  -  Radiational 

Downward  (upward)  turbulent  surface  heat  flux 
Cooling  (warming)  ABL  and  surface 

Type  A  -  Advective 

Downward  (upward)  turbulent  surface  heat  flux 
Warming  (cooling)  ABL  and  surface 
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Type  R  is  identified  as  "radiational"  because  changes  in  Tsfc  not  caused 

by  sensible  heat  fluxes  are  closely  linked  to  changes  in  radiation  conditions.  The 
only  other  factor  affecting  Tsfc  is  conductive  heat  flux  through  the  snow/ice.    This 

can  be  considered  a  passive  process  because  snow  and  ice  conditions  do  not  change 
much  on  the  time  scales  of  a  few  hours  being  considered  here.  Type  R  ABLs  are 
common  over  land  where  diurnal  variations  in  temperature  are  large,  particularly  in 
the  summer  when  solar  radiation  is  intense  and  the  soil  is  dry.  Because  advection  of 
soil  or  snow  is  insignificant,  pure  Type  R  ABLs  are  by  definition  caused  by  a  one- 
dimensional  process. 

Changes  in  radiation  conditions  in  the  Arctic  winter  are  strongly  related 
to  cloud  conditions.  These  radiational  changes  are  an  order  of  magnitude  less  than 
diurnal  solar  radiation  changes  at  lower  latitudes,  where  Type  R  ABLs  are  common. 
However,  it  will  be  shown  that  the  winter  Arctic  ABL  often  has  Type  R 
characteristics.  This  is  because  the  temperatures  of  the  dry  snow  surface  and 
shallow  ABL  quickly  adjust  to  radiation  changes. 

Type  A  ABLs  are  caused  by  any  process,  other  than  turbulent  surface 
heat  flux,  that  changes  the  temperature  within  the  ABL.  Note  that  this  is  termed 
"advective"  in  Table  6.  This  is  because  advection  is  assumed  to  be  the  dominant 
non-radiative  process  controlling  the  structure  of  the  ABL.  Type  A  ABLs  occur 
over  open  water  regions  with  SST  gradients  and  anywhere  following  an  atmospheric 
frontal  passage. 

2.       Scale  Analysis 

To  evaluate  the  factors  which  change  ABL  temperature,  consider  the 
following  equation  for  the  time  rate  of  change  of  the  vertically-integrated  potential 
temperature  of  the  ABL. 
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The  overbar  represents  a  vertical  integral  through  the  ABL: 
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time  derivative 

advection  (u  =  vector  wind  in  ABL) 

surface  flux  ( <  >  represents  time  average) 

entrainment  flux 

radiation  up  from  surface 

radiation  into  surface 

radiation  up  from  ABL 

radiation  down  into  ABL 

latent  heat  of  evaporation 

latent  heat  of  sublimation 
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Measurements  or  estimates  from  the  CEAREX  drift,  during  periods 
when  reliable  data  were  available,  were  used  to  derive  ranges  of  the  values  for  each 
of  the  terms  in  Equation  (5)  and  are  given  in  Table  7.  The  time  derivative  (Term 
1),  advection  (Term  2)  and  cloud  top  cooling  (Terms  7  and  8)  dominate  the  mean 
temperature  equation.  Surface  fluxes  (Term  3)  and  entrainment  (Term  4)  are  an 
order  of  magnitude  smaller  and  clear  air  radiation  divergence  (Terms  5,  6,  7  and  8) 
and  latent  heat  (Terms  9  and  10)  are  almost  another  two  orders  of  magnitude 
smaller.  The  smallness  of  the  latent  heat  term  justifies  ignoring  latent  heat  when 
considering  the  effect  of  clouds  on  ABL  temperature.  This  is  based  on  periods 
when  ABL  temperatures  are  below  -20  C.  As  shown  by  Curry  (1986),  the  latent 
heat  term  can  be  important  in  the  summer  when  ABL  temperatures  are  near  0  C. 

3.       Graphical  Representation  of  ABL  Temperature  and  Surface  Heat  Flux 
Events. 

Only  terms  1  and  3  in  Equation  (5)  were  directly  measured;  the  others 

were  estimated  for  typical  Arctic  conditions.  A  plot  of  the  measured  terms,  Term  1 

vs.    Term    3,    provides    a    visualization    of   factors    affecting    ABL    temperature 

(diagrammed  in  Figure  18).  The  vertical  axis  represents  surface  turbulent  heat  flux 

while   the   horizontal   axis   represents  time   change   in  ABL  temperature.     This 

represents   a   heat   flux-temperature    change    or   F-C   diagram.      A   steady-state 

equilibrium  (no  turbulent  surface  heat  flux)  situation  is  represented  by  the  origin  in 

Figure  18.    The  dashed  lines  in  Figure  18  represent  lines  of  constant  Term  2  + 

Term  4  through  9.    These  lines  will  be  termed  "isoadvects"  because  advection  is 

usually  the  dominant  term,  although  factors  other  than  advection  may  be  important 

for  individual  cases. 
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TABLE  7 
SCALE  ANALYSIS  OF  ABL  TEMPERATURE  EQUATION 

(units  are  Wm"2) 


Term 


Percentile 

5  95 


Median  Magnitude 


1. 

-121 

86 

2. 

-110 

100 

3. 

-12 

11 

4. 

0 

10 

5.  +  6.  +  7.  +  8. 

(Clear) 

-1 

1 

5.  +  6.  +  7.  +  8. 

(ABL  Cloud) 

-90 

0 

9.+ 10. 

-3 

3 

34 

35 

4 

3 

<1 

35 
1 


59 


o 

E 
3 


m 

m 

x 


£1 
-( 
tfl 
C 

ai 

(0 


lee 

I  '   '  '  '  i  y '  '   ■ 

, , , , ,  i^, , , , 

,     ,     ,      ,     y      ,      ,     , 

,D 

,B      : 

75 

/ 

/             . 

y 

Radiational 

>^ 

50 

:         / 

s 

25 

/.I 

Aduectiue 

P   ;: 

\A 

0 

y'                     ^ 

-25 

-/ 

i 

Adject  lue/ 

/ 

/ 

/ 

-50 

—  ■ 

/ 

Radiational 

/               : 

S                         i    ; 

-75 

/                  >- 

r  4 

\'P 

/ 
/ 

/ 

/ 

/ 

/ 

y           \    - 

100 

-/ 

S 

s 

A               \  ~ 

/ 1 1 1 1 1 1 1 1 1 1 

lii     i  ^    I    i     iii 

^^.^A    ,    ,    ,    , 

i  ,   ,   1-r,   i   ,   ,   ,   .  i 

100     -75      -50      -25       0       25       50       75      100 
Normalized  Temperature  Change  (Ulm~2) 


a  e 


at  pcv 


Figure  18  F-C  Diagram  Showing  Characteristics  of  Different  Types  of  Sensible 
Heat  Flux  Events.  The  X-Axis  represents  the  local  time  derivative  of  temperature 
times  ABL  depth  and  volumetric  heat  capacity  (Wm  2)  and  the  Y-axis  represents 
sensible  heat  flux  (Wm2).  Pure  Type  R  events  fall  along  line  AB  while  pure  Type 
A  events  fall  along  line  segments  CD  or  EF. 
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Pure  Type  R  situations  (i.e.  Terms  1  =  Term  3)  would  have  points  only 
along  the  zero  isoadvect,  which  is  represented  by  line  AB  in  Figure  18.  A  cloud- 
clearing  event  occurring  in  a  situation  which  was  previously  in  equilibrium  would 
cause  the  representative  location  in  Figure  18  to  move  from  the  origin  to  Point  A  as 
the  negative  surface  fluxes  respond  to  the  radiation  change.  Eventually,  the  cooling 
of  the  ABL  will  counteract  the  fluxes  caused  by  the  new  radiation  conditions  and  the 
representative  location  will  move  back  to  the  origin.  Similarly,  a  cloud-forming 
event  will  move  the  representative  location  toward  Point  B  and  then  back  to  the 
origin.  Note  that  for  pure  Type  R  situations,  the  cloud-clearing  or  cloud  forming 
must  occur  above  the  ABL;  clouds  within  the  ABL  can  affect  ABL  temperature 
directly,  a  Type  A  process. 

A  pure  Type  A  situation  (surface  temperature  driven  entirely  by  the 
turbulent  heat  flux)  would  begin  at  a  point  on  the  X-axis  because  at  first  the  ABL 
temperature  will  not  have  changed  and  the  fluxes  will  be  zero.  For  example,  at  the 
onset  of  steady  cold  advection,  the  representative  location  would  be  on  the  X-axis  at 
the  value  of  the  advection,  such  as  Point  C.  Eventually  the  cooling  ABL  will 
generate  positive  fluxes  and  the  isoadvect  will  be  followed  until  Point  D  is  reached. 
At  this  time,  the  ABL  has  cooled  enough  so  that  the  positive  heat  flux  is  balanced  by 
the  advection  and  a  steady-state  is  reached.  Alternatively,  a  warm  advection  event 
would  cause  the  representative  location  to  move  from  Point  E  to  Point  F.  Note  that 
line  AB  is  entirely  within  the  upper  right  or  lower  left  quadrants,  where  the  signs  of 
Term  1  and  Term  3  are  the  same.  Lines'CD  and  EF  are  in  the  upper  left  and  lower 
right  quadrants  where  the  signs  of  Term  1  and  Term  3  are  opposite.  This  is 
consistent  with  the  earlier,  more  general  definitions  of  Type  R  and  Type  A 
situations  which  were  defined  in  Table  6.   "Pure"  Type  R  and  "pure"  Type  A  events 
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were  shown  in  Figure  18.  In  nature,  completely  "pure"  situations  almost  never  occur 
and  advection  is  not  constant.  The  more  general  definitions  for  Type  A  or  Type  R 
categories  will  be  used  for  this  study.  Therefore  any  situation  with  changing 
temperatures  and  surface  heat  fluxes  must  be  either  Type  A  or  Type  R. 

F.       EXAMPLES  FROM  CEAREX  DRIFT 

Term  1  vs.  Term  3  from  Equation  (5)  for  the  entire  CEAREX  drift  period 
when  surface  flux  measurements  were  available  is  plotted  in  Figure  19.  The  data 
are  averaged  over  one  hour  intervals.  The  vertical  scale  in  Figure  19  is  exaggerated 
compared  to  Figure  18  so  that  the  isoadvects  are  almost  vertical.  The  majority  of 
the  points  are  clustered  toward  the  center  where  fluxes  and  temperature  changes 
are  small.  There  is  too  much  noise  (randomness)  in  the  data  to  detect  any 
discernable  patterns  in  the  central  data  cluster. 

The  large  heat  flux  events  (high  y-axis  magnitude  points  in  Figure  19)  are 
examined  more  closely  because  these  have  the  potential  to  strongly  influence  the 
effect  of  atmospheric  stability  on  wind  stress.  Unlike  the  center  of  the  scatterplot, 
the  outer  points  in  Figure  19  show  some  definite  trends.  The  largest  negative  heat 
fluxes  (below  -15  Win'")  usually  occur  during  periods  of  dropping  temperatures. 
This  is  consistent  with  a  cloud  clearing  event  which  cools  the  ABL  by  first  cooling 
the  surface.  Type  R  events  are  more  common  than  Type  A  events  for  explaining 
large  negative  heat  fluxes.  This  study  will  concentrate  on  Type  R  events,  since  these 
are  the  situations  when  clouds  have  the  strongest  effects  on  wind  stress. 

A  good  example  of  a  Type  R  event  occurred  from  1800  8  October  to  1200  9 
October  (Figure  20).  At  the  beginning  of  this  period,  the  representative  locations 
on  the  F-C  diagram  were  near  the  origin,  but  then  they  moved  down  and  to  the  left. 
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Figure  19    F-C  Diagram  of  Entire  CEAREX  Drift.    Heat  flux  scale  exaggerated 
compared  to  Figure  16. 
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Figure  20  F-C  Diagram  of  Clearing  Event  From  CEAREX  Drift  1800  8  Oct  to  1200 
9  October.  The  axes  scaling  is  different  compared  to  Figure  18. 
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Longwave  radiation  and  sky  observations  (not  shown)  clearly  indicate  that  the 
observed  decrease  in  temperature  and  in  negative  sensible  heat  flux  was  associated 
with  a  change  in  cloud  conditions.  This  is  a  Type  R  case  since  the  points  are  near 
the  zero  isoadvect. 

Most  points  in  the  lower  left  quadrant  of  Figure  19  are  to  the  left  of  the  zero 
advect.  This  means  that  the  ABL  cooling  cannot  be  entirely  explained  by  a  one- 
dimensional,  Type  R,  surface  heat  flux  argument  for  most  cases.  The  one- 
dimensional  reasoning  would  require  that  cloud  clearing  occur  everywhere  at  once. 
This  does  not  happen;  the  cloud  edge  moves  horizontally,  usually  in  the  same 
direction  as  the  air  parcels  in  the  ABL.  Therefore,  the  cooling  that  occurs  after 
cloud  clearing  is  more  intense  because  the  air  parcels  have  been  exposed  to  the 
clear  conditions  and  cold  surface  for  a  longer  period  than  the  time  after  the  cloud 
clearing  event  at  any  particular  location  on  the  surface. 

Large  upward  (positive)  sensible  heat  fluxes  are  equally  divided  between  Type 
A  and  Type  R  events  (top  of  Figure  19  to  the  left  and  right  respectively).  Cloud- 
forming  events  do  not  seem  to  dominate  the  positive  upward  heat  flux  cases  to  the 
extent  that  cloud  clearing  dominates  negative  heat  flux  events. 

The  largest  temperature  increases  (to  the  right  of  150  Wm  in  Figure  19)  are 
all  associated  with  negative  fluxes,  a  Type  A  situation.  But  large  temperature 
increases  (to  the  left  of  -150  Wm"2)  can  occur  with  both  upward  and  downward  heat 
fluxes.  Again,  this  indicates  that  cloud-clearing  events  are  a  relatively  more 
important  cause  for  surface  heat  fluxes  than  cloud-forming  events.  This  means  that 
the  magnitude  of  the  ABL  warming  associated  with  cloud  formation  is  less  than  the 
magnitude  of  the  ABL  cooling  associated  with  cloud  clearing. 
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Two  reasons  explain  the  difference  in  surface  heat  flux  response  between 
cloud-forming  and  cloud-clearing  situations.  The  first  reason  is  that  cloud-forming 
events  often  involve  ABL  clouds.  The  radiational  warming  is  counteracted  by  cloud 
top  cooling  so  that  there  are  large  positive  fluxes  at  the  surface  but  the  ABL 
temperature  does  not  warm.  Cloud-clearing  has  no  counteracting  radiational 
warming  in  the  ABL.  The  second  reason  is  that  the  depth  of  the  ABL  is  usually 
lower  for  cloud-clearing  than  cloud-forming  events.  Shallow  ABLs  respond  more 
quickly  to  surface  heat  fluxes  than  do  deep  ABLs. 

Several  conclusions  concerning  the  effect  of  clouds  on  ABL  processes  can  be 
derived  from  the  CEAREX  drift  data  discussed  in  this  sub-section.  Type  R  ABL 
events  occurred  at  about  the  same  rate  as  type  A  events  when  all  data  are 
considered.  Therefore,  cloud  condition  changes  appear  to  be  as  important  as  other 
factors  in  affecting  surface  heat  fluxes  when  all  cases  are  considered.  Large 
negative  heat  fluxes  are  particularly  well-correlated  with  cloud  condition  changes. 

This  section  demonstrates  that  clouds  are  often  important  to  surface  heat 
fluxes.  Therefore,  much  of  the  wind  stress  variation  caused  by  surface  heat  fluxes 
over  sea  ice  or  ABL  stability  effects  can  be  attributed  to  cloud  condition  changes. 

G.      CONCLUSIONS  ON  CLOUD  EFFECTS  ON  ATMOSPHERIC 
THERMODYNAMIC  STRUCTURE 

This  examination  shows  the  strong  association  between  cloud  conditions  and 

the  thermodynamic  structure  of  the  lower  atmosphere  over  sea  ice  in  winter. 

Observed  data  show  that  clouds  have  a  large  effect  on  surface  temperature  which  in 

turn  is  closely  linked  to  the  ABL  properties  by  turbulent  fluxes.  Clouds  prevent  the 

formation  of  surface-based  inversions  and  are  associated  with  higher  inversion  bases 

and  multiple-level  inversions.    The  presence  of  low  clouds  leads  to  no  correlation 
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between  the  inversion  height  and  wind  speed.  During  the  CEAREX  drift,  clouds 
had  as  great  an  effect  on  surface  layer  stability  as  all  other  factors  combined.  Cloud 
clearing  events  were  especially  effective  at  generating  relatively  large  negative 
surface  heat  fluxes  and  cooling  the  ABL. 
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IV.  THE  EFFECT  OF  CLOUDS  ON  SURFACE  LAYER  STABILITY,  THE 
SNOW/ICE  LAYER  AND  WIND  STRESS 


A.       DEFINING  SURFACE  LAYER  STABILITY  EFFECTS 

The  quadratic  geostrophic  drag  coefficient,  CG  ,  parameterizes  the  effect  of 
surface  and  ABL  physics  on  wind  stress.  The  physical  processes  which  affect  CG~ 
can  be  divided  into  two  categories  (1)  surface  layer  effects  and  (2)  outer  ABL  layer 
effects.  The  quadratic  geostrophic  drag  coefficient  can  be  expressed  as 


c2     =       (- 


G 


Q    CLG  .  (6) 


where  U10  is  the  wind  speed  10  meters  above  the  surface,  which  is  near  the  top  of 

the  surface  layer.  Surface  layer  effects  on  wind  stress  are  parameterized  by  the 
surface  drag  coefficient,  Cd,  while  outer  ABL  processes  are  parameterized  by  the 
quadratic  reduction  factor,  CLG. 

Surface  heat  fluxes  cause  surface  layer  stability  effects  which  change  the  value 
of  Cd.  The  changes  in  Cd  caused  by  clouds  will  be  analyzed  in  this  chapter.  The  next 

chapter  will  include  the  effects  of  clouds  on  the  reduction  factor. 

B.       PURPOSE  OF  SURFACE  LAYER  STABILITY  STUDIES 

Discussion  in  Chapter  III  addressed  how  changes  in  cloud  conditions  caused 
heat  fluxes  and  stratification  at  the  surface.     Here,  these  heat  fluxes  will  be 
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analytically  and  numerically  modeled  assuming  pure  Type  R  events,  i.e.  no 
advection,  entrainment  or  diabatic  ABL  cooling  other  than  surface  fluxes.  One 
purpose  of  this  section  is  to  determine  analytical  expressions  for  important  time 
scales  related  to  Type  R  heat  flux  events.  Another  purpose  is  to  determine  the  most 
important  parameters  for  predicting  heat  flux  and  wind  stress  on  various  time  scales 
over  sea  ice. 

OG  modeled  the  effect  of  an  instantaneous  change  in  downward  radiation  on 
Tair  and  surface  heat  flux  (Figure  21)  to  describe  changes  in  conditions  above  the 

ABL  from  cloudy  to  clear  and  back  to  cloudy.  They  did  not  consider  the  wind  stress 
implications  of  the  radiational  changes.  Heat  flux  values  from  Figure  21  indicated 
that  there  were  periods  (the  first  two  hours  of  day  2)  when  turbulence  and  wind 
stress  would  be  almost  completely  suppressed  and  other  periods  (first  hour  of  day  6) 
when  it  would  be  enhanced  by  30%.  These  conclusions  were  based  on  surface  layer 
Monin-Obukhov  similarity  theory  (MO).  Inclusion  of  ABL  changes  would  magnify 
the  surface  layer  stability  effects.  Although  the  OG  case  (Figure  21)  had  low  wind 
speeds  and  therefore  may  not  be  significant  to  wind  stress  variations,  it  shows  that 
for  certain  conditions,  changes  in  cloud  cover  can  have  large  relative  effects  on  wind 
stress  for  short  time  periods. 

C.       CONDUCTION  OF  HEAT  THROUGH  SNOW/ICE  (ANALYTICAL 
MODEL) 

A  skin  depth,  dskin,  is  defined  as  the  effective  or  e-folding  distance  that 

longwave  radiation  penetrates  the  snow/ice  surface.  If  there  is  an  instantaneous 
change  in  downward  surface  radiation,  Frad,  from  a  previously  equilibrium  (no  net 
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Figure  21  Model  Simulation  Time  Series  of  (a)  Surface  Temperature,  (b) 
Difference  between  Surface  and  10  meter  Air  Temperature,  (c)  Surface  Sensible 
Heat  Flux.  A  change  from  cloudy  to  clear  occurs  at  2  days  and  becomes  overcast 
again  after  6  days.  From  Overland  and  Guest  (1991). 
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flux)  situation,  then  the  surface  temperature,  Tsfc,  cannot  change  temperature  any 
faster  than 


at 


^skin^-skin    >  (') 


where  Csldn  is  the  volumetric  heat  capacity  of  the  snow  or  ice  at  the  surface. 

The  heat  conduction  within  the  snow  or  ice  cover  is  governed  by  the  one- 
dimensional  heat  flux  equation  (Semtner,  1976) 

cTx    -    KZ>  (8> 

where  k  is  the  thermal  conductivity  of  the  snow  or  ice  and  C  is  the  volumetric  heat 
capacity.  In  the  first  few  seconds  following  a  change  in  radiation  conditions,  the 
surface  temperature  has  not  changed  enough  to  generate  sensible  turbulent  fluxes. 
Therefore,  all  the  radiational  imbalance  not  stored  in  the  skin  layer  must  be 
counteracted  by  a  conductive  flux,  Fc,  from  the  skin  layer  to  the  snow  just  below. 

Frad     "     FC     "     ksfc\az/sfc  (9) 

If  it  is  assumed  that  C  and  k  are  constant  to  infinite  depth  and  d  ,  ■  is 
infinitesimal,  the  only  parameters  governing  the  value  of  Tsfc  are  Frad  (Js~  m  ),  C 
(Jm"3K_1),  k  (Jm'V1^1)  and  the  time  since  the  change  in  radiation  conditions,  t  (s). 
From  dimensional  analysis 
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(if 


Tsfc(t)   =  Tsfc(0)  -  nFrad\^/    .  (10) 

The  dimensionless  constant  (n)  was  determined  to  be  equal  to  1.13  based  on 
results  from  a  numerical  model  similar  to  Semtner  (1976),  with  1000  computational 
levels  in  the  snow  layer  and  a  time  step  of  0.01  seconds. 

A  time   scale,   tskin,   can  be   derived  from  Equations   (7)   and   (10)  which 

represents  the  time  required  to  heat  the  skin  layer  enough  so  that  all  the  radiational 

5Tsfc 
heat  imbalance  is  conducted  to  the  lower  layers  and  no  longer  depends  on 

storage  of  heat  in  the  skin  layer: 


(ndskin)2 
Urn     =        „,.       ,  (11) 


K 

where  usnaw  =  "77  is  the  thermal  diffusivity.  At  approximately  t  <  tskin  Equation  (7) 

is  applicable  while  at  t  >  tskin  Equation  (10)  applies.  Note  that  tskin  is  independent 
of  Frad.  Using  a  value  of  1  mm  for  dskin  and  u  =  1.0  x  10  which  are  characteristic 
values  for  feathery  new  snow  (Stull,  1988),  tskjn  has  the  value  of  12.8  seconds.  For  an 
ice  surface,  tskin  equals  1.4  seconds.  Old  snow  has  intermediate  values.  In  nature, 
downward  longwave  radiation  changes  usually  takes  a  few  minutes;  therefore  the 
heat  storage  represented  by  Equation  (7)  is  not  significant  for  downward  radiation. 
Equation  (7)  would  be  more  applicable  to  downward  solar  shortwave  radiation, 
which  is  not  present  in  the  winter.  Shortwave  radiation  can  change  within  a  few 
seconds  as  the  sun  goes  behind  a  cloud  and  penetrates  several  centimeters  into  the 
snow/ice. 


72 


Equation  (10)  is  not  valid  after  the  time,  tturb,  when  Tsfc  has  changed  enough  so 
that  the  surface  turbulent  heat  flux,  Fsens,  is  significant,  assumed  to  be  when  the 
sensible  heat  flux  is  10%  of  the  initial  radiational  imbalance 

F 

■^  rad  .    _x 

F^    >     lo"   •  (12) 

The  turbulent  heat  flux  can  be  approximated  by 

Fsens     =      ATPcpU10CH,  (13) 

where  p  is  the  density  of  air,  cp  is  the  specific  heat  of  air  at  constant  pressure,  AT  is 
the  difference  in  potential  temperature  between  the  surface,  Tsfc,  and  a  reference 
height  near  the  top  of  the  surface  layer,  T10.  U10  is  the  wind  speed  at  the  reference 
height,  and  CH  is  the  heat  transfer  coefficient.  Assuming  that  Equation  (10)  is  valid 
and  T10  remains  equal  to  Tsfc  (t  =  0)  when  t  <  =  tturb,  an  approximate  expression 
for  tturb  can  be  derived  using  Equations  (10),  (12)  and  (13): 

Ck 
lturb      -      (iOPcpU10CHn)2  •  (14) 

This  time  scale  represents  the  time  required  for  the  turbulent  heat  fluxes  to  become 
significant  after  a  change  in  radiation  conditions.  Some  values  for  various  snow/ice 
types  are  listed  in  Table  8.   The  low  values  of  tturb  explain  the  previously  described 

measurements  of  the  close  coupling  between  the  turbulent  heat  fluxes  and  radiation. 
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TABLE  8 
VALUES  OF  Tturb  (seconds)  FOR  U10  =  5  ms"1,  CH  =  1.0  x  10"3 


New  feathery  snow  (Stull,  1988)  0.7 

Old  packed  snow  (Stull,  1988)  44.8 

Ice  (Stull,  1988)  614.0 

Snow  (Semtner,  1976)  33.9 

Ice  (Semtner,  1976)  604.0 

CEAREX  drift  top  3  cm  (estimated)  20.0 
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E.       CONDUCTION  OF  HEAT  THROUGH  SNOW/ICE  (NUMERICAL 
MODEL) 

1.  Model  Description 

The  analytical  models  described  above  are  useful  for  determining  simple 
formulae  for  time  scales  associated  with  radiation  changes  and  surface  temperature 
and  can  be  accurate  for  predicting  Tsfc  for  the  first  few  seconds  or  minutes. 
However,  a  numerical  multi-level  model  is  better  suited  for  predicting  Tsfc  when  the 
radiational  change  is  "felt"  at  lower  levels  in  the  ice/snow,  where  k  and  C  may  not 
be  constant. 

The  numerical  snow/ice  thermodynamic  model  used  was  based  on 
Semtner  (1976).  This  model  numerically  solves  Equation  (8)  for  each  level.  The 
surface  temperature,  Tsfc,  is  determined  by  assuming  a  balance  of  net  radiation, 

turbulent  sensible  and  latent  heat  fluxes,  and  heat  conduction  through  the  ice  at  the 
surface.  The  value  of  Tsfc  for  the  next  time  step  is  based  on  analytical  expressions 
for  the  partial  derivatives  of  each  of  the  fluxes  with  respect  to  a  change  in  Tsfc. 
Semtner  (1976)  prescribed  the  turbulent  fluxes  since  he  was  studying  long  term 
processes  associated  with  ice  growth  and  climate.  For  this  study,  time  dependent 
expressions  for  the  effect  of  turbulent  fluxes  on  Tsfc  have  been  added. 

2.  Model  Time  Scales 

The  purpose  of  this  sub-section  is  to  provide  a  guide  for  choosing  enough 
vertical  layers  so  that  turbulent  fluxes  can  be  accurately  modeled.  Appropriate  time 
scales  associated  with  various  snow/ice  model  parameters  are  described.  These 
time  scales  can  be  used  to  determine  how  many  layers  are  needed  in  the  snow/ice 
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model  to  accurately  predict  turbulent  surface  fluxes.  A  forward  differencing  scheme 
with  time  step  At,  and  model  layer  depth  d,  should  be  numerically  stable  if 

d2 
At   <    ~  .  (15) 

Model  predicted  Tsfc  changes  instantly  in  response  to  changes  in  surface  heat 
balance.  The  analytical  formulae  for  ATsfc  (Equations  (7)  and  (10))  show  that 
actual  changes  in  T\fc  cannot  be  instantaneous.  This  is  because  the  top  of  the  snow 

conducts  and  stores  some  heat.  This  change  is  not  resolved  by  the  numerical  model. 
Therefore,  the  Tsfc  change  after  a  radiation  change  requires  an  amount  of  time, 
tstart,  which  may  be  different  than  the  model  time  step,  At.  For  Semtner's  (1976) 
model  with  prescribed  turbulent  fluxes,  this  time  scale  can  be  approximated  by 


d2 
t  =     -7—  (16) 

siartlO  n~4  V 


"  tstartlO  <   tturtv 

The  time  scale,  tslartll,  would  represent  the  time  required  for  the 
snow/ice  to  absorb  the  amount  of  heat  which  the  model  artificially  created  in  one 
time  step. 


_dl 

tstartll  ~        g  „  (1 ') 


The  version  of  the  snow/ice  model  used  for  this  study  included  the  effects  of 
atmospheric  sensible  and  latent  heat  fluxes  in  the  prediction  of  Tsfc.  In  this  case 
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tstartl2     ~  v  '  (1°) 


dFa    . 

where  Tf-  is  the  predicted  change  in  total  atmospheric  heat  flux  due  to  a  change  in 


dFa 
surface  temperature.    This  equation  reverts  to  Equation  (16)  when  "ZjT"  is  set  to 

Q  1  <.f„ 


zero.  Equation  (18)  is  only  valid  when  tstart  <  tturb. 

A  more  accurate  estimate  of  tstart,  which  can  be  applied  to  cases  when 

Start  Murb>  IS 


i      =  t      MM  .  (i9) 

startl3  start  12    I      p       J       '  v       ' 


where    F      is  a  weighted   average   total   atmospheric  heat  flux   (radiation  and 

turbulence  terms),  which  was  determined  based  on  the  best  fit  to  the  "true"  tstart 
defined  below.  tstartl3  is  always  greater  than  tstartl2  because  turbulent  fluxes  into  the 
atmosphere  buffer  the  change  in  surface  temperature. 

A  summary  of  the  values  of  the  surface  temperature  jump  in  the  initial 
time  step,  ATsfcinit,  tstartl2,  tslartl3  and  a  "true"  tstart  for  various  snow/ice  conditions  is 
presented  in  Table  9.  The  value  of  the  "true"  tstart  is  derived  from  the  amount  of 
time  required  to  have  Tsfc  change  by  ATsfcinjt  based  on  a  1000  snow  layer  model  with 
a  slab  atmosphere.  The  snow,  ice  and  ABL  parameters  used  for  the  "standard  case" 
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TABLE  9 
TEMPERATURE  JUMP,  ATsfcinit,  AND  VARIOUS  PREDICTIONS  OF  THE 

TIME  ELAPSED  IN  THE  FIRST  STEP*  OF  A  NUMERICAL  SNOW/ICE 

MODEL 


Number  of  snow  layers  in  model 
2  10  30  1000 


Standard  Case  (old  snow)  tturb  =  44.8  sees 


ATsfcinjt 

4.80 

4.07 

1.835 

0.773 

0.027 

(Q 

ts,ar.i2  (sees) 

563 

405 

82.4 

14.6 

0.017 

Start  13 

5014 

2133 

148.7 

18.4 

0.017 

Start 

3197 

1838 

148.0 

18.3 

0.017 

New  feathery  snow  tturb  =  0.7  sees 
ATsfcinit  5.76  5.68  5.09  4.05  0.37 


Startl2 

16.7 

16.3 

13.1 

8.2 

0.078 

Startl3 

407 

354 

151 

43.0 

0.078 

Start 

861 

643 

280 

55.0 

0.070 

Snow  replaced  by  ice  tturb  =  614  sees 
ATsfcinit  2.36  1.475  0.370  0.129  0.004 


Start  12 

2510 

985 

61.8 

7.5 

0.007 

Start  13 

5553 

1559 

68.7 

8.1 

0.007 

Start 

4875 

1491 

68.5 

7.6 

0.007 

after  an  instantaneous  change  in  radiation  conditions 
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are  described  in  the  next  subsection.  The  standard  case  represents  typical  central 
Arctic  winter  conditions  with  old  packed  snow.  For  comparison,  a  case  with  Cs  and 
k,.  characteristic  of  new  feathery  snow  (Stull,  1988)  and  a  case  with  Cs  and  k&  equal 
to  ice  values  are  also  shown  in  Table  9.  The  snow  is  assumed  to  be  0.2  m  deep; 
divide  this  into  the  number  of  snow  layers  in  the  model  to  get  the  model  layer 
thickness,  d. 

This  table  shows  that  Equation  (18)  is  accurate  when  tstart  <  tturb  and 

Equation  (19)  is  reasonably  accurate  for  all  cases  except  the  one  and  two  snow-layer 
models.  The  one  snow-layer  model  shown  in  Table  9  has  the  same  parameters  as 
the  OG  case.  Note  that  the  temperature  change  predicted  in  the  first  time  step 
takes  3197  seconds  or  53  minutes,  which  is  much  larger  than  the  numerical  model 
time  step.  In  order  to  fully  resolve  the  period  when  surface  turbulent  fluxes  are 
important,  the  normal  model  time  step,  At,  should  be  smaller  than  tturb.    During 

these  first  50  minutes  the  turbulent  fluxes  have  already  reached  their  maximum 
values. 

The  results  of  this  sub-section  showed  that,  for  snow  conditions  typical  of 
the  CEAREX  drift,  the  Semtner  model  with  one  snow  layer,  as  in  OG,  did  not  have 
enough  vertical  resolution  to  predict  accurately  surface  temperature  and  surface 
fluxes  during  the  first  few  minutes  after  a  change  in  radiation  conditions.  The  first 
time  step  of  the  30  layer  model  required  only  18  seconds;  this  is  small  enough  to 
provide  accurate  representation  of  turbulent  fluxes  and  changes  in  Cd  during  the 

time  period  immediately  after  a  radiation  change.  The  30  layer  version  with 
Equation  (19)  used  to  predict  tstart  was  used  to  derive  the  results  presented  below. 
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3.       Snow/Ice  Model  Predictions  of  Surface  Layer  Stability  Effects  on  Wind 
Stress 

An  analysis  is  performed  in  this  subsection  on  how  surface  layer  stability 

affects  wind  stress.  The  snow/ice  model  simulates  the  time  variation  and  magnitude 

of  wind  stress  following  a  cloud   condition  change.      The   turbulent  heat  and 

momentum  fluxes  into  the  atmosphere  are  based  on  transfer  coefficients  which  have 

been  corrected  for  stability  using  MO  scaling.    The  expressions  for  the  fluxes  are 

implicit  so  that  an  iterative  numerical  procedure  is  required.    The  ABL  above  10 

meters  is  assumed  to  mix  heat  instantly  to  a  prescribed  ABL  depth,  h.    The  wind 

speed  at  10  meters  is  fixed.    These  last  two  assumptions  are  unrealistic  but  the 

purpose  is  to  analyze  only  surface  layer  effects;  later  an  ABL  model  will  be  used  to 

examine  how  the  outer  layer  (the  ABL  above  the  surface  layer)  affects  h,  U10  and 

wind  stress. 

a.       Standard  (Reference)  Case 

A  standard  case  is  defined  based  on  conditions  which  are  typical  of 
the  CEAREX  drift  period.  The  values  of  various  parameters  which  were  used  for 
the  standard  case  are  shown  in  Table  10.  The  initial  snow  and  ice  temperatures 
were  determined  by  running  the  snow/ice  model  for  60  days  from  an  initial  profile 
of  0  C.  This  simulates  the  change  from  summer  to  early  winter  conditions. 

For  the  initialization  run,  the  turbulent  fluxes  were  zero.  The 
downward  longwave  radiation  was  prescribed  to  be  -221.3  Wm  ,  corresponding  to 
the  "previous  sky  temperature"  value  of  -23.19  C  given  in  Table  10.  This  value  was 
chosen  because  the  surface  temperature  becomes  exactly  -20.0  C  after  60  days.  At 
this  surface  temperature,  the  upward  radiation  is  232.9  Wm  .  The  net  radiation 
loss  at  the  surface  is  counteracted  by  a  conductive  heat  flux  through  the  snow  of  11.6 

-2  -2 

Wm   .    This  compares  well  with  the  value  of  12.9  Wm     that  was  estimated  by 
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Parameter 


TABLE  10 
STANDARD  CASE  VALUES 

Symbol  Value  Units 

Constants  (varied  for  non-standard  cases) 


Wind  speed  U 

ABL  depth  h 

Initial  net  radiat.  F 

Snow  depth  ds 

Snow  conductivity  ks 

Snow  heat  capacity  C 

Ice  depth  dj 

Ice  conductivity  k; 

Ice  heat  capacity  C 

Snow  density  ps 


10 


rad 


6 

ms~ 

50 

m 

-80 

Wm'2 

0.2 

m 

0.310 

Jm'V1^1 

6.9  x 

105 

Jm^K"1 

2.0 

m 

2.03 

-1   -1t^-1 

m  s   K 

1.88  > 

.  106 

Jm"3K_1 

330 

Kern-3 

Initial  value  of  time-dependent  parameters 


sfc 


Snow  surface  temp.   T 
ABL  and  10  m  temp.<910  0ABL 

Snow/ice  temps         T(z) 
Previous  sky  temp      T, 


skyl 


-20 

C 

-20 

C 

(see  text) 

-23.19 

c 

Other  constants 

2.3  x  10"3 

dimensionless 

1.5  xlO"3 

dimensionless 

1.4 

kgm"3 

1005 

J  kg^K"1 

-1.7 

C 

Neutral  drag  coeff.  Cdnl0 

Neutral  heat  coeff.  Chnl0 

Air  density  pair 

Air  heat  capacity  C 

Ocean  temp.  Tsea 

Surface  emissivity  t  1.00 

MO  constants  (Guest  and  Davidson,  1987;  Large  and  Pond,  1981,  1982) 
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Maykut  (1982)  for  Arctic  pack  ice  greater  than  0.8  thickness  in  October.  After  60 
days  of  constant  conditions,  the  upper  part  of  the  ice  floe  is  essentially  in 
equilibrium  so  that  the  conductive  fluxes  are  nearly  constant  with  depth,  z,  and  Tsfc 

is  nearly  constant  with  time. 

After  the  initialization  period,  the  downward  radiation  is  instantly 
changed  by  an  amount  represented  by  Frad,  which  is  -80  Wm"2  for  the  standard  case. 

The  standard  case  represents  an  instant  change  from  cloudy  to  clear  conditions. 
The  cloud  layer  is  above  the  ABL  so  that  it  does  not  cool  the  ABL;  the  only  cloud 
effect  is  on  downward  radiation  reaching  the  surface.  The  change  in  radiation  of  -80 
Wm  "  and  the  snow  and  ice  parameters  are  the  same  as  used  by  OG.  The  transfer 
coefficients  represent  values  for  central  Arctic  pack  ice  in  undisturbed  areas  during 
early  winter  based  on  measurements  (OG,  Guest  and  Davidson,  1991a,  and 
unpublished  data). 

The    surface    temperature,    Tsfc,    ABL    temperature,    9ABL    and 
turbulent  sensible  heat  flux,  Fsens,  for  the  24-hour  period  following  a  change  in 

downward  radiation  for  the  standard  case  is  shown  in  Figures  22  and  23. 
Immediately  following  the  radiation  change,  Tsfc  decreases  rapidly  as  the  surface 
experiences  a  net  radiation  loss  (Figure  23).    After  48  seconds  (t  =    tturb)  Tsfc  has 

dropped  enough  so  that  the  turbulent  sensible  heat  flux  is  able  to  start  decreasing 
the  ABL  temperature.  After  about  0.9  hours  01O  decreases  at  the  same  rate  as  Tsfc 
and  the  turbulent  heat  flux,  Fsens,  is  at  a  minimum  (maximum  magnitude).  This  time 
will  be  represented  as  tmax.  After  tmax,  Tsfc  decreases  slower  than  6ABL  so  that  the 
difference  between  them  becomes  less  and  the  magnitude  of  Fsens  decreases,  but 
remains  above  10  Wm    for  10  hours. 
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Figure  22  Slab  Model  Standard  Case  Simulation  Time  Series  of  10  Meter  Air 
Temperature  C  (small  squares),  Surface  Temperature  C  (plusses)  and  Surface 
Sensible  Heat  Flux  Wm'  (asterisks).  The  time  scale  on  this  and  the  following 
figures  in  this  chapter  represents  the  time  after  an  instantaneous  change  in 
downward  radiation. 
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Figure  23  Same  as  Figure  21a  but  Logarithmic  Time  Scale. 
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At  all  times  the  radiation  deficit  at  the  surface,  Rnetsfc,  must  equal 
the  conductive  heat  flux  through  the  snow,  Fsnow,  minus  the  air-snow  turbulent  flux, 
Fsens.  During  the  first  few  seconds  after  the  radiation  change  the  radiation  deficit  is 
counteracted  mostly  by  a  conductive  heat  flux  through  the  snow,  Fsnow  (Figure  24). 
At  time  tmax,  about  55%  of  the  radiation  deficit  is  counteracted  by  Fsnow  while  the 

rest  goes  into  the  atmosphere  (Figures  24  and  25).  The  sensible  turbulent  heat  flux 
remains  below  -10  Wm  for  ten  hours.  After  a  day,  Fsens  has  become  small  and 
Rnetsfc  is  again  almost  totally  counteracted  by  Fsnow.  The  magnitudes  of  both  Fsens 
and  Rnetsfc  will  continue  to  decrease  for  several  days  as  the  surface  cools  off,  but  will 
eventually  approach  an  equilibrium  as  Tsfc  becomes  just  few  degrees  warmer  than 

1  sky 

The  results  from  the  standard  case,  single-layer  ABL,  multi-layer 
snow/ice  model  show  that  surface  layer  fluxes  cause  a  maximum  30%  change  in 
wind  stress  due  to  surface  layer  stability  changes.  Unlike  what  is  implied  by  models 
without  multiple  snow  layers  such  as  used  by  OG,  the  effect  on  wind  stress  is  not 
instantaneous;  some  time  is  required  to  change  the  surface  temperature  before 
stability  is  affected. 

Because  the  atmospheric  model  is  simplistic,  the  actual  magnitude 
of  the  wind  stress  change  in  nature  might  be  different.  However,  the  purpose  of  this 
section  is  to  examine  the  relative  importance  of  various  parameters  in  affecting 
surface  layer  stability.  Inclusion  of  complicated  ABL  physics  would  mask  the 
surface  layer  stability  effects  which  are  the  focus  of  this  chapter.  The  important 
results  here  are  in  the  comparison  of  the  effect  of  various  parameters,  not  in  the 
quantitative  effects  on  wind  stress  during  an  actual  radiation  change  event. 
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Figure  24  Slab  Model  Standard  Case  Simulation  Time  Series  of  Net  Longwave 
Radiation  (plusses),  Conductive  Heat  Flux  at  top  of  Snow  (asterisks)  and  Surface 
Sensible  Heat  Flux  (small  squares).  The  time  scale  is  logarithmic. 
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Figure  25  Slab  Model  Standard  Case  Simulation  Time  Series  of  Conductive  Heat 
Flux  at  Top  of  Snow  Divided  by  Net  Longwave  Radiation  (small  squares)  and 
Negative  Surface  Sensible  Heat  Flux  divided  by  Net  Longwave  Radiation  (large 
squares).  This  figure  shows  how  the  energy  deficit  created  by  the  longwave  cooling 
at  the  surface  is  partitioned  between  turbulent  fluxes  with  the  atmosphere  and 
conductive  fluxes  through  the  snow. 
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In  the  following  subsections,  the  results  from  the  standard  case  will 
be  compared  with  cases  with  different  initial  radiation  imbalances,  wind  speeds, 
ABL  depths,  snow  type  and  snow  depth.  The  following  figures  indicate  the  stability 
wind  stress  effect  with  the  ratio  of  the  10-meter  drag  coefficient  to  the  neutral  10- 
meter  drag  coefficient,  7; —  or  Cdl0  .    Because  the  ten  meter  wind  speed  is  fixed, 

MlnlO 

this  ratio  is  equivalent  to  the  change  in  wind  stress  from  the  initial  neutral  value.  In 
each  of  the  following  cases  one  parameter  is  varied,  while  all  the  other  parameters 
are  set  to  the  standard  case  values. 

b.       Initial  Radiation  Imbalance,  Wind  Speed  and  ABL  Depth  Effects 

The  standard  case  represented  a  change  from  totally  clear  to 
overcast  conditions  above  the  ABL.  Changes  in  the  reverse  direction,  clear  to 
cloudy,  are  just  as  common  as  the  standard  case.  Also,  partial  changes  in  cloud 
cover  are  often  observed  in  the  Arctic.  These  cases  are  modeled  by  setting  the 
value  of  the  initial  radiation  imbalance,  Frad,  to  80  Wm"    (clouding),  40  Wm~ 

(partial  clouding),  -40  Wm"  (partial  clearing)  or  -80  Wm"  (clearing,  standard  case). 
The  surface  layer  effects  of  these  radiation  changes  on  wind  stress  are  shown  in 
Figures  26  and  27.    Note  that  for  the  same  magnitude  Frad,  clearing  produces  a 

greater  effect  than  clouding.  This  is  because  MO  theory  predicts  that  the  heat  flux 
into  the  ABL  will  be  suppressed  by  the  stable  conditions  following  clearing  and 
enhanced  by  the  unstable  conditions  following  clouding.  Therefore,  01O  is  slower  to 
follow  Tsfc  and  the  time  of  maximum  stability  effect,  tmax,  is  greater  for  clearing  cases 
compared  to  clouding  cases.  The  slower  response  of  01O  during  clearing  causes 
larger  surface  layer  temperature  gradients  and  greater  stability  effects. 
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Figure  26      Slab  Model  Simulation  Time  Series  of  Drag  Coefficient  Ratio  after 
Different  Cloud  Change  Events. 
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Figure  27     Same  as  Figure  24  hut  Logarithmic  Time  Scale. 
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Different  wind  speeds  were  tested  in  another  experiment  (Figures 
28  and  29).  The  largest  relative  (Figures  28  and  29)  and  absolute  (Figure  30) 
changes  in  wind  stress  after  a  radiation  change  event  occur  during  the  lightest  winds. 
The  maximum  wind  stress  reduction  was  only  9%  for  the  9  ms"1  wind  case, 
compared  to  a  30%  reduction  for  the  standard  6  ms"1  case  and  a  58%  reduction  for 
the  5  ms  case.  At  wind  speeds  below  4.5  ms"1,  the  wind  stress  becomes  zero  after 
about  an  hour  (not  shown)  because  MO  theory  predicts  complete  suppression  of 
turbulence  and  heat  and  momentum  fluxes.  The  CEAREX  drift  period  from  2 
November  1988  to  5  November  1988  (Figure  5)  is  an  example  of  a  situation  when 
fluxes  have  been  completely  suppressed  during  light  winds.  The  results  show  that 
surface  layer  stability  effects  on  wind  stress  are  important  in  the  central  Arctic 
during  average  wind  conditions,  i.e.    U10   =   5  ms"1.  During  high  wind  periods,  U10 

>  10  ms"  ,  the  surface  layer  stability  effect  is  almost  negligible. 

The  ABL  depth,  h,  is  important  to  surface  layer  stability  changes 
after  a  radiation  change  (Figures  31  and  32).  The  temperature  of  the  ABL,  01O, 
responds  slower  to  changes  in  the  surface  temperature,  Tsfc,  when  h  is  large.    In 

nature,  wind  stress  effects  for  an  initial  large  h  are  less  than  what  is  shown  in  Figures 
31  and  32  because  during  these  stable  periods  a  new  low-level  inversion  would  form 
and  h  would  become  much  less  after  a  few  minutes. 

c.       Snow  Characteristics  Effects 

Snow  characteristics  effects  must  also  be  considered.  For  such 
examinations,  different  standard  case  initializations  were  necessary  because 
different  snow/ice  depths  or  conductivities  result  in  different  equilibrium 
temperature  profiles.  New  60  day  initializations  were  based  on  the  same  downward 
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Figure  28     Slab  Model  Simulation  Time  Series  of  Drag  Coefficient  Ratio  For 
Different  Wind  Speeds. 
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Figure  29     Same  as  Figure  28  but  Logarithmic  Time  Scale. 
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Figure  30  Slab  Model  Simulation  Time  Series  of  u,  for  Different  Wind  Speeds. 
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Figure  31     Slab  Model   Simulation  Time  Series  of  Drag  Coefficient  Ratio  for 
Different  ABL  Depths. 
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Figure  30  Slab  Model  Simulation  Time  Series  of  u»  for  Different  Wind  Speeds. 
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Figure  31     Slab   Model  Simulation  Time  Series  of  Drag  Coefficient  Ratio  for 
Different  ABL  Depths. 
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Figure  33  Temperature  (X-Axis)  as  a  Function  of  Snow  and  Ice  Depths  (Y-Axis) 
after  60-day  Slab  Model  Simulation.  Different  simulations  for  various  snow  types 
and  thicknesses  are  shown. 
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Figure  34  Errors  Resulting  from  Using  Too  Few  Ice  Levels  in  the  Numerical 
Snow/ice  Model.  Slab  model  drag  coefficient  ratio  time  series  results  using  snow 
depth  of  5  cm  or  1/4  or  the  standard  case  values.  Two  different  ABL  depths,  h,  and 
two  different  numbers  of  ice  layers  are  shown. 
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layers.  Using  more  than  30  ice  layers  results  in  virtually  identical  results  as  the  30 
ice  layer  version  shown  in  Figure  34;  30  ice  layers  is  enough  resolution  for  accurate 
modeling  of  all  cases.  The  error  was  much  less  for  snow  depths  of  0.2  m  or  more. 
When  this  error  was  significant,  the  number  of  ice  layers  was  increased  from  2  to  30. 
The  surface  of  the  pack  ice  can  consist  of  ice  or  various  types  of 
snow.  Different  surfaces  are  modeled  by  changing  the  volumetric  heat  capacity,  Cs, 
and  thermal  conductivity,  ks,  of  the  upper  0.2  m  of  the  floe  (Figures  35  and  36).  The 
values  Cs  and  ks  used  for  the  different  surfaces  are  the  same  as  used  for  Table  10 
with  an  additional  category  represented  by  C&  and  ks  values  between  the  old  packed 
and  new  feathery  values,  representing  an  intermediate  snow  type.  The  results  in 
Figures  35  and  36  show  that  the  wind  stress  after  a  radiation  change  is  strongly 
dependent  on  snow  type.  Light  new  snow  is  a  better  insulator  and  has  less 
volumetric  heat  capacity  than  denser  old  snow;  both  characteristics  cause  Tsfc  to 

change  more  rapidly  and  generate  larger  surface  layer  stability  effects  for  light  snow 
vs.  dense  snow.  Volumetric  heat  capacity,  Cs,  is  proportional  to  density,  ps.  Thermal 
conductivity,  ks,  is  proportional  to  psnow  (Stull,  1988).  Initially,  the  change  in  Tsfc  is 
therefore  proportional  to  psnow ~  ,  using  Equation  (10).  After  time  tturb,  the 
difference  in  Tsfc  changes  between  ice  types  is  less  than  implied  by  Equation  (10), 
but  still  significant  enough  to  cause  large  differences  in  the  maximum  magnitudes  of 
heat  flux  and  wind  stress  reduction  stability  effect.  For  example,  after  one  hour, 
which  is  approximately  tmax,  wind  stress  for  ice,  intermediate  snow  and  old  packed 

snow  has  been  reduced  by  12%,  44%  and  63%,  respectively,  compared  to  30%  for 
the  standard  case. 

The  insulating  characteristics  of  snow  allow  significant  stability 
effects  on  wind  stress  to  occur  over  the  central  Arctic.    Liquid  water  has  a  greater 
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volumetric  heat  capacity  than  snow  and  a  much  greater  thermal  conductivity  (due  to 
turbulence);  therefore  the  stability  effects  that  have  been  discussed  could  not  occur 
over  open  ocean  regions. 

The  snow  cover  depth,  ds,  of  pack  ice  can  vary  from  no  snow  to  tens 

of  centimeters.  The  surface  layer  stability  effect  on  wind  stress  for  three  different 
snow  depths  and  no  snow  is  shown  in  Figures  37  and  38.  For  the  standard  case,  with 
ABL  depth  equal  to  50  meters,  there  is  virtually  no  difference  between  the  0.2  m 
deep  snow  and  the  0.8  m  deep  snow.  The  0.05  m  thick  snow  case  shows  small 
differences  from  the  thicker  cases  after  0.5  hours.  The  tmax  occurs  25%  sooner  and 

there  is  about  a  20%  less  stability  effect  after  one  hour  for  the  0.05  m  snow  case 
compared  to  the  standard  0.2  m  case.  These  are  minor  relative  differences 
compared  to  the  solid  ice  case  which  had  much  less  stability  effect  during  the  first  4 
hours  than  all  the  snow  cases. 

<L       Physical  Time  Scales 

When  the  height  of  the  ABL  is  increased  to  1000  m,  there  is  more 
difference  in  the  results  for  the  different  snow  depths  (Figure  39)  than  the  h  =  50 
case  (Figure  36).  Initially,  the  stability  effects  are  identical  as  long  as  some  snow  is 
present.  After  1.5  hours,  the  0.05  snow  thickness  case  has  10%  less  stability  effect 
than  the  0.2  m  or  0.8  m  snow  thickness  cases.  Note  that  this  time  scale  (1.5  hours) 
was  identical  for  both  h  =  50  m  and  h  =  1000  m  cases.  This  is  also  the  time  after 
which  the  2  ice  layer  model  begins  to  show  some  error  (Figure  34).  This  time 
represents  the  time  elapsed,  tdeep,  before  the  ice  layer  below  the  snow  can  "feel"  the 

different  radiation  conditions  and  begin  to  affect  the  temperature  within  and  on  top 
of  the  snow  layer.  The  value  of  tdeep  depends  on  the  thermal  diffusivity,  ^snow,  and 
depth,  dsnow,  of  the  snow  layer.  From  dimensional  analysis 
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Figure  36    Same  as  Figure  35  but  Logarithmic  Time  Scale. 
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volumetric  heat  capacity  than  snow  and  a  much  greater  thermal  conductivity  (due  to 
turbulence);  therefore  the  stability  effects  that  have  been  discussed  could  not  occur 
over  open  ocean  regions. 

The  snow  cover  depth,  ds,  of  pack  ice  can  vary  from  no  snow  to  tens 

of  centimeters.  The  surface  layer  stability  effect  on  wind  stress  for  three  different 
snow  depths  and  no  snow  is  shown  in  Figures  37  and  38.  For  the  standard  case,  with 
ABL  depth  equal  to  50  meters,  there  is  virtually  no  difference  between  the  0.2  m 
deep  snow  and  the  0.8  m  deep  snow.  The  0.05  m  thick  snow  case  shows  small 
differences  from  the  thicker  cases  after  0.5  hours.  The  tmav  occurs  25%  sooner  and 

UlcU 

there  is  about  a  20%  less  stability  effect  after  one  hour  for  the  0.05  m  snow  case 
compared  to  the  standard  0.2  m  case.  These  are  minor  relative  differences 
compared  to  the  solid  ice  case  which  had  much  less  stability  effect  during  the  first  4 
hours  than  all  the  snow  cases. 

d.       Physical  Time  Scales 

When  the  height  of  the  ABL  is  increased  to  1000  m,  there  is  more 
difference  in  the  results  for  the  different  snow  depths  (Figure  39)  than  the  h  =  50 
case  (Figure  36).  Initially,  the  stability  effects  are  identical  as  long  as  some  snow  is 
present.  After  1.5  hours,  the  0.05  snow  thickness  case  has  10%  less  stability  effect 
than  the  0.2  m  or  0.8  m  snow  thickness  cases.  Note  that  this  time  scale  (1.5  hours) 
was  identical  for  both  h  =  50  m  and  h  =  1000  m  cases.  This  is  also  the  time  after 
which  the  2  ice  layer  model  begins  to  show  some  error  (Figure  34).  This  time 
represents  the  time  elapsed,  tdeep,  before  the  ice  layer  below  the  snow  can  "feel"  the 

different  radiation  conditions  and  begin  to  affect  the  temperature  within  and  on  top 
of  the  snow  layer.  The  value  of  tdeep  depends  on  the  thermal  diffusivity,  ^snow,  and 
depth,  dsnow,  of  the  snow  layer.  From  dimensional  analysis 
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Figure  37       Slab  Model  Simulation  Time  Series  of  Drag  Coefficient  Ratio  for 
Different  Snow  Depths. 
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where  Q  is  a  dimensionless  constant.  The  value  of  Q  depends  on  what  magnitude 

of  ice  layer  effect  is  deemed  significant.  If  significant  is  defined  as  a  1%  effect  on 
the  wind  stress,  then  the  constant  C2  has  a  value  of  about  0.3.  If  significant  is 
defined  as  a  10%  effect  on  the  wind  stress,  then  the  constant  C2  has  a  value  of  about 
one  (see  Figure  39).  Analytical  (Equation  (20)  with  C^  =  1.0)  and  numerical 
(Figure  39)  results  using  both  imply  that  tdeep  is  approximately  equal  to  1.5  and  24 
hours  for  dsnow  =  0.05  m  and  0.2  m  respectively. 

Because  tdeep  equals  24  hours  for  0.2  m  thick  snow  and  384  hours 
for  0.8  m  thick  snow,  a  significant  difference  between  the  two  cases  should  be 
noticeable  after  24  hours.  There  is  an  obvious  difference  at  24  hours  for  the  h  = 
1000  m  case  (Figure  39)  but  it  is  less  obvious  for  the  h  =  50  m  case  (Figure  36).  The 
reason  that  the  effect  of  depth  is  not  very  noticeable  for  the  h  =  50  m  case  is 
because  after  24  hours  the  ABL  has  already  adjusted  to  the  changes  in  surface 
temperature  and  heat  fluxes  are  small,  i.e.  t  >  >  tmax.    The  stability  effect  has 

become  small,  so  relative  differences  of  10%  correspond  to  a  small  absolute 
difference  in  stability  effect.  In  order  for  the  snow  depth  to  have  an  effect  on  wind 
stress,  tdeep  must  be  less  than  an  atmospheric  time  scale,  tABL,  which  represents  the 

amount  of  time  required  before  the  atmosphere  has  adjusted  to  the  change  in 
surface  temperature.  According  to  the  simple  ABL  slab  model,  at  any  time,  t,  after 
a  surface  radiation  change,  the  following  equality  will  exist 
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Figure  39  Slab  Model  Simulation  Time  Series  of  Drag  Coefficient  Ratio  for 
Different  Snow  Depths  Using  an  ABL  Depth  of  1000  m  or  20  Times  The  Standard 
Case  Value. 
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where  C2  is  a  dimensionless  constant.  The  value  of  C2  depends  on  what  magnitude 

of  ice  layer  effect  is  deemed  significant.  If  significant  is  defined  as  a  1%  effect  on 
the  wind  stress,  then  the  constant  C2  has  a  value  of  about  0.3.  If  significant  is 
defined  as  a  10%  effect  on  the  wind  stress,  then  the  constant  C2  has  a  value  of  about 
one  (see  Figure  39).  Analytical  (Equation  (20)  with  C2  =  1.0)  and  numerical 
(Figure  39)  results  using  both  imply  that  tdeep  is  approximately  equal  to  1.5  and  24 
hours  for  dsnow  =  0.05  m  and  0.2  m  respectively. 

Because  tdeep  equals  24  hours  for  0.2  m  thick  snow  and  384  hours 
for  0.8  m  thick  snow,  a  significant  difference  between  the  two  cases  should  be 
noticeable  after  24  hours.  There  is  an  obvious  difference  at  24  hours  for  the  h  = 
1000  m  case  (Figure  39)  but  it  is  less  obvious  for  the  h  =  50  m  case  (Figure  36).  The 
reason  that  the  effect  of  depth  is  not  very  noticeable  for  the  h  =  50  m  case  is 
because  after  24  hours  the  ABL  has  already  adjusted  to  the  changes  in  surface 
temperature  and  heat  fluxes  are  small,  i.e.  t  >  >  tmax.    The  stability  effect  has 

become  small,  so  relative  differences  of  10%  correspond  to  a  small  absolute 
difference  in  stability  effect.  In  order  for  the  snow  depth  to  have  an  effect  on  wind 
stress,  td      must  be  less  than  an  atmospheric  time  scale,  tABL,  which  represents  the 

amount  of  time  required  before  the  atmosphere  has  adjusted  to  the  change  in 
surface  temperature.  According  to  the  simple  ABL  slab  model,  at  any  time,  t,  after 
a  surface  radiation  change,  the  following  equality  will  exist 
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where  aTabl  is  the  change  in  ABL  potential  temperature,  and  Fsens  is  the  average 
sensible  heat  flux.     By  picking  appropriate  scales  for  ATABL  and  Fsens  a  non- 

dimensional  time,  t  ,  can  be  derived.  Since  the  ABL  temperature  change  is  dictated 
by  the  change  in  downward  radiation,  a  reasonable  scale  for  ATABL  is  ATsky  which 

expresses  the  radiation  change  in  terms  of  change  in  sky  temperature.    A  scale  to 
use  for  Fsens  is  the  maximum  sensible  heat  flux,  Fsensmax.    Using  these  scales  the 

nondimensional  time  can  be  defined  as 


(22) 


(23) 

At     t    =  1  the  ABL  has  adjusted  to  the  change  in  surface  conditions.  This  occurs 
at    *  =  Ubl-  For  the  standard  case,  tABL  equals  15  hours,  which  corresponds  well 
with  the  time  period  over  which  wind  stress  is  affected  by  surface  stability  after  a 
change  in  radiation  conditions. 

The  model  results  show  that  the  maximum  sensible  heat  fluxes 
occur  at  approximately 


1 
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The  dimensional  time  of  the  maximum  sensible  fluxes,  tmax,  therefore  occurs  at 

t 
I 

15 


ABL 
tmax     ~       ^S.      ■  (25) 


108 


change  represents  a  time  change  different  from  the  model  time  step.  An  expression 
for  this  time  period  was  derived. 

Model  simulations  show  that,  unlike  results  from  single  snow  layer 
surface  models  such  as  that  used  by  OG,  the  surface  fluxes  require  some  time  to 
"kick  in"  after  a  cloud  cover  change.  Analytical  expressions  for  this  time  period,  as 
well  as  expressions  for  the  time  of  maximum  surface  layer  stability  effect  and  the 
entire  time  period  over  which  surface  layer  stability  effects  are  important  were 
derived.  Assuming  constant  ABL  depth,  the  surface  layer  stability  change  is 
inversely  related  to  wind  speed,  as  expected.  At  wind  speeds  greater  that  9  ms"  the 
surface  layer  stability  effect  is  negligible.  The  specified  value  of  ABL  depth  can 
strongly  affect  the  surface  flux  changes.  The  temperature  of  a  deep  ABL  cannot 
respond  quickly  to  changes  in  surface  temperature  and  therefore  flux  magnitudes 
and  stability  effects  are  larger  and  stay  large  for  longer  periods  than  for  shallow 
ABLs. 

The  thermodynamic  characteristics  of  an  ice  floe  surface  can  vary  from 
mostly  solid  ice  to  light  feathery  snow.  The  type  of  surface  strongly  affects  surface 
fluxes  following  a  cloud  cover  change.  However,  the  depth  of  the  snow  layer  on  top 
of  the  pack  ice  may  or  may  not  affect  surface  fluxes.  Analytical  expressions  based 
on  snow  and  ABL  characteristics  can  be  used  to  determine  when  snow  depth  is  an 
important  parameter. 
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TABLE  11. 
SUMMARY  OF  TIME  SCALES  RELATED  TO  EFFECTS  OF  h  AND  ds 

SURFACE  LAYER  STABILITY 


ON 


Condition 


Influence  on  wind  stress 


1 


t   <  tturb 

*     <  *deep 

t     <  th 

1  '■max 

t     >  tABL 

*deep  >     *ABL 


*ABL     >     ^deep 


>     t. 


ldeep 


<     t. 


small  stability  effects 

dsnowhasno  effect 

h  has  no  effect 

time  of  maximum  effect 

again  small  stability  effects 

dsnow  never  has  an  effect 

dsnow  has  small  effect  on  stability 


dsnow  has  large  effect  on  stability 


i 
"no"  means   <  1%  effect 

"small"  means  >  1%  and  <  10%  effect 

"large"  means  >  10%  effect 

Assuming      tdeep  <   t  <   tABL 
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change  represents  a  time  change  different  from  the  model  time  step.  An  expression 
for  this  time  period  was  derived. 

Model  simulations  show  that,  unlike  results  from  single  snow  layer 
surface  models  such  as  that  used  by  OG,  the  surface  fluxes  require  some  time  to 
"kick  in"  after  a  cloud  cover  change.  Analytical  expressions  for  this  time  period,  as 
well  as  expressions  for  the  time  of  maximum  surface  layer  stability  effect  and  the 
entire  time  period  over  which  surface  layer  stability  effects  are  important  were 
derived.  Assuming  constant  ABL  depth,  the  surface  layer  stability  change  is 
inversely  related  to  wind  speed,  as  expected.  At  wind  speeds  greater  that  9  ms"  the 
surface  layer  stability  effect  is  negligible.  The  specified  value  of  ABL  depth  can 
strongly  affect  the  surface  flux  changes.  The  temperature  of  a  deep  ABL  cannot 
respond  quickly  to  changes  in  surface  temperature  and  therefore  flux  magnitudes 
and  stability  effects  are  larger  and  stay  large  for  longer  periods  than  for  shallow 
ABLs. 

The  thermodynamic  characteristics  of  an  ice  floe  surface  can  vary  from 
mostly  solid  ice  to  light  feathery  snow.  The  type  of  surface  strongly  affects  surface 
fluxes  following  a  cloud  cover  change.  However,  the  depth  of  the  snow  layer  on  top 
of  the  pack  ice  may  or  may  not  affect  surface  fluxes.  Analytical  expressions  based 
on  snow  and  ABL  characteristics  can  be  used  to  determine  when  snow  depth  is  an 
important  parameter. 
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V.  CASE  STUDIES  OF  THE  EFFECT  OF  CLOUDS  ON  WIND  STRESS 

The  influence  of  clouds  on  the  wind  stress  resulting  from  their  effect  on 
stability  dependent  surface  layer  drag  coefficient  was  examined  in  the  previous 
chapter.  In  that  examination,  the  relationship  between  downward  radiation  changes 
due  to  cloud  cover  changes  and  the  surface  stability  changes  was  considered  for 
Type  R  situations.  Because  the  wind  speed  is  the  primary  factor  in  the  wind  stress, 
the  effect  of  clouds  on  surface  wind  speeds  as  well  as  stability  requires  examination, 
based  on  Arctic  boundary  layer  characteristics.  Clouds  can  affect  the  surface  wind 
vector  by  being  radiatively  cooled  at  the  top  and  by  influencing  the  entrainment  of 
air  into  the  ABL.  Entrainment  and  cloud  top  cooling  within  the  ABL  are  non-Type 
R  phenomena. 

Unlike  the  simple  Monin-Obukhov  surface  layer  theory  used  to  determine 
radiation  effects  on  surface  layer  stability  in  Chapter  IV,  the  combined  effect  of 
these  processes  is  too  complicated  to  describe  with  analytical  equations  or  simple 
numerical  models.  Therefore,  this  chapter  will  examine  the  combined  effect  of  all 
physical  processes  which  may  link  wind  stress  to  cloud  conditions  by  examining  the 
results  of  numerical  ABL  simulations,  with  cloud  and  radiation  physics  included, 
and  in  situ  observations  from  the  CEAREX  drift. 

Three  case  studies  are  presented  which  use  numerical  model  results, 
initialized  and  verified  using  CEAREX  observations,  to  quantify  the  various  cloud 
effects  on  wind  stress.  The  case  studies  will  be  used  to  show  realistic  examples  of 
how  various  cloud  conditions  and  changes  in  cloud  conditions  can  influence  wind 
stress  in  the  Arctic.    The  purpose  of  the  case  studies  presented  here  is  to  use  the 
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longwave  radiation  with  model  predictions  based  on  a  concurrent  rawinsonde 
sounding. 

In  this  study,  the  original  model  was  found  to  significantly  underpredict 
the  downward  longwave  radiation  at  the  surface  during  clear  winter  Arctic 
conditions.  When  Thompson  and  Burk  (1990)  used  this  radiation  code  to  model  an 
Arctic  front,  they  were  forced  to  set  the  surface  temperature  to  a  constant  value 
since  the  incorrect  radiation  was  causing  a  snow/ice  relaxation  scheme  to  predict 
surface  temperature  values  that  were  too  low.  The  problem  only  occurred  at 
temperatures  below  0  C  and  only  became  obvious  when  the  model  was  compared  to 
the  CEAREX  drift  measurements.  For  this  study,  the  radiation  code  was  modified 
(Burk,  personal  communication)  and  the  model  radiation  values  compared  well  with 
the  clear  sky  measurements  during  cold  clear  weather  on  the  average.  There  was 
considerable  scatter  remaining. 

The  radiation  model  also  underpredicted  the  absorption/emission  of 
longwave  radiation  in  Arctic  winter  clouds.  Arctic  clouds  contain  a  relatively  small 
amount  of  particulate  water  but  are  observed  to  be  quite  black.  The  model  is  based 
on  lower  latitude  measurements,  as  are  all  radiation  schemes,  and  does  not  account 
for  ice  particles  and  variations  in  particle  size  spectra  which  may  occur  in  the  Arctic. 
Information  on  the  characteristics  of  cloud  particulates  was  not  available 
concurrently  with  the  data  used  for  this  dissertation;  therefore  a  detailed  study  of 
the  radiation  characteristics  of  Arctic  clouds  was  not  attempted.  By  increasing  the 
absorptivity  of  particulate  water  in  the  clouds  by  a  factor  of  30,  a  good  match 
(average  difference  3.5  Wm  )  between  measurements  of  downward  longwave 
radiation  and  model  predictions  during  cloudy  conditions  was  achieved.  The  match 
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potential  temperature,  0,  (2)  total  water  content,  q,  (3)  turbulent  kinetic  energy, 
TKE  (4)  the  variance  of  0,  (5)  the  variance  of  qt  (6)  the  cross  correlation  of  0,  and  q, 
(7)  the  east  component  of  mean  momentum,  u,  and  (8)  the  north  component  of 
mean  momentum,  v.  Other  second-moment  equations  are  solved  with  diagnostic 
equations.  The  turbulence  closure  occurs  at  level  2.5  and  level  3  using  Mellor  and 
Yamada  (1974)  terminology.  The  second-order  closure  physics  allows  vertical 
virtual  temperature  fluxes  to  occur  even  when  the  gradient  is  zero  (i.e.  a  mixed 
layer)  as  long  as  turbulence  is  present. 

2.  Cloud  Model 

By  using  thermodynamic  turbulence  statistics  generated  by  the  model, 
the  method  of  Sommeria  and  Deardorff  (1977)  is  used  to  predict  the  cloud  cover  at 
each  level.  With  no  turbulence,  the  cloud  cover  instantly  changes  from  zero  to  total 
when  the  mean  humidity  reaches  saturation.  But  with  increasing  turbulence,  there 
is  an  expanding  mean  humidity  range  around  100%  relative  humidity  where  partial 
cloudiness  is  predicted. 

When  partial  cloudiness  exists  at  several  levels,  random  overlapping  is 
assumed  (Harshvardhan  et  al,  1987;  Burk  and  Thompson,  1989).  The  radiation  at 
each  level  is  partitioned  into  the  clear  and  cloudy  air  components. 

3.  Radiation  Model 

Radiative  cooling/heating  is  based  on  a  broadband  radiative  scheme 
with  transmission  functions  at  each  level  determined  from  the  amount  of  water 
vapor,  liquid  water,  carbon  dioxide  and  aerosol  (Oliver  et  al.,  1978;  Lewellen  et  al., 
1976).  The  same  radiation  scheme  has  been  tested  for  a  variety  of  conditions, 
although  not  under  Arctic  conditions  to  my  knowledge.  The  CEAREX 
measurements  provided  a  good  opportunity  to  compare  measurements  of  downward 
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longwave  radiation  with  model  predictions  based  on  a  concurrent  rawinsonde 
sounding. 

In  this  study,  the  original  model  was  found  to  significantly  underpredict 
the  downward  longwave  radiation  at  the  surface  during  clear  winter  Arctic 
conditions.  When  Thompson  and  Burk  (1990)  used  this  radiation  code  to  model  an 
Arctic  front,  they  were  forced  to  set  the  surface  temperature  to  a  constant  value 
since  the  incorrect  radiation  was  causing  a  snow/ice  relaxation  scheme  to  predict 
surface  temperature  values  that  were  too  low.  The  problem  only  occurred  at 
temperatures  below  0  C  and  only  became  obvious  when  the  model  was  compared  to 
the  CEAREX  drift  measurements.  For  this  study,  the  radiation  code  was  modified 
(Burk,  personal  communication)  and  the  model  radiation  values  compared  well  with 
the  clear  sky  measurements  during  cold  clear  weather  on  the  average.  There  was 
considerable  scatter  remaining. 

The  radiation  model  also  underpredicted  the  absorption/emission  of 
longwave  radiation  in  Arctic  winter  clouds.  Arctic  clouds  contain  a  relatively  small 
amount  of  particulate  water  but  are  observed  to  be  quite  black.  The  model  is  based 
on  lower  latitude  measurements,  as  are  all  radiation  schemes,  and  does  not  account 
for  ice  particles  and  variations  in  particle  size  spectra  which  may  occur  in  the  Arctic. 
Information  on  the  characteristics  of  cloud  particulates  was  not  available 
concurrently  with  the  data  used  for  this  dissertation;  therefore  a  detailed  study  of 
the  radiation  characteristics  of  Arctic  clouds  was  not  attempted.  By  increasing  the 
absorptivity  of  particulate  water  in  the  clouds  by  a  factor  of  30,  a  good  match 
(average  difference  3.5  Wm  )  between  measurements  of  downward  longwave 
radiation  and  model  predictions  during  cloudy  conditions  was  achieved.  The  match 
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is  good  because  both  adjusted  model  and  observations  have  virtually  black  cloud 
bases  if  the  cloud  is  thicker  than  100  m. 

During  clear  sky  conditions,  the  model  predictions  of  downward 
longwave  radiation  were  less  accurate,  with  an  average  difference  of  22  Wm  ,  even 
after  the  modifications.  It  is  doubtful  that  other  radiation  schemes  could  be 
significantly  more  accurate  during  clear  skies.  This  is  because  large  (>  30  Wm  ) 
differences  in  measured  downward  radiation  occurred  between  soundings  which  had 
very  similar  humidity  and  temperature  profiles.  No  radiation  model  could  predict 
these  differences  based  on  the  rawinsonde  profile  alone. 

Two  explanations  for  the  large  variations  in  measured  downward 
longwave  radiation  for  similar  rawinsonde  profiles  are  (1)  Varying  amounts  of  small 
ice  crystal  concentrations  and  (2)  Different  temperature/humidity  characteristics 
above  the  highest  rawinsonde  point.  OG  and  Curry  (1983)  show  how  ice  crystals  or 
"diamond  dust"  are  important  to  radiative  transfer  in  the  Arctic.  Variable 
concentrations  of  diamond  dust  were  personally  observed  during  the  CEAREX  drift 
period.  Assuming  ice  crystal  concentrations  are  also  variable  throughout  the  entire 
atmosphere,  it  is  reasonable  to  expect  that  the  longwave  radiation  fluctuates  due  to 
the  ice  crystals.  The  additional  downward  longwave  radiation  at  the  surface  due  to 
ice  crystals  throughout  the  troposphere  can  vary  from  15  to  80  Wm  (Curry  et  al, 
1989). 

The  model  radiation  predictions  at  the  surface  during  clear  weather  were 
quite  sensitive  to  the  specification  of  temperature  and  humidity  in  the  stratosphere, 
above  the  top  level  of  most  of  the  rawinsonde  profiles.  The  temperature  at  the  100 
mb  level  was  adjusted  so  that  the  snow/ice  model  produced  the  measured  clear  sky 
surface  temperature.    The  adjustment  accounts  for  the  unknown  amounts  of  ice 
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longwave  radiation  over  pack  ice  in  winter  (Chapter  IV).   Three  choices  for  initial 
surface  temperature  are 

(Tl)  the  value  produced  by  the  snow/ice  model  using  the 

measured  radiation, 

(T2)  the  value  produced  by  the  snow/ice  model  using  the 

radiation  value  produced  by  the  model  radiation  code 
and  the  measured  temperature  and  humidity  profiles  or 

(T3)  the  measured  value. 

These  correspond  to  the  following  values  of  downward  radiation: 

(Rl)  the  measured  value, 

(R2)  the  value  given  by  the  radiation  code 

(R3)  the  value  required  by  the  snow/ice  model  to  produce 

the  measured  surface  temperature. 

Because  the  measured,  radiation  code  and  ice/snow  model  values  differ, 
the  longwave  radiation  and  surface  temperature  values  must  be  paired,  i.e.  Rl  must 
be  paired  with  Tl,  R2  with  T2  etc.,  for  the  initialization  to  have  an  energy  balance  at 
the  surface.  Once  the  model  run  begins,  the  surface  temperature  is  determined  by 
the  same  method  that  was  used  to  determine  T2  during  the  initialization  procedure. 
If  either  Tl  or  T3  is  used  at  the  start  of  the  model,  the  model  surface  temperature 
will  quickly  jump  toward  T2,  again  because  of  an  energy  imbalance  at  the  surface. 
Therefore  the  only  pair  which  will  not  create  a  rapid  change  in  surface  temperature 
due  to  a  surface  energy  imbalance  at  the  start  of  the  model  run  is  R2  and  T2.  This 
is  unsatisfactory  because  the  (model-generated)  surface  temperature  is  determined 
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half  of  the  surface  geostrophic  value.  For  this  case,  a  low-level  thermal  wind  was 
assumed,  and  the  geostrophic  winds  above  h  were  set  to  the  true  winds  measured 
with  a  rawinsonde. 

Advection  could  be  simulated  by  having  tendency  terms  for  the 
prognostic  variables.  All  tendencies  were  set  to  zero  in  the  case  studies  shown  here. 
An  exception  was  that  sometimes  a  sudden  moisture  tendency  was  imposed  to  force 
clouds  to  form  or  dissipate.  Any  moisture  tendencies  used  will  be  discussed  in  the 
individual  case  study  discussions. 

5.       Model  Initialization 

The  1-D  ABL  model  is  initialized  from  a  rawinsonde  or  model  profile  of 
temperature,  humidity  and  geostrophic  wind  vector  and  a  modeled  or  specified 
surface  temperature.  But  a  classic  problem  arises:  the  values  of  different  measured 
parameters  are  not  in  equilibrium  with  the  model  physics.  Usually  this  problem  is 
approached  by  having  an  adjustment  or  spin-up  period  before  the  actual  model  run; 
this  method  was  used  for  several  parameters.  But  for  comparison  between  the 
modeled  effect  of  clouds  on  wind  stress  and  measurements,  the  spin-up  should  not 
substantially  change  the  specification  of  clouds,  the  initial  surface  temperature  and 
the  low  level  temperature  profile.  The  CEAREX  drift  measurements  provided 
further  information,  such  as  measurements  of  longwave  radiation,  wind  stress  and 
heat  flux,  which  should  allow  the  initializations  to  match  reality  to  a  greater  extent. 
However,  the  detailed  CEAREX  data  provides  the  dilemma  of  having  too  much 
initialization  information. 

To  illustrate  this  problem,  consider  the  initialization  of  surface 
temperature.      The   surface   temperature   is   strongly   related   to   the   downward 
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longwave  radiation  over  pack  ice  in  winter  (Chapter  IV).   Three  choices  for  initial 
surface  temperature  are 

(Tl)  the  value  produced  by  the  snow/ice  model  using  the 

measured  radiation, 

(T2)  the  value  produced  by  the  snow/ice  model  using  the 

radiation  value  produced  by  the  model  radiation  code 
and  the  measured  temperature  and  humidity  profiles  or 

(T3)  the  measured  value. 

These  correspond  to  the  following  values  of  downward  radiation: 

(Rl)  the  measured  value, 

(R2)  the  value  given  by  the  radiation  code 

(R3)  the  value  required  by  the  snow/ice  model  to  produce 

the  measured  surface  temperature. 

Because  the  measured,  radiation  code  and  ice/snow  model  values  differ, 
the  longwave  radiation  and  surface  temperature  values  must  be  paired,  i.e.  Rl  must 
be  paired  with  Tl,  R2  with  T2  etc.,  for  the  initialization  to  have  an  energy  balance  at 
the  surface.  Once  the  model  run  begins,  the  surface  temperature  is  determined  by 
the  same  method  that  was  used  to  determine  T2  during  the  initialization  procedure. 
If  either  Tl  or  T3  is  used  at  the  start  of  the  model,  the  model  surface  temperature 
will  quickly  jump  toward  T2,  again  because  of  an  energy  imbalance  at  the  surface. 
Therefore  the  only  pair  which  will  not  create  a  rapid  change  in  surface  temperature 
due  to  a  surface  energy  imbalance  at  the  start  of  the  model  run  is  R2  and  T2.  This 
is  unsatisfactory  because  the  (model-generated)  surface  temperature  is  determined 
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with  a  different  method  from  the  rawinsonde  air  temperature.  This  will  usually 
generate  artificial  turbulent  surface  heat  fluxes,  again  resulting  in  a  rapid  jump  in 
the  surface  temperature,  this  time  toward  the  ABL  temperature. 

The  only  way  to  prevent  an  unwanted  jump  in  surface  temperature  at  the 
start  of  the  model  run  is  to  specify  radiation  and  snow,  surface  and  ABL 
temperatures  that  are  in  equilibrium.  Three  methods  to  accomplish  this  are: 

(1)  alter  the  snow/ice  model, 

(2)  alter  the  radiation  model  or 

(3)  alter  the  low-level  measured  temperature  profile  used 
for  the  initialization 

The  equilibrium  surface  temperature  is  not  greatly  sensitive  to  changes 
in  snow  depth  or  characteristics,  as  shown  in  Chapter  IV.  Large,  unrealistic 
alterations  to  the  snow  characteristics  would  be  required  to  make  the  snow/ice 
model  produce  the  measured  surface  temperature,  eliminating  choice  (1). 

Further,  it  is  desirable  to  preserve  the  measured  thermal  structure  of  the 
lower  atmosphere  as  much  as  possible  since  this  is  expected  to  be  an  important  link 
between  clouds  and  wind  stress.  Therefore  (3)  is  eliminated  to  have  the  thermal 
structure  altered  by  modeled  cloud  and  ABL  physics,  not  artificially. 

This  leaves  (2)  altering  the  radiation  model,  as  the  best  choice.  It  is 
realistic  because  the  snow/ice  model  is  quite  sensitive  to  changes  in  longwave 
radiation,  so  that  drastic  changes  are  not  required  to  produce  the  measured  surface 
temperature.   Also  it  was  shown  earlier  how  the  radiation  model  has  considerable 
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B.       OCT  22  CASE  STUDY  -    CLEARING  EVENT 

Three  case  studies  were  chosen  to  illustrate  cloud  effects  on  wind  stress.  The 
first  two  cases  represent  cloud  forming  and  cloud  clearing,  respectively,  above  the 
ABL.  The  third  case  study  involves  cloud  formation  and  dissipation  within  the  ABL. 

The  first  case  study  at  the  CEAREX  drift  location  used  a  rawinsonde  profile 
at  2300  October  21  to  initialize  the  model  at  t  =  0  hours  (Figure  40).  The  case 
study  ends  at  1500  October  22,  model  t  =  16  hours.  All  the  case  studies  presented 
here  are  based  on  16  hour  model  runs.  Some  moisture  at  cloud  level  was  artificially 
removed  from  the  model  at  t  =  2  hours  in  order  to  simulate  the  observed  cloud 
dissipation  which  occurred  during  the  first  three  hours  of  this  case  study  period. 

During  this  period,  overcast  conditions  changed  to  clear  and  the  surface 
temperature  underwent  a  sharp  decrease  (Figure  41b).  The  measured  radiation 
value  of  almost  270  Wm~2  at  t  =  0  dropped  to  about  160  Wm'2  by  t  =  4  hours  as  the 
clouds  dissipated  (Figure  41a).  Profile  measurements  of  surface  heat  flux  (not 
shown)  averaged  -20  Wm"",  indicating  a  stable  Type  R  event.  This  was  one  of  the 
larger  magnitude  heat  flux  events  measured  with  the  profile  tower  during  the 
CEAREX  drift.  This  situation  is  similar  to  the  standard  case  modeled  in  Chapter 
IV  except  that  the  surface  wind  speed  was  a  little  higher,  about  6  ms"  instead  of  5 
ms"1  and  h  adjusted  to  124  m  instead  of  the  prescribed  50  m  in  the  standard  case. 

The  model  surface  heat  fluxes  were  -30  Wm"2  to  -50  Wm"2.  Despite  the  larger 
predicted  heat  loss  compared  to  the  measurements,  the  model  air  temperature 
decreased  at  only  about  one-half  the  rate  of  the  measurements,  although  the  shape 
of  the  temperature  curves  in  Figure  41  are  similar,  resembling  a  negative 
exponential  decrease  with  time  (e"  ).  The  measurements  indicated  cold  advection 
which  magnified  the  temperature  decrease  and  somewhat  negated  the  large  heat 


123 


In  the  final  hours  of  the  snow/ice  model  initialization,  the  value  of  the 
atmosphere-snow  total  heat  flux  (net  radiation  and  turbulent  fluxes)  was  chosen  by 
trial  and  error  which  would  result  in  the  exact  surface  temperature  which  was 
measured  at  the  start  of  the  case  study  period.  This  prevents  the  type  of  energy 
imbalance  at  the  start  of  the  model  run  described  above. 

The  ABL  model  performs  a  half-hour  turbulence  spin-up  by  keeping  the 
mean  prognostic  variables  fixed  while  allowing  the  second  order  variables  and  the 
TKE  to  adjust  from  an  initial  value  of  zero,  except  the  TKE,  which  starts  at  a 
neutral  atmosphere  value  (Burk  and  Thompson,  1989).  Time  equals  zero  after  this 
spin-up. 

The  recommended  method  for  initializing  the  mean  wind  velocity 
components  is  to  set  the  true  u  and  v  wind  components  to  the  geostrophic  value  at 
all  levels  and  let  the  model  spin-down  the  winds  after  the  start  of  the  model  run 
(Burk,  personal  communication).  This  prevents  inertial  oscillations  and  the 
adjustment  usually  requires  only  two  or  three  hours.  The  problem  with  this  method 
is  that  at  the  start  of  the  model  run,  the  surface  wind  speed  is  too  high  and  may  alter 
the  structure  of  the  ABL  before  adjustment  occurs.  This  problem  is  more 
pronounced  during  common  low  Arctic  inversion  situations. 

For  the  case  studies,  a  preliminary  run  was  used  to  provide  the  wind 
profiles  for  the  start  of  the  actual  case  study  run  and  adjusted  by  trial  and  error  to 
prevent  inertial  oscillations.  These  were  compared  with  spin-down  versions  of  the 
same  cases.  Generally,  the  cloud  effects  were  similar  using  either  method  so  this 
was  not  a  crucial  issue  for  studying  the  effect  of  clouds  on  wind  stress. 
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B.   OCT  22  CASE  STUDY  -  CLEARING  EVENT 

Three  case  studies  were  chosen  to  illustrate  cloud  effects  on  wind  stress.  The 
first  two  cases  represent  cloud  forming  and  cloud  clearing,  respectively,  above  the 
ABL.  The  third  case  study  involves  cloud  formation  and  dissipation  within  the  ABL. 

The  first  case  study  at  the  CEAREX  drift  location  used  a  rawinsonde  profile 
at  2300  October  21  to  initialize  the  model  at  t  =  0  hours  (Figure  40).  The  case 
study  ends  at  1500  October  22,  model  t  =  16  hours.  All  the  case  studies  presented 
here  are  based  on  16  hour  model  runs.  Some  moisture  at  cloud  level  was  artificially 
removed  from  the  model  at  t  =  2  hours  in  order  to  simulate  the  observed  cloud 
dissipation  which  occurred  during  the  first  three  hours  of  this  case  study  period. 

During  this  period,  overcast  conditions  changed  to  clear  and  the  surface 
temperature  underwent  a  sharp  decrease  (Figure  41b).  The  measured  radiation 
value  of  almost  270  Wm  at  t  =  0  dropped  to  about  160  Wm  by  t  =  4  hours  as  the 
clouds  dissipated  (Figure  41a).  Profile  measurements  of  surface  heat  flux  (not 
shown)  averaged  -20  Wm  ",  indicating  a  stable  Type  R  event.  This  was  one  of  the 
larger  magnitude  heat  flux  events  measured  with  the  profile  tower  during  the 
CEAREX  drift.  This  situation  is  similar  to  the  standard  case  modeled  in  Chapter 
IV  except  that  the  surface  wind  speed  was  a  little  higher,  about  6  ms"  instead  of  5 
ms"1  and  h  adjusted  to  124  m  instead  of  the  prescribed  50  m  in  the  standard  case. 

The  model  surface  heat  fluxes  were  -30  Wm""  to  -50  Wm"'.  Despite  the  larger 
predicted  heat  loss  compared  to  the  measurements,  the  model  air  temperature 
decreased  at  only  about  one-half  the  rate  of  the  measurements,  although  the  shape 
of  the  temperature  curves  in  Figure  41  are  similar,  resembling  a  negative 
exponential  decrease  with  time  (e"  ).  The  measurements  indicated  cold  advection 
which  magnified  the  temperature  decrease  and  somewhat  negated  the  large  heat 
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Figure  40  Model  (solid  line)  and  Measured  (dashed  line)  Profiles  of  Potential 
Temperature  for  the  October  21-22  Case  Study.  The  lines  labelled  "1"  are  at  model 
t  =  0  hours  and  the  lines  labeled  "2"  are  at  I  =   12  hours. 
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Figure  42  Model  Output  and  Measured  Values  of  Various  Parameters  Related  to 
Wind  Stress  for  the  October  22  Case  Study.  The  top  panel  (a)  shows  the  model 
output  for  the  normalized  quadratic  reduction  factor  (solid  line),  the  normalized 
drag  coefficient  (dotted  line)  and  the  normalized  quadratic  geostrophic  drag 
coefficient  (dashed  line).  The  middle  panel  (b)  also  shows  the  model  output 
normalized  quadratic  reduction  factor  (solid  line)  and  the  measured  values  (x's). 
Both  panels  have  a  straight  solid  line  to  use  as  a  reference  for  the  average  value  of 
1.0.  The  bottom  panel  (c)  shows  model  output  (solid  line)  and  measured  (x's) 
turning  angle.  See  text  for  definitions  of  these  parameters. 
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flux  into  the  surface.  The  cold  advection  extended  to  at  least  2000  m,  as  indicated 
by  the  measured  line  2  in  Figure  40,  which  has  clearly  shifted  to  the  left  of  line  1  at 
all  levels.  The  presence  of  cold  advection  can  be  explained  by  the  air  within  the 
ABL,  after  the  clearing  event,  being  exposed  to  clear  skies  and  cooling  surface 
conditions  for  a  much  longer  time  than  suggested  by  the  local  conditions. 

Despite  not  accounting  for  cold  advection,  the  model  does  a  fairly  good  job  of 
reproducing  the  internal  boundary  layer  which  forms  after  the  clearing  event.  After 
twelve  hours,  both  the  model  and  measured  ABL  heights  were  about  100  m  (Figures 
40  and  41).  The  measured  ABL  height  shown  in  Figure  40  was  based  on  the  height 
of  the  lowest  temperature  inversion.  The  model  ABL  height  shown  in  Figure  41  was 
based  on  the  level  where  the  TKE  reaches  10%  of  its  near-surface  value. 

Measured  and  modeled  parameters  relating  to  wind  stress  are  shown  in  Figure 

42.    One  parameter  plotted  in  Figure  42a  is  the  normalized  quadratic  wind  speed 

* 
ratio,  CUG   (solid  line).  This  is  defined  as 

M 

QjG      ~      O  "347         •  (    °' 


u10 


The  0.347  normalization  factor  is  the  average  value  of  \Tj~  J  measured  during  the 
CEAREX  drift.  The  value  of  the  geostrophic  wind  speed,  U  ,  was  determined  from 

an  array  of  pressure  buoys  which  surrounded  the  Polarbjoern  during  the  drift  phase. 
The  geostrophic  wind  was  estimated  to  be  accurate  to  within  1.5  ms-1.  This 
parameter  represents  a  non-dimensional  outer  layer  effect  or  wind  speed  effect. 
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Figure  42  Model  Output  and  Measured  Values  of  Various  Parameters  Related  to 
Wind  Stress  for  the  October  22  Case  Study.  The  top  panel  (a)  shows  the  model 
output  for  the  normalized  quadratic  reduction  factor  (solid  line),  the  normalized 
drag  coefficient  (dotted  line)  and  the  normalized  quadratic  geostrophic  drag 
coefficient  (dashed  line).  The  middle  panel  (b)  also  shows  the  model  output 
normalized  quadratic  reduction  factor  (solid  line)  and  the  measured  values  (x's). 
Both  panels  have  a  straight  solid  line  to  use  as  a  reference  for  the  average  value  of 
1.0.  The  bottom  panel  (c)  shows  model  output  (solid  line)  and  measured  (x's) 
turning  angle.  See  text  for  definitions  of  these  parameters. 
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Another  value  in  Figure  42a  is  the  drag  coefficient  ratio  or  normalized  surface 
drag  coefficient,  Cdl0  (small  dashes)  defined  as 


u,^2 


*  ft  'Ul(T  Cd10 

C<no     =      2.3  xlO"3  Cdnl0      •  (29) 

This  is  the  same  parameter  used  in  Chapter  IV  to  describe  surface  layer  stability 

effect. 

The  final  parameter  in  Figure  42a  is  the  normalized  quadratic  geostrophic 
* 
drag  coefficient,  CG  (dashed  line)  defined  as 


u 


c°*  "  ttrf  •  (30) 


_4 
where  8.2  x  10     is  the  average  measured  value  during  the  CEAREX  drift.    This 

parameter  represents  a  non-dimensional  wind  stress. 

Because 

C       -    C  C  \3i) 

^G       _     ^dlO       MjG         ' 


the  relation  between  wind  stress  and  the  geostrophic  forcing  is  a  product  of  the 

*  * 

surface  layer  effect,  Cdl0  ,  and  an  outer  layer  effect  CUG  .    This  is  the  normalized 

version  of  Equation  (6). 

Although  surface  stress  measurements  were  not  usually  available,  there  were 
always  continuous  measurements  of  the  wind  speed  ratio,  CUG  .  They  appear  in 
Figure  42b  along  with  a  repeat  of  the  model  CUG  value  plotted  in  Figure  42a. 
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to  keep  up  and  a  downward  turbulent  heat  flux  is  generated,  resulting  in  a  strong 
surface  stability  effect. 

During  the  third  regime  (t  >  5  hours)  the  wind  stress  changes  slightly  but 
remains  significantly  less  than  during  regime  1.  The  surface  temperature  begins  to 
drop  at  a  slower  rate  allowing  the  10  meter  and  ABL  air  temperature  to  catch  up 
and  form  a  less  stable  internal  layer  below  the  new  capping  inversion.  When  this 
occurs,  CUG    starts  to  level  off.    The  mixed  layer  and  surface  layer  continue  to 

become  less  stable,  but  the  entire  surface/ABL  system  is  cooling  faster  than  the  free 
atmosphere  (Figure  40).  The  effect  on  wind  stress  is  that  during  regime  3,  Cdl0  will 
become  closer  to  one  (neutral)  but  CUG  and  hence  CG  will  remain  "permanently" 
lower  than  average.  This  is  because  the  ABL  during  regime  3  is  shallower  than 
during  regime  1  and  the  continuously  strengthening  capping  inversion  limits  further 
growth  or  contact  with  higher  momentum  air  in  the  free  atmosphere. 

The  measured  geostrophic  wind  speed  and  direction  were  based  on  one  hour 

averages.  The  resulting  measured  CVG   and  a  for  each  hour  are  plotted  on  Figures 

* 
42b  and  42c.    The  measured  CUG    does  not  agree  well  with  the  model  results. 

* 
Although  there  is  a  drop  in  measured  CL,G    of  about  40%  two  hours  after  the 

clearing  event  (t  =  4  hours),  there  is  also  an  unexpected  increase  at  t  =  2  hours  and 
after  t  =  13  hours.  All  the  modeled  CVG  values  were  higher  than  the  observations 

after  the  clearing  event. 

The  measured  a  is  also  poorly  correlated  with  the  model  values.  The 
observations  show  a  significant  trend  toward  greater  turning  angle,  which  might  be 
expected  as  the  internal  stable  layer  strengthens.  However  the  model  predicts  only 
small  changes  in  a. 
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event,  and  that  the  wind  stress  would  be  30%  less  at  t  =  tmax  due  to  the  surface  layer 

stability  effect. 

Wind  speed,  U10,  was  kept  constant  for  the  Chapter  IV  standard  case.    But 
here  it  is  shown  that  changes  in  U10  due  to  cloud  events  are  as  important  as  surface 

stability  effects,  thus  the  need  for  complete  ABL  physics  model  studies.     The 
combined  effect  of  Cdl0  and  CUG  ,  as  represented  by  CG   (dashed  line)  is  to  reduce 

the  modeled  wind  stress  by  50%  one  hour  after  the  clearing  event  and  20%  to  30% 

for  the  rest  of  the  run.  If  the  model  is  run  without  removing  the  clouds  (not  shown) 

* 
the  value  of  CG    maintains  a  constant  value  of  about  1.05;  therefore  the  clearing 

event  has  had  a  large  effect  on  wind  stress,  according  to  the  model. 

*  *  * 

The  time  series  of  model  CG  ,  Cdl0    and  CUG    (Figure  42a)  reveals  three 

"regimes"  which  illustrate  the  effect  of  boundary  layer  structure  and  cloud  clearing 

on  wind  stress.    In  the  first  regime  (0  <  t  <  2)  there  is  a  a  slightly  stable  surface 

layer  below  a  mixed  layer  capped  by  a  cloudy  inversion  at  300  m  (Figure  40,  line  1). 

The  lack  of  heat  fluxes  or  unusual  ABL  structure  results  in  average  wind  stress 

* 
values,  i.e.  CG   is  about  one.  (Here  t  refers  to  the  time  of  model  run,  not  time  after 

cloud  cover  change  as  in  Chapter  IV.) 

In  the  next  regime  (2  <  t  <  5),  the  wind  stress  is  undergoing  rapid  changes  in 
response  to  the  clearing  event.  After  a  very  rapid  initial  decrease  before  t  =  tmax  the 

wind  stress  more  gradually  returns  toward  a  value  of  about  75%  of  the  average 

value  at  about  t  =  6  hours.  During  this  time  period  a  strong  surface-based  inversion 

has  formed  (not  shown).    The  effect  of  the  surface -based  inversion  is  to  prevent 

* 
momentum  from  reaching  the  surface  layer  so  that  CUG   is  reduced.    At  the  same 

time  the  surface  temperature  is  dropping  so  fast  that  the  air  temperature  is  not  able 
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to  keep  up  and  a  downward  turbulent  heat  flux  is  generated,  resulting  in  a  strong 
surface  stability  effect. 

During  the  third  regime  (t  >  5  hours)  the  wind  stress  changes  slightly  but 
remains  significantly  less  than  during  regime  1.  The  surface  temperature  begins  to 
drop  at  a  slower  rate  allowing  the  10  meter  and  ABL  air  temperature  to  catch  up 
and  form  a  less  stable  internal  layer  below  the  new  capping  inversion.  When  this 
occurs,  CUG    starts  to  level  off.    The  mixed  layer  and  surface  layer  continue  to 

become  less  stable,  but  the  entire  surface/ABL  system  is  cooling  faster  than  the  free 
atmosphere  (Figure  40).  The  effect  on  wind  stress  is  that  during  regime  3,  Cdl0  will 
become  closer  to  one  (neutral)  but  CUG  and  hence  CG  will  remain  "permanently" 
lower  than  average.  This  is  because  the  ABL  during  regime  3  is  shallower  than 
during  regime  1  and  the  continuously  strengthening  capping  inversion  limits  further 
growth  or  contact  with  higher  momentum  air  in  the  free  atmosphere. 

The  measured  geostrophic  wind  speed  and  direction  were  based  on  one  hour 

* 
averages.  The  resulting  measured  CVG   and  a  for  each  hour  are  plotted  on  Figures 

* 
42b  and  42c.    The  measured  CUG    does  not  agree  well  with  the  model  results. 

Although  there  is  a  drop  in  measured  CUG    of  about  40%  two  hours  after  the 

clearing  event  (t  =  4  hours),  there  is  also  an  unexpected  increase  at  t  =  2  hours  and 
after  t  =  13  hours.  All  the  modeled  CVG  values  were  higher  than  the  observations 

after  the  clearing  event. 

The  measured  a  is  also  poorly  correlated  with  the  model  values.  The 
observations  show  a  significant  trend  toward  greater  turning  angle,  which  might  be 
expected  as  the  internal  stable  layer  strengthens.  However  the  model  predicts  only 
small  changes  in  a. 
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In  summary,  in  this  case  study  the  model  qualitatively  simulated  the 
thermodynamics  of  a  Type  R  event  well.  The  actual  cooling  was  greater  because  the 
model  did  not  include  the  cold  advection  which  occurred  in  nature,  but  the  general 
pattern  of  the  surface  temperature  decrease,  negative  surface  heat  fluxes  and 
internal  boundary  layer  formation  was  comparable.  The  model  generated  an  ABL 
depth,  h,  similar  to  the  measured  value  after  12  hours.  The  model  predicted  a  major 
(50%)  short  term  reduction  in  wind  stress  and  a  longer  term  significant  (20%) 
reduction.  The  short  term  model  results  were  consistent  with  the  results  of  the 
simple  thermodynamics  approach  used  in  Chapter  IV.  However,  the  measured  wind 
speed  ratio  and  turning  angle  showed  some  large  variations  which  may  have  been 
related  to  the  clearing  event,  but  they  are  not  reproduced  well  by  the  model. 

C.       OCT  30  CASE  STUDY  -  CLOUD  FORMS  ABOVE  ABL 

This  case  study  investigates  the  response  of  the  ABL  and  wind  stress  to  cloud 
formation  just  above  the  ABL.  The  model  is  initialized  (model  t  =  0)  with  the  2300 
October  29  rawinsonde  profile  (Figure  43,  line  1).  The  measured  geostrophic  wind 
speed  was  between  9  and  10  ms"  throughout  the  period.  "Very  clear"  conditions 
were  observed  six  hours  prior  to  the  study  period  and  for  the  first  four  hours  of  the 
case  study.  At  t  =  5  hours  mid-level  clouds  covered  the  sky.  By  t  =  7  hours  the 
clouds  had  lowered  to  just  above  the  ABL.  Associated  with  the  cloud  lowering  was 
an  increase  in  the  measured  downward  radiation  (Figure  44a).  At  t  =  14.5  hours 
the  clouds  were  observed  to  clear,  resulting  in  lower  radiation.  The  clearing  event 
at  this  time  was  not  modeled. 

The  model  created  clouds  in  a  layer  650  to  1000  m  above  the  surface  at  t  =  7 
hours,  causing  a  more  sudden  jump  in  radiation  than  measured,  but  at  about  the 
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Figure  43  Potential  Tempera  lure  Profiles  for  the  October  30  Case  Study.  The 
model  (solid  line)  and  measured  (dashed  line)  initial,  1 2  hour  and  16  hour  profiles 
are  labeled  "l",  "2",  and  "3". 
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Figure  44  Thermodynamic  Variable  Times  Scries  lor  October  30  Case  Study.  Same 
point  label  and  line  definitions  ;is  Figure  4  I,  with  some  changes  in  the  temperature 
and  radiation  Y-axis  scales. 
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correct  time  (Figure  44a).  The  measured  air  temperature  increased  by  8  C 
compared  to  the  model  2  C  during  the  eight  hour  period  following  the  cloud 
formation  (Figure  44b).  Toward  the  end  of  the  run,  the  measured  air  temperature 
decreased  3  C,  in  association  with  the  clearing  event.  The  model  did  not  have 
clearing  and  continued  to  increase  the  air  temperature.  The  measured  ABL  at  t  = 
12  hours  was  twice  as  deep  as  the  model  simulation  (Figure  44b)  while  another 
sounding  after  the  clearing  at  t  =  16  hours  was  just  slightly  higher  than  the  model 
value. 

The  measured  potential  temperature  profiles  show  a  warming  at  all  levels 
below  1400  m  during  the  first  12  hours  (Figure  43).  The  final  (t  =  16  hours) 
measured  profile  shows  that  a  level  of  cooler  air  below  250  m  has  appeared.  Most 
of  these  temperature  changes  must  be  due  to  advection,  since  the  surface  and 
radiative  heating/cooling  are  not  great  enough  to  explain  the  measured  temperature 
deviations.  The  model  profiles,  which  have  no  advection,  show  the  formation  of  a 
mixed  layer  in  the  lower  200  m  and  some  cooling  of  this  air.  Similar  to  the  previous 
case  study  period,  the  radiational  effect  of  the  sky  condition  change  on  air 
temperature  is  magnified  by  advection.  This  is  a  Type  R  event  in  the  general  sense. 

The  model  normalized  wind  stress,  CG  ,  is  below  average  during  the  clear 

period  before  t  =  6  hours  (Figure  45a).  When  the  cloud  forms,  wind  stress 
increases  by  almost  50%  in  the  one  hour  and  then  decreases  somewhat.  Later  at  t 
=  13  hours  the  wind  stress  increases  due  to  an  increase  in  the  ABL  depth.  Most  of 
the  changes  in  wind  stress  are  due  to  changes  in  wind  speed  (CUG  )  not  changes  in 
surface  drag  coefficient  (CU10  ). 

The  measured  wind  speed  ratio  was  near  average  for  the  first  nine  hours 
(Figure  45b).    Between  t  =  9  hours  and  t  =  13  hours  the  measured  wind  speed 
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Figure  45  Dynamic  Variable  Time  Series  for  October  30  Case  Study.  Same  point 
label  and  line  definitions  as  Figure  42.  Note  that  unlike  Figure  42,  the  range  of  the 
middle  plot  (b)  is  larger  the  range  of  the  top  plot  (a). 
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almost  doubled  while  the  geostrophic  wind  remained  constant,  causing  CU10  to 
increase  to  2.5  times  the  average  value.  The  measured  wind  speed  changes  had  a 
similar,  but  greatly  exaggerated  pattern  compared  to  the  model  output. 

The  measured  turning  angle,  a,  sharply  decreased  after  the  clouds  formed 
(Figure  45c).  The  model  predicted  a  slight  decrease  in  a  at  this  time  (t  =  6  hours). 
The  model  a  was  much  more  constant  throughout  the  run  than  the  measured  values. 

In  summary,  this  case  study  illustrates  how  upper  level  clouds  can  have  a  large 
effect  on  wind  stress  by  destroying  the  surface  inversion  or  stable  layer  which  is 
common  during  clear  skies  in  the  Arctic  winter.  During  these  moderate  (5  -  10 
ms  )  winds  the  surface  layer  stability  effect  was  small,  but  the  measured  and 
modeled  wind  stress  increased  by  50%  and  250%  respectively  due  to  an  increase  in 
wind  speed  after  cloud  formation.  The  large  increase  in  measured  wind  speed  six 
hour  after  the  initial  cloud  formation  (Figure  45b)  may  have  occurred  when  a  mixed 
layer  underneath  the  lowering  cloud  layer  first  reached  the  ground.  A  cloud 
lowering  process  was  not  modeled.  The  measured  profile  shows  a  complicated 
situation  with  variable  advection  at  different  levels  and  large  changes  in  wind  stress 
that  are  not  entirely  understandable. 

D.      OCT  21  CASE  STUDY  -  CLOUD  FORMATION  AND  DISSIPATION  IN 
THEABL 

This  case  study  is  initialized  with  the  1100  21  October  rawinsonde  sounding. 

The  sounding  revealed  a  surface  inversion  (Figure  46)  with  surface  geostrophic 

winds  of  12.5  ms"1  but  measured  winds  of  around  6  to  7  ms"1  above  and  below  the 

ABL.    The  wind  just  above  the  ABL  is  assumed  to  be  geostrophic  so  there  is  a 

strong  gradient  in  the  geostrophic  wind  which  could  be  caused  by  a  sloping 

inversion.  This  opposing  thermal  wind  is  included  in  the  model  for  the  entire  run. 
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Figure  46    Potential  Temperature  Profiles  for  the  October  21  Case  Study.    Same 
label  and  line  definitions  as  Figure  40. 
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Skies  were  clear  until  t  =  5  hours,  at  which  time  the  measured  and  modeled 
radiation  increased  sharply  as  cloud  formed  at  the  top  of  ABL  (Figure  47a).  At  t  = 
11  hours  the  clouds  were  observed  to  dissipate.  Associated  with  the  cloud  formation 
and  dissipation  was  an  increase  and  decrease  of  measured  temperature  of  6  C 
(Figure  47b).  The  model  temperature  has  the  same  pattern,  but  with  one  third  of 
the  amplitude.  The  measured  ABL  height  after  12  hours  was  twice  the  model 
height  (Figure  47c).  But  the  soundings  (Figure  46)  reveal  a  similar  pattern  in  the 
temperature  structure  between  the  measurements  and  the  model.  Both  have  an 
upper  level  mixed  layer  capped  by  an  inversion  associated  with  the  cloud  layer. 
Below  the  mixed  layer  the  rawinsonde  profile  has  a  stable  layer  while  the  model  had 
another  inversion/mixed  layer  combination.  In  other  similar  case  studies  not  shown 
here,  any  clouds  which  formed  in  the  ABL  tended  to  separate  from  the  ABL  and 
increase  in  height.  The  clouds  were  kept  low  in  this  case  study  by  imposing  a 
negative  moisture  tendency  just  above  the  ABL.  Still,  layer  separation  occurred. 

The  model  wind  stress  was  lower  than  average  during  the  initial  five-hour 
clear  period,  due  equally  to  low  normalized  wind  speed  and  drag  coefficient  (Figure 
48a).  Immediately  after  cloud  formation  the  wind  stress  increased  by  about  15%, 
mostly  because  of  the  surface  layer  stability  effect.  The  wind  speed,  as  shown  by 
CUG  ,  does  not  increase  by  much  during  the  first  two  hours  after  cloud  formation  (5 

<  t  <  7  hours)  because  the  upper  level  air  that  is  being  entrained  does  not  have  the 
usual  high  momentum  because  of  the  opposing  thermal  wind.  After  t  =  8  hours  the 
surface  layer  stability  is  nearly  neutral  but  the  wind  speed  is  enhanced  causing  the 
wind  stress  to  become  20%  greater  than  normal.  At  t  =  13  hours  the  model  clouds 
are  removed  and  wind  stress  drops  by  40%  due  to  both  wind  speed  (outer  layer)  and 
surface  layer  effects. 
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Figure  47  Thermodynamic  Variable  Times  Series  lor  October  21  Case  Study.  Same 
point  label  and  line  definitions  as  Figure  41,  with  some  changes  in  the  temperature 
and  radiation  Y-axis  scales. 
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Figure  48    Dynamic  Variable  Time  Series  for  October  21  Case  Study.    Same  point 
label  and  line  definitions  as  Figure  42. 
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The  measured  CUG  was  lower  than  normal  before  the  cloud  formation 
(Figure  48b)  and  increased  by  about  50%  and  dropped  at  the  end  of  the  run.  This 
pattern  matched  the  model  quite  well.  The  modeled  and  measured  a  were  quite 
high,  as  expected,  during  the  five  hour  initial  clear  period  (Figure  48c).  After  the 
clouds  formed,  the  measured  a  decreased  to  a  very  small  value  while  the  model  a 
decreased  only  a  little. 

In  summary,  the  surface  layer  stability  in  this  case  study  followed  the  same 
pattern  as  the  previous  case  studies  and  Chapter  IV  cases,  with  a  rapid  deviation 
immediately  after  the  sky  condition  change,  and  then  a  slow  decrease  toward 
neutral.  The  surface  layer  stability  effect  was  not  great  due  to  the  relatively  high 
winds.  The  major  factor  causing  the  modeled  changes  in  wind  stress  was  the  wind 
speed  variation.  Unlike  the  earlier  cases,  the  modeled  wind  speed  was  not  greatly 
affected  by  the  cloud  formation  until  a  few  hours  later.  The  model  appeared  to 
reproduce  the  wind  stress  effects  of  the  cloud  formation  and  clearing  quite  well. 
Here,  the  effect  of  clouds  was  significant  even  during  higher  wind  speeds. 

E.       OTHER  MODEL  SIMULATION  RESULTS 

Many  other  situations  related  to  cloud  effects  were  simulated  with  the 
ABL/Snow/Ice  model.  Some  general  results  from  these  simulations  will  be 
summarized  without  showing  each  case.  The  response  of  wind  stress  to  cloud  effects 
depended  on  the  dynamic  and  thermodynamic  structure  of  the  atmosphere,  as  well 
as  the  ice/snow  characteristics  discussed  in  Chapter  IV. 

An  important  dynamic  feature  is  the  thermal  wind.  The  model  showed  that 
presence  of  a  thermal  wind  can  counteract  or  amplify  cloud  effects  depending  on 
the  orientation  of  the  thermal  wind  and  whether  the  ABL  is  growing  or  shrinking. 


142 


Simulations  of  strong  surface-based  inversions,  which  often  form  under  clear 
skies  during  light  winds,  resulted  in  the  lowest  normalized  wind  stress  values.  When 
clouds  form  during  these  conditions,  the  relative  change  in  wind  stress  is  large 
because  the  normalized  wind  stress  changes  from  a  very  low  value  to  an  average  or 
above  average  value. 

The  largest  effect  on  wind  stress  occurred  when  there  was  a  very  low  stable 
ABL  with  a  deep  neutral  layer  above.  If  cloud  formation  can  warm  the  cold  air  and 
deepen  the  ABL  to  the  top  of  the  neutral  layer,  wind  stress  increases  by  a  factor  of 
two  or  three. 

As  a  test  of  long  term  cloud  effects,  the  model  was  run  for  periods  of  several 
days.  It  was  hard  to  draw  conclusions  with  the  1-D  model  because  the  long  term 
effects  of  the  clouds  were  dependent  on  the  specification  of  temperature  and 
humidity  advection.  Without  any  warm  advection,  a  cloud  layer  near  in  the  lower 
atmosphere  continues  to  cool  to  unrealistically  low  temperatures  or  dissipate  after  a 
few  days.  The  location  of  the  clouds  was  a  crucial  factor  in  determining  the  long 
term  wind  stress  effect.  In  the  long  term,  the  effect  of  the  surface  fluxes  becomes 
small  but  the  depth  of  the  ABL  is  an  important  factor.  If  the  clouds  are  well  above 
the  ABL,  their  effect  on  wind  will  become  small  after  the  initial  surface  flux  event. 
But  if  there  is  a  cloud  layer  just  above  the  ABL  it  may  be  cool  enough  to  be 
incorporated  into  the  ABL,  thus  increasing  h  and  the  wind  stress. 

F.       CASE  STUDY  CONCLUSIONS 

The  numerical  studies  show  that  clouds  can  have  a  substantial  effect  on  wind 
stress  by  influencing  the  surface  layer  wind  speed  and  stratification.  With  moderate 
surface  winds,  a  cloud  formation/dissipation  event  will  change  wind  stress  by  at  least 
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40%  temporarily  and  about  15%  permanently.  The  verification  measurements 
showed  amplified  temperature  responses  in  some  cases,  caused  by  advection.  The 
normalized  verification  wind  stress  had  large  deviations  that  appeared  to  be  related 
to  the  cloud  events  but  did  not  match  the  model  well  in  most  cases.  The  deviations 
in  measured  turning  angles  were  poorly  modeled. 

The  ABL  in  the  Arctic  is  rarely  in  a  steady-state  situation.  Horizontal  changes 
in  upwind  surface  roughness,  cloud  cover,  temperature,  geostrophic  forcing,  etc. 
change  on  time  scales  shorter  than  the  model  runs,  which  assume  constant  forcing 
conditions.  Therefore  verification  of  the  model  with  the  measurements  was  not 
expected. 

The  examination  based  on  measurements  and  model  predictions  shows  that, 
all  other  factors  being  equal,  clouds  have  a  significant  effect  on  wind  stress.  The 
surface  layer  effect  is  quite  straightforward  and  the  complete  physics  model  gives 
similar  results  as  the  Chapter  IV  thermodynamic  model.  The  outer  layer  (wind 
speed  ratio)  effect  is  often  greater  than  the  surface  layer  effect,  especially  during 
higher  winds.  This  effect  depends  on  a  number  of  factors  such  as  thermal  wind, 
inertial  effects,  ABL  depth,  surface  fluxes  and  cloud  location. 

Verification  of  the  cloud  effect  is  difficult  because  the  measured  wind  stress 
has  a  large  amount  of  unaccounted  variability.  The  model  predicts  stress  values 
which  seem  to  be  reasonable  averages  of  several  measurements,  but  may  not 
compare  well  with  one  particular  measurement. 
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VI.  STATISTICAL  RELATIONSHIPS  BETWEEN  CLOUDS  AND  WIND  STRESS 

The  effect  of  clouds  on  wind  stress  was  examined  for  a  few  short  periods  in  the 
previous  chapter.  To  quantify  average  cloud  effects,  the  relationship  between 
clouds  and  wind  stress  using  all  valid  geostrophic  wind  measurements  from  the 
CEAREX  drift  will  be  studied. 

The  relationship  between  ABL  structure  and  clouds  was  quantified  in  Chapter 
III.  This  chapter  continues  that  line  of  reasoning  to  explain  the  effect  of  clouds  on 
wind  stress  by  first  describing  the  relationship  between  wind  stress  and  bulk  ABL 
properties.  Then  the  direct  correlation  between  observed  cloud  cover  or  relative 
humidity  and  wind  stress  will  be  analyzed. 

Because  the  surface  heat  flux  and  momentum  measurements  were  usually  not 

u, 
available,  all  geostrophic  drag  coefficients,  CG  =  77",  in  this  chapter  assume  a 

constant  surface  drag  coefficient,  Cdl0.  Therefore,  variations  of  CG  are  only  due  to 
ABL-influenced  variations  in  the  ratio  of  surface  wind  speed  to  geostrophic  wind 
speed.  The  wind  vector  was  determined  from  buoy  array. 

A.   WIND  STRESS  AND  ABL  STRUCTURE 

Clouds  are  expected  to  affect  wind  stress  by  changing  the  structure  of  the 
ABL,  which  in  turn  affects  the  transfer  of  momentum  to  the  surface.  The 
relationship  between  ABL  structure  and  clouds  was  quantified  in  Chapter  III.  In 
this  section,  the  relation  between  ABL  structure  and  wind  stress  will  be  verified  for 
winter  Arctic  conditions.  A  study  (Guest,  1992)  of  all  factors,  not  just  clouds, 
affecting  CG    will  be  summarized  because  ABL,  including  surface  layer,  structure 
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provides  the  mechanism  by  which  clouds  affect  wind  stress.    ABL  structure  must 
have  an  effect  on  wind  stress  in  order  for  clouds  to  affect  wind  stress. 

Using  the  CEAREX  measurements,  OD  defined  a  static  stability  measure 
based  on  the  temperature  difference,  A6,  and  distance,  z^,  between  the  surface  and 

900  mb  level 

2  g  A# 

N «  =  tZ  ■  (31) 

The  following  empirical  equation  explains  17%  of  the  variation  in  geostrophic 
drag  coefficient,  CG,  assuming  constant  Cdnl0. 


N       \4 

JN900    v^ 


CG   =0.037-  0.0083  \^J       .  (32) 

An  expanded  CEAREX  data  set,  compared  to  OG,  was  created  for  this  study. 
Data  from  the  entire  CEAREX  drift  period  September  17  -  November  17  were  used 
except  for  some  obviously  spiky  geostrophic  winds,  very  low  winds  and  a  few 
outliers.  Many  different  stability  parameters  were  tested  for  correlation  with  CG. 
The  method  based  on  the  N^  parameter  explained  14.3%  of  the  CG  variance  using 
the  expanded  data  set  (Figure  49).  This  was  as  good  or  better  than  any  other  tested 
stability  parameters. 

The  best  predictor  for  the  turning  angle  was  based  on  a  measure  of  the 
strength  of  the  lowest  part  of  the  inversion  defined  as 


Ninverei0n    =Ln(f^)     .  (33) 
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Figure  49  The  Effect  of  Nratio  =  (N^/0.24)4  on  Geostrophic  Drag  Coefficient 
During  the  CEAREX  Drift.  The  solid  line  is  the  linear  regression  fit,  the  inner 
dashed  lines  are  the  95%  confidence  intervals  of  the  regression  and  the  outer 
dashed  lines  are  the  95%  data  range  limits. 
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where  Ad  is  the  difference  in  temperature  between  the  inversion  base  height,  Zi5  and 
100  meters  higher  or  Zj  +  100.  This  parameter  accounts  for  23.3%  of  the  variance 
in  turning  angle  (Figure  50).  Stronger  inversions  are  associated  with  greater  turning 
angles. 

Inversion  base  height,  Zi?  is  also  correlated  with  the  wind  stress  vector  because 
higher  inversion  bases  are  associated  with  increased  CG  (Figure  51)  and  decreased 
turning  angle  (Figure  52).  The  largest  CG  and  turning  angle  effects  are  seen  for  the 
surface-based  inversions,  Zj  =  0,  cases. 

A  statistical  test  to  prove  that  Zj  is  related  to  wind  stress  is  performed  because 

of  the  large  scatter  in  Figures  51  and  52.  The  ABL  data  with  classic  inversions  are 
divided  into  two  groups,  one  group  with  surface  based  inversions  (Zi  =  0  m)  and 
the  other  group  with  elevated  inversions  (Z,  >  100  m).  The  surface  inversion  group 
has  an  average  CG  and  turning  angle  of  0.0237  and  27.8  degrees,  while  the  elevated 
inversion  group  has  values  of  0.00313  and  21.3  degrees,  representing  a  40%  average 
difference  in  scalar  wind  stress.  One-sided  t-tests  showed  that  the  two  groups  are 
different  in  CG  and  turning  angle  with  greater  than  95%  confidence.    While  the 

elevated  inversion  cases  occur  during  clear  and  cloudy  weather,  almost  all  of  the 
surface  inversion  cases  occur  during  clear  weather.  When  a  surface-based  inversion 
is  present,  atmospheric  stratification  has  the  greatest  effect  on  the  measured  CG;  the 

absence  of  clouds  is  required  for  these  situations  to  occur. 

B.      WIND  STRESS  AND  CLOUDS  AND  MOISTURE 

The  previous  section  demonstrated  that  a  statistical  relationship  exists 
between  ABL  structure  and  wind  stress  vector.  This  section  will  look  at  the  direct 
statistical  relationship  between  wind  stress  and  low  cloud  amount,  NL.  For  this,  the 
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Figure  50    The  Effect  of  Nim.c„inn  on   Turning  Angle,  a.    Same  line  definitions  as 
Figure  49. 
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Figure  51  The  Relation  Between  Inversion  Base  Height,  Z>,  find  Geostrophic  Drag 
Coefficient,  C0,  During  the  CRARRX  Drift.  Only  classic  ABLS  (see  Chapter  IV  for 
definition)  are  shown. 
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Figure  52  The  Relation  Between  Inversion  Base  Height,  Zj,  and  Turning  Angle,  a, 
During  the  CEAREX  Drift. 
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CEAREX  data  were  divided  into  a  group  with  no  low  clouds  (NL  =  0  octals)  and  a 
group  with  low  overcast  (NL  =  8  octals). 

The  frequency  histograms  of  CG  for  the  two  groups  reveals  that  there  is  not  a 

clear  separation  in  the  distributions  (Figure  53),  but  that  the  overcast  group  has  a 
higher  average  CG  than  the  clear  group,  a  value  of  0.0311  vs.  0.0280,  respectively, 

representing  a  19%  difference  in  average  normalized  stress.  A  one-sided  t-test 
shows  that  there  is  a  significant  difference  between  the  two  groups  at  the  90% 
confidence  level. 

The  average  turning  angles  for  the  clear  vs.  low  overcast  group  are  25.0  and 
18.7  respectively.  The  difference  is  significant  to  the  95%  confidence  level. 
However,  there  is  considerable  scatter  shown  in  Figure  54;  therefore  virtually  any 
turning  angle  is  possible  for  any  cloud  coverage. 

Another  parameter  related  to  clouds  and  longwave  radiation  is  the  humidity. 
The  relation  between  relative  humidity  at  several  different  levels  and  CG  was  tested. 
The  relative  humidity  at  Z,,  RH(Z,)  as  the  most  highly  correlated  with  wind  stress 
(Figure  55).  An  analysis  of  the  rawinsonde  measurements  of  relative  humidity,  RH, 
showed  that  values  of  greater  than  98%  are  never  measured;  this  is  apparently  the 
limit  of  the  hygrometer  on  the  rawinsondes.  If  RH  was  between  93%  and  98%, 
clouds  may  or  may  not  be  present.  Clouds  are  rarely  observed  if  RH  is  below  90%. 
Therefore,  it  is  difficult  to  be  sure  if  a  cloud  layer  exists  from  a  moist  (RH  >  93%) 
rawinsonde  sounding,  but  a  dry  sounding  always  indicates  no  cloud  is  present.  For 
the  results  of  this  dissertation,  cloud  amount  was  based  on  human  observations,  but 
often  a  rawinsonde  or  model  output  is  the  only  moisture  or  cloud  information 
available.  The  points  to  the  left  of  Figure  55  were  mostly  from  clear  periods  and 
have  generally  low  CG  values.   The  points  to  the  right  were  from  clear  and  cloudy 
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Figure  53  Distribution  of  C(i  for  Clear  (upper)  and  Overcast  (lower,  upside  clown) 
Cases.  The  X-Axis  represents  different  bins  of  C(i  values,  while  the  y-Axis  shows 
the  number  of  occurrences  (frequency). 
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Figure  54    The  Relation  Between  Amount  of  Low  Cloud  Coverage  and  Turning 
Angle. 
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Figure  55  The  Relation  Between  Relative  Humidity  at  the  Height  of  the  Inversion 
Base  and  Geostrophic  Drag  Coefficient. 
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periods  and  have  a  large  range.  Ninety  per  cent  of  the  greatest  CG  values  occurs  at 
RH  >  93%. 

A  similar  relation  between  relative  humidity  and  turning  angle  exists  (Figure 
56).  The  average  turning  angle  is  smaller  for  high  humidity  and  the  lowest  turning 
angles  occur  during  high  humidity  periods. 

C.       CONCLUSIONS  ON  STATISTICAL  RELATIONSHIPS 

These  results  show  a  statistically  significant  relationship  exists  between  ABL 
structure  and  wind  stress,  as  quantified  by  the  stress  parameters  CG  and  a.    The 

static  stability  based  on  temperature  difference  between  the  surface  and  900  mb, 
NgQQ,  from  OD  provided  the  best  fit  for  the  CG  values,  of  the  several  stability 
parameters  tested.  A  function  based  on  the  strength  of  the  lower  inversion,  Ninversion, 
was  the  best  predictor  for  turning  angle.  Surface-based  inversions  had  substantially 
lower  average  normalized  stress  and  greater  turning  angles,  a,  than  elevated 
inversions.     The  variations  in  average  stress  parameters  as  a  function  of  Z;  is 

negligible  when  considering  elevated  inversions  only. 

There  was  a  19%  difference  in  average  normalized  stress  between  clear  and 
overcast  conditions  but  only  6%  of  the  variance  in  CG  is  explained  by  cloud 

coverage.  These  results  represent  the  wind  speed  (outer  layer)  effect  only.  The 
actual  wind  stress  is  also  affected  by  surface  layer  stability  effects.  Because  surface 
layer  stability  effects  usually  augment  the  wind  speed  in  affecting  wind  stress,  as 
shown  in  Chapter  V,  the  average  difference  in  actual  wind  stress  was  larger  than 
20%. 

Correlation  does  not  prove  cause  and  effect;  there  is  the  possibility  that  the 
value  of  CG  and/or  a  is  physically  affecting  the  formation  of  clouds,  instead  of  vice 
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Figure  56  The  Relation  Between  Relative  Humidity  at  the  Height  of  the  Inversion 
Base  and  Turning  Angle. 
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versa,  or  that  some  third  factor  affects  both  clouds  and  wind  stress.  Cloudiness  is  a 
function  of  moisture,  subsidence,  surface  fluxes  and  other  4-dimensional  factors. 
Wind  stress  has  no  effect  on  clouds  above  the  ABL  and  is  only  a  contributing  factor 
for  ABL  clouds,  so  there  is  not  a  strong  physical  link  by  which  changes  in  CG  and/or 

a  can  affect  cloud  formation.  The  geostrophic  wind  speed  is  the  most  dominant 
third  factor  which  causes  a  spurious  correlation  to  exist  between  clouds  and  wind 
stress.  An  increase  in  pressure  gradient  increases  wind  stress  and  deepens  the  ABL, 
making  cloud  formation  more  likely.  But  by  using  a  normalized  wind  stress 
parameter  such  as  CG,  the  wind  speed  effect  on  wind  stress  is  removed.  Other  third 

factors  are  discussed  in  the  final  conclusions. 

Surface-based  inversions  are  associated  with  suppressed  wind  stress,  on  the 
average.  The  surface  base  inversions  virtually  only  occur  during  clear  weather.  But 
clear  weather  does  not  guarantee  the  existence  of  a  surface  inversion.  Therefore, 
although  one  can  predict  wind  stress  is  likely  to  be  suppressed  if  a  surface-based 
inversion  exists,  one  cannot  predict  that  clear  skies  will  always  suppress  wind  stress. 

Another  relationship  was  that  virtually  all  of  the  most  enhanced  CG  values 
occurred  when  RH(Z,)  was  high,  but  high  RH(Z,)  conditions  were  also  associated 
with  normal  CG  values.  So  enhanced  CG  can  be  used  to  predict  high  RH(Zj).  But 
the  reverse  process  is  not  true;  high  RH(Zj)  does  not  guarantee  suppressed  CG.  In 
other  words,  high  RH(Z,)  is  a  necessary  but  not  sufficient  condition  for  large 
enhancement  of  CG  and  low  turning  angle. 

The  enhanced  and  normal  wind  stress  values  which  occur  when  RH(Zj)  is 
almost  100%  can  be  explained  as  follows.    A  larger  CG  difference  is  expected  to 

occur  between  clear  and  cloudy  conditions  than  between  different  RH  values. 
Because  some  of  the  high  RH(Z,)  values  occurred  during  clear  weather,  this  would 
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mix  low  CG  values  from  the  clear  periods  with  the  high  CG  values  from  the  cloudy 
periods.  On  the  other  hand,  the  low  RH  values  were  virtually  all  clear  sky  cases, 
therefore  no  large  enhanced  values  occurred  when  RH  was  below  90%. 

The  best  fit  for  average  CG  as  a  function  of  RH(Zj)  is  not  linear  as  shown  in 
Figure  55.  Most  of  the  change  in  average  CG  occurs  between  93%  to  98%  relative 
humidity.  An  explanation  for  this  could  be  that  much  of  the  variation  apparently 
related  to  RH  is  actually  related  to  the  amount  of  cloud  cases  within  each  RH  bin. 
The  most  changes  in  cloudiness  conditions  between  RH  bins  occurs  above  93%. 

Another  factor  to  consider  is  cloud  history.  Clear  air  which  has  a  high 
humidity  may  have  been  cloudy  recently  and  the  ABL  structure  created  by  the 
previous  cloud  could  still  affect  wind  stress.  Thus,  some  of  the  clear  cases  had  a  CG 

which  had  earlier  been  enhanced  by  clouds.  But  below  90%,  the  clear  cases 
probably  had  been  clear  for  some  time,  so  there  are  no  enhanced  CG  cases. 

Possibly  because  humidity  is  an  indicator  of  cloud  history,  which  may  in  some 
cases  be  more  important  than  current  cloud  conditions,  RH(Zj)  seems  to  be  a  better 
predictor  for  CG  than  current  cloud  conditions. 
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VII.  CONCLUSIONS 

A.       SUMMARY  OF  FINDINGS 

This  study  shows  that  clouds  affect  wind  and  wind  stress  during  the  central 
Arctic  dark  season.  The  significance  of  clouds  is  supported  by  measurements  and 
model  predictions.  However,  a  lower  than  expected  correlation  between  cloud 
coverage  and  wind  stress  exists  when  the  entire  data  set  is  considered. 

Model  predictions  and  observed  results  show  that  clouds  have  a  strong  effect 
on  AEL  structure.  They  also  show  that  ABL  structure  is  linked  to  normalized  wind 
stress  as  represented  by  CG.    A  physical  link  between  clouds  and  wind  stress  has 

been  demonstrated. 

The  effect  of  clouds  on  wind  stress  is  a  function  of  time  after  the  last  cloud 
formation  or  dissipation  event.  Therefore,  existing  cloud  conditions  are  not  as 
related  to  wind  stress  as  the  previous  cloud  changes  or  cloud  history  in  an  air  mass. 
Correlations  between  relative  humidity  at  the  inversion  base,  RH(Zj),  and  wind 

stress  were  higher  than  clouds  vs.  wind  stress  correlations.  This  is  further  evidence 
that  cloud  history  is  important. 

Two  regimes  of  cloud  effects  on  wind  stress  were  identified.  The  first  regime, 
which  lasts  for  a  few  hours,  consists  of  a  sudden  change  in  downward  radiation  at 
the  surface  due  to  a  cloud  formation/dissipation  event  which  changes  the  surface 
temperature  rapidly,  so  that  surface  heat  fluxes  are  generated.  This  causes  surface 
layer  stability  effects  which  change  the  value  of  Cdl0  and  wind  stress.    The  surface 

heat   flux   affects   the   outer  layer   (ABL  above   the   surface   layer)   by   causing 

* 
entrainment  and  modifying  turbulence.  This  changes  surface  wind  speed  and  CUG 
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on  the  same  order  of  magnitude  as  the  Cdl0.  The  combination  of  these  effects  would 
cause  a  wind  stress  change  of  about  40%  maximum.  This  maximum  occurs  about 
one  hour  after  the  cloud  change  for  typical  Arctic  winter  conditions.  The  wind 
stress  changes  during  the  first  regime  depend  primarily  on  wind  speed,  thermal 
wind,  ABL  depth,  magnitude  of  radiation  change,  snow  age  and,  sometimes,  snow 
depth. 

The  second  regime  begins  when  heat  fluxes  become  small,  typically  4  to  24 
hours  after  a  radiation  change.  Surface  heat  fluxes  are  small,  but  the  atmospheric 
thermodynamic  structure  has  undergone  changes  which  permanently  affect  wind 
stress.  A  typical  second  regime  cloud  effect  is  15%.  This  varies  depending  on 
whether  the  cloud  layer  exists  within  or  above  the  ABL.  If  the  cloud  layer  is  within 
the  ABL,  the  cloud  layer  is  dynamically  coupled  to  the  surface  and  the  cloud  effect 
is  greater  than  for  upper  level  clouds.  The  initial  lower  atmospheric  structure  is 
important.  Under  certain  conditions,  the  initial  ABL  structure  will  be  radically 
changed  by  cloud  formation.  For  example,  if  a  weak  surface  inversion  is  capped  by 
a  thick  mixed  layer,  clouds  may  trigger  a  large  deepening  event  and  increase  wind 
stress  more  than  average.  A  cloud  clearing  event  when  the  ABL  is  deep  may  cause 
the  formation  of  an  internal  ABL  and  suppress  wind  stress.  If  the  ABL  is  already 
shallow  before  a  clearing  event  or  already  deep  before  a  clouding  event,  then  the 
second  regime  effects  of  clouds  on  wind  stress  will  be  small. 

The  wind  velocity  at  the  surface  is  forced  not  only  by  surface  geostrophic  wind 
but  by  the  ABL  integrated  geostrophic  wind.  A  thermal  wind  near  the  top  of  the 
ABL,  common  with  sloping  inversions,  will  affect  wind  stress  depending  on  whether 
ABL  depth  is  increasing  or  decreasing  with  time  and  on  the  orientation  and 
magnitude  of  the  thermal  wind,  relative  to  the  wind  stress  vector. 
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B.       SIGNIFICANCE  OF  RESULTS 

Clouds  have  been  shown  to  have  a  significant  but  not  dominant  effect  on  wind 
stress.  The  most  important  factor  affecting  wind  stress  is  the  geostrophic  wind 
vector.  Cloud  effects  are  second  order,  with  about  the  same  order  of  magnitude  as 
surface  roughness  changes  and  thermal  wind  effects.  Changes  in  wind  stress  due  to 
any  of  these  second  order  effects  are  significant,  as  shown  by  models,  but  difficult  to 
detect  from  measurements.  For  example,  model  results  here  and  by  Brown  (1981) 
show  that  including  thermal  wind  effects  should  result  in  changes  in  wind  stress  of 
20%  or  more.  However,  no  thermal  wind  effect  was  found  using  the  CEAREX  data 
(it  was  estimated  from  the  buoy  array,  which  had  temperature  as  well  as  pressure 
sensors).  This  is  an  example  of  a  proven  effect  not  being  readily  evident  from  the 
measurements.  Similarly,  the  cloud  effect  is  large  in  model  predictions  but  less 
obvious  in  the  noisy  statistical  studies. 

The  reason  for  all  the  "noise",  or  apparently  random  variations,  in  wind  stress, 
even  after  accounting  for  geostrophic  wind  vector  and  all  other  known  effects 
(clouds,  thermal  wind,  measurement  error,  etc.)  is  probably  due  to  inertial  effects. 
The  dynamic  forcing  in  the  Arctic  is  continually  changing  and  the  air  parcels  require 
some  time  to  reach  a  force  balance.  All  the  model  runs  assumed  constant  forcing 
and  therefore  did  not  simulate  this  continual  adjustment  process. 

Several  problems  remain  in  understanding  the  relation  between  clouds  and 
air-ice-sea  interactions.  The  model  results  were  usually  based  on  cloud  formation  or 
dissipation  events  which  were  artificially  forced.  The  physics  of  cloud  formation  in 
the  winter  is  not  well  understood.  The  presence  of  ice  crystals  will  continue  to  cause 
problems  with  verification  and  initialization  of  radiation  models  for  Arctic  regions. 
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C.      APPLICATIONS  OF  RESULTS 

Wind  stress  has  special  importance  in  Arctic  regions  because  it  is  the  primary 
driving  force  for  pack  ice  motion.  Ocean/ice  features  such  as  leads,  polynyas,  ice 
edges,  acoustic  noise  and  Ekman  pumping  are  affected  by  variations  in  wind  stress. 
Cloud  conditions  influence  these  variations.  Therefore,  studies  of  these  features 
should  consider  the  effect  of  clouds  on  wind  stress. 

An  objective  is  to  provide  operational  guidance  for  prediction  of  wind  stress  in 
the  Arctic  for  use  as  input  for  dynamic  ice  and  upper  ocean  models.  The  problem 
with  achieving  this  is  that  CG  varies  by  50%.  Current  cloud  state  explains  a  small 
fraction  of  the  measured  variability  in  CG,  and  is  therefore  an  incomplete  wind 

stress  predictor.  Nevertheless,  wind  stress  forecasts  should  be  improved  by  an 
average  of  about  10%  if  CG  and  a  are  based  on  amount  of  low  clouds  rather  than 

kept  fixed.  The  following  values  are  recommended  for  the  central  Arctic  during 
dark  seasons,  for  totally  cloudy  and  clear  situations.  (Partly  cloudy  skies  have 
intermediate  values.) 

Cloudy Clear 

CG  0.0311  0.0280 

a  25.0  degrees  18.7  degrees 

These  values  could  be  used  to  provide  estimates  of  wind  stress  as  input  into  dynamic 
ice  and  ocean  models  when  the  surface  horizontal  pressure  field  and  cloud  amount 
are  known,  but  no  upper  level  information  is  available.  This  is  the  typical  situation 
in  the  Arctic  Ocean,  where  buoys  provide  surface  pressure  information  and  satellites 
reveal  clouds,  but  detailed  information  on  low  level  stratification  is  unavailable. 
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If  stratification  information  is  available,  then  the  methods  suggested  in 
Chapter  V  and  by  OD  provide  somewhat  better  predictions  than  the  above  values. 
Because  of  the  lack  of  routine  detailed  meteorological  measurements  in  the  Arctic 
Ocean,  it  is  doubtful  that  3-dimensional  numerical  ABL  models  can  provide  more 
accurate  wind  stress  predictions  than  these  simple  analytical  relationships. 

This  research  has  shown  that  numerical  models  of  radiation  and  other  physical 
processes  must  be  tuned  to  Arctic  conditions.  Phenomena  unique  to  high  latitudes, 
such  as  low  level  ice  crystals,  cause  effects  which  will  not  be  properly  modeled  if  the 
model  formulations  are  based  on  mid-latitude  verifications. 

Climate  studies  must  consider  the  effect  of  future  changes  in  cloud  coverage 
on  air/ice/sea  interactions.  The  results  show  that  clouds  increase  surface  forcing  in 
marine  Arctic  regions.  This  can  change  ice  coverage  and  may  affect  global  climate. 
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